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Abstract 

The purpose of this review is to provide a comprehensive pedagogical introduction into Keldysh 
technique for interacting out-of-equilibrium fermionic and bosonic systems. The emphasis is 
placed on a functional integral representation of underlying microscopic models. A large part 
of the review is devoted to derivation and applications of the non-linear cr-model for disor- 
dered metals and superconductors. We discuss such topics as transport properties, mesoscopic 
effects, counting statistics, interaction corrections, kinetic equation, etc. The chapter devoted 
to disordered superconductors includes Usadel equation, fluctuation corrections, time-dependent 
Ginzburg-Landau theory, proximity and Josephson effects, etc. (This review is a substantial ex- 
tension of ar Xiv:cond-mat/04 1 2296 1 ) 

Keywords: Keldysh technique; Green functions; kinetic equation; non-linear sigma model; 
mesoscopic systems; fluctuating superconductors. 
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1 Introduction 



1.1 Motivation and outline 

This review is devoted to the Keldysh formahsm for the treatment of out-of-equihbrium interacting 
many-body systems. The technique takes its name from the 1964paper of L. V. Keldysh fT\. Among 
earUer approaches that are closely related to the Keldysh technique, one should mention Konstanti- 
nov and Perel ||2l, Schwinger Q, Kadanoff and Baym H, and Feynman and Vernon jSj. Classi- 
cal counterparts of the Keldysh technique are extremely useful and interesting on their own right. 
These include the Wild diagrammatic technique |6]| and Matrin-Siggia-Rose (MSR) method |7j for 
stochastic systems (see also related work of DeDominicis ID). 

There is a number of presentations of the method in the existing literature 191 [TOl fTTl [121 [131 fT4l 
[Tsl . The emphasis of this review, which is a substantially extended version of Les Houches Session 
LXXX/ lectures IH, is on the functional integration approach. It makes the structure and the internal 
logic of the theory substantially more clear and transparent. We focus on various applications of the 
method, exposing connections to other techniques such as the equilibrium Matsubara method lllTl 
[Tsl and the classical Langevin and Fokker-Planck (FP) equations |[T9l[20l . The major part of the 
review is devoted to a detailed derivation of the non-linear cr-model (NLSM) 1211 [22l [23l [24l . 
which is probably the most powerful calculation technique in the theory of disordered metals and 
superconductors. This part may be considered as complimentary to earlier presentations of the 
replica l26l[27l[28l [2911301 and the supersymmetric OTl [32l l33ll versions of the cr-model. 

Our aim is to expose the following applications and advantages of the Keldysh formulation of 
the many-body theory. 

• Treatment of systems away from thermal equilibrium, either due to the presence of external 
fields, or in a transient regime. 

• An alternative to replica and supersymmetry methods in the theory of systems with quenched 
disorder 

• Calculation of the full counting statistics of a quantum observable, as opposed to its average 
value or correlators. 

• Treatment of equilibrium problems, where Matsubara analytical continuation may prove to be 
cumbersome. 

Our intent is not to cover all applications of the technique that have appeared previously in the 
literature. We rather aim at a systematic and self-contained exposition, helpful for beginners. The 
choice of cited literature is therefore very partial and subjective. It is mainly intended to provide 
more in-depth details about the chosen examples, rather than a comprehensive literature guide. 

The outline of the present review is as follows. We first introduce the essential elements of the 
Keldysh method: concept of the closed contour Sec. II. 21 Green's functions, ext., starting from a sim- 
ple example of non-interacting system of bosons, Sec.|2l and fermions, Sec.jSj Boson interactions, 
the diagrammatic technique and quantum kinetic equation are discussed in Sec. [3] Section |4| is de- 
voted to a particle in contact with a dissipative environment (bath). This example is used to establish 
connections with the classical methods (Langevin, Fokker-Planck, Martin-Siggia-Rose) as well as 
with the equilibrium Matsubara technique. Non-interacting fermions in presence of quenched dis- 
order are treated in Sec. |6] with the help of the Keldysh non-linear cr-model. It is generalized to 
include Coulomb interactions in Sec.|7]and superconducting correlations in Sec. [8] All technicalities 
are accompanied by examples of applications, intended to illustrate various aspects of the method. 
We cover spectral statistics in mesoscopic samples, universal conductance fluctuations (UCFs), shot 
noise and full counting statistics of electron transport, interaction corrections to the transport coef- 
ficients in disordered metals and superconductors. Coulomb drag, etc. We also devote much atten- 
tion to derivations of effective phenomenological models, such as Caldeira-Leggett, time dependent 
Ginzburg-Landau (TDGL), Usadel, etc., from the microscopic Keldysh formalism. 
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1.2 Closed time contour 



Consider a quantum many-body system governed by a (possibly time-dependent) Hamiltonian H(t). 
Let us assume that in the distant past t - -oo the system was in a state, specified by a many-body 
density matrix p(-oo). The precise form of the latter is of no importance. It may be, e.g., the 
equilibrium density matrix associated with the Hamiltonian H{-oo). The density matrix evolves 
according to the Von Neumann equation dip{t) = -i[H{t),p{t)], where we set = 1. It is formally 
solved by p(t) - 'ZV, _oop(-oo)['ZV,,-oo] ' - 'Ut,-<x>p{-°°)'li-oo,t, where the evolution operator is given 
by the time-ordered exponent; 

iA,j. = T exp J-/ J' H{t)A^ = Hm g-'HOftg-'-ffC-^,)* _ . . ^-m'+6,)5,^ 

where an infinitesimal time-step is 6t - {t— t')/N. 

One is usually interested in calculations of expectation value for some observable O (say density 
or current) at a time f, defined as 

<^^^^> - W ^ T^Tr{'//_.,,d'l/._.p(-oo)} , (2) 

where the traces are performed over the many-body Hilbert space. The expression under the last 
trace describes (read from right to left) evolution from t - -oo, where the initial density matrix is 
specified, forward to f, where the observable is calculated, and then backward to t - -co. Such 
forward-backward evolution is avoided in the equilibrium by a specially designed trick. 

Let us recall how it works, for example, in the zero-temperature quantum field theory IfTSl . 
The latter deals with the expectation values of the type (GS|(9|GS) - {0|'Z/_co,r(5^,,_co|0), where 
|GS) = 'i//,_oo|0) is a ground state of full interacting system. The only time dependence allowed for 
the Hamiltonian is an adiabatic switching of interactions on and off in the distant past and future, 
respectively. The evolution operator therefore describes the evolution of a simple non-interacting 
ground state |0) toward |GS) upon adiabatic switching of the interactions. Now comes the trick: one 
inserts the operator 'ZY+oo,-oo in the left-most position to accomplish the evolution along the entire 
time axis. It is then argued that (0|'Z/+oo,-oo - (0|e'^. This argument is based on the assumption that 
the system adiabatically follows its ground state upon slow switching of the interactions "on" and 
"off' in the distant past and future, respectively. Therefore, the only result of evolving the non- 
interacting ground-state along the entire time axis is acquiring a phase factor e'^. One can then 
compensate for the added evolution segment by dividing this factor out. As the result: {GS|<9|GS) = 
(0|'Z/+oo,/(57Y, _co|0)/e'^ and one faces description of the evolution along the forward time axis without 
the backward segment. However, it comes with the price: one has to take care of the denominator 
(which amounts to subtracting of the so-called disconnected diagrams). 

Such a trick does not work in a non-equilibrium situation with a truly time-dependent Hamil- 
tonian. If the system was driven out of equilibrium, then the final state of its evolution does not 
have to coincide with the initial one. In general, such a final state depends on the peculiarities of the 
switching procedure as well as on the entire history of the system. Thus, one can not get rid of the 
backward portion of the evolution history contained in (|2]i. Schwinger [S) was the first to realize that 
this is not an unsurmountable obstacle. One has to accept that the evolution in the non-equilibrium 
quantum field theory takes place along the closed time contour. Along with the conventional forward 
path, the latter contains the backward path. In this way one avoids the need to know the state of the 
system at f = +oo. 

It is still convenient to extend the evolution in (|2l) to f = -Hoo and back to t. It is important to 
mention that this operation is identical and does not require any additional assumptions. Inserting 
'Ut,+oo^+oa,i = 1 to the left of O in (|2]), one obtains 

0{t)) = — 7TTr{'Z/-co.+co'l/+c«,,mY,,-c.p(-'x))} . (3) 
Tr{p(-oo)) 
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Figure 1 : The closed time contour C. Dots on the forward and the backward branches of the contour 
denote discrete time points. 

Here we also used that, according to the Von Neumann equation, the trace of the density matrix is 
unchanged under the unitary evolution. As a result, we have obtained the evolution along the closed 
time contour C depicted in Figure [1] 

The observable O is inserted at time f, somewhere along the forward branch of the contour. Note 
that inserting the unit operator 'ZY, +co'i/+oo,/ = 1 to the right of (3, we could equally well arrange to 
have an observable on the backward branch of the contour As we show below, the most convenient 
choice is to take a half-sum of these two equivalent representations. The observable may be also 
generated by adding to the Hamiltonian a source term Ho{t) = H(t) + Ori(t)/2, where the plus 
(minus) signs refer to the forward (backward) parts of the contour One needs to calculate then 
the generating functional Z[ri] defined as the trace of the evolution operator along the contour C 
with the Hamiltonian Hoit)- Since the latter is non-symmetric on the two branches, such a closed 
contour evolution operator is not identical to unity. The expectation value of the observable may 
then be generated as the result of functional differentiation {0{t)) - 6Z[Ti]/6ri{t)\,j^Q. We first omit 
the source term and develop a convenient representation for the partition function 

Tr{'l/cp(-oo)) 

Z[0] ^ ' = 1 , (4) 

Tr{p(-oo)} 

where 'Uc - 'U-cx,,+o/U +00-00 - 1 . The source term, breaking the forward-backward symmetry, will 
be discussed at a later stage. Note that since Z[0] - 1, the observable may be equally well written 
in the form, which is more familiar from the equilibrium context: {0{t)) - S In Z[ri]/Sri{t)\ij^o. The 
logarithm is optional in the theory with the closed time contour 

The need to carry the evolution along the two-branch contour complicates the non-equilibrium 
theory in comparison with the equilibrium theory. The difficulties may be substantially reduced by 
a proper choice of variables based on the forward-backward symmetry of the theory. There are also 
good news: there is no denominator e'^, unavoidably present in the single-branch contour theory. 
(One should not worry about Tr{p(-oo)) in (|4|. Indeed, this quantity refers entirely to f = -00, 
before the interactions were adiabatically switched "on". As a result, it is trivially calculated and 
never represents a problem.) The absence of the denominator dramatically simplifies description of 
systems with the quenched disorder It is the denominator, e'^, which is the main obstacle in per- 
forming the disorder averaging of the expectation values of observables. To overcome this obstacle 
the rephca 1251 l26l |271 and the supersymmetry IISTl [32l tricks were invented. In the closed time 
contour theory the denominator is absent and thus there is no need in any of these tricks. 

2 Bosons 

2.1 Partition function 

Let us consider the simplest many-body system: bosonic particles occupying a single quantum state 
with energy wq. Its secondary quantized Hamiltonian has the form 

// = wo P b , (5) 
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where and b are bosonic creation and annihilation operators with the commutation relation [b,b'] - 
1. Let us define the partition function as 



Z = 



Tr{p) 



(6) 



If one assumes that all external fields are exactly the same on the forward and backward branches 
of the contour, then 'Z/ c = 1 and, therefore, Z = 1 . The initial density matrix p - p{H) is some 
operator-valued function of the Hamiltonian. To simplify the derivations one may choose it to be the 
equilibrium density matrix, po = exp{-j6(H - nN)] - exp{-y6(cLio - nWb]. Since arbitrary external 
perturbations may be switched on (and off) at a later time, the choice of the equilibrium initial 
density matrix does not prevent one from treating non-equilibrium dynamics. For the equilibrium 
initial density matrix one finds 



Tr{po) = £e-'^<'^°-''"' = [l-p(c^o)]- 



(7) 



where p((jJo) = e"^*"»"''\ An important point is that, in general, Tr{p) is an interaction- and disorder- 
independent constant. Indeed, both interactions and disorder are supposed to be switched on (and 
off) on the forward (backward) parts of the contour sometime after (before) f = -oo. This constant 
is, therefore, frequently omitted without causing a confusion. 

The next step is to divide the C contour into {2N - 2) time-steps of length 6,, such that - 
= and = t^+i = +oo as shown in Figure[T] One then inserts the resolution of unity in the 
over-complete coherent state basi^ llMl 



1 



d(Re0/)d(Im0;) 



(8) 



at each point j - 1,2,..., 2N along the contour. For example, for N = 3 one obtains the following 
sequence in the expression for TrCZ/ cPo) (read from right to left): 



<06l'Z/-5,l<^5><05l'W-*,l<^4)<04|i|03)<03|1/+<5,l02)<02l'Z/+5,l01><0llpo|06), 



(9) 



where 'ZYj-^^ is the evolution operator during the time interval (5, in the positive (negative) time direc- 
tion. Its matrix elements are given by: 



(0,+, \€(^s,\<f>j) = * |(1 + iH{bKb)6r\ <f>j) 



(10) 



where the approximate equalities are valid up to the linear order in 5,. Obviously this result is not 
restricted to the toy example Q, but holds for any normally ordered Hamiltonian. Note that there 
is no evolution operator inserted between and t^+i- Indeed, these two points are physically indis- 
tinguishable and thus the system does not evolve during this time interval. Employing the following 
properties of coherent states: {4>\<p'} - exp{00') along with (0|e"^^'^°"'^'* - exp {00'p(wo)), and 
collecting all of the matrix elements along the contour, one finds for the partition function Q, 



Z = 



Tr{po 



//n 



d(Re0;)d(Im0;) 



exp 



2^ 



V j,j'=i 



(11) 



'The Bosonic coherent state {(p) (((f>\), parameterized by a complex number 0, is defined as a right (left) eigenstate 
of the annihilation (creation) operator: b\<l>) = ip\ip) ((<l>\b' = {(p\(j>). Matrix elements of a normally ordered opera- 
tor, such as Hamiltonian. take the form {ip\H(b^ ,h)\(l>') = H(<fi,<fi')(<fi\<l)'). The overlap between two coherent states is 
(ip\ip') = exp{<i><p']. Since the coherent state basis is over-complete, the trace of an operator. A, is calculated with the weight: 
Tr{Al = JT-^ ff £/(Re(/.)c/(Im<^) (<l>\A\<p). 
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where the 2N x 2N matrix iG . , stands for 



-1 




1-/1 -1 




\-h -1 




1 


-1 




1+/7 -1 




\+h -1 



(12) 



and h = iojoSt- It is straightforward to evaluate the determinant of such a matrix 



Det[/G-'] =(-l)2'^-p(^^o)(l-n 



1 - p{loq) e 



(wo(5,)-(iV-I) 



1 - p((yo) , 



(13) 



where one used that <5f ^ if ^ oo (indeed, the assumption was 6,f'^ const). Employing 
the fact that the Gaussian integral in (fTTl i is equal to the inverse determinant of ;G"' matrix, see 
Appendix IaI along with one finds 



Det-'[/G-'] ^ 
Trjpo) 



(14) 



as it should be, of course. Note that keeping the upper-right element of the discrete matrix (fT2l i is 
crucial to maintaining this normalization identity. 

One may now take the limit — > oo and formally write the partition function in the continuum 
notations, — > 0(f), as 



Z= Jd[00] exp(/5[0,0]) = jD[#]expJ/£dr[0(f)G-V(O] 
where according to (fTTT )- ( fT2l i the action is given by 



2^ 



with (Jfj = tj 



±6,. Thus, the continuum form of the operator G ' is 



idi 



■ Wo . 



(15) 



(16) 



(17) 



It is important to remember that this continuum notation is only an abbreviation that represents the 
large discrete matrix (fT2] i. In particular, the upper-right element of the matrix (the last term in (fTSIl). 
which contains the information about the distribution function, is seemingly absent in the continuum 
notation ( fTTI ). 

To avoid integration along the closed time contour, it is convenient to split the bosonic field 0(f) 
into the two components 0+(f) and 0-(f) that reside on the forward and the backward parts of the 
time contour, respectively. The continuum action may be then rewritten as 



CO 



df [0+(f)(/5, - wo)0+(f) - 0-(f)(/5, - wo)0-(O] , 



(18) 



where the relative minus sign comes from the reversed direction of the time integration on the back- 
ward part of the contour. Once again, the continuum notations are somewhat misleading. Indeed, 
they create an undue impression that 0+(f) and 0-(f) fields are completely uncorrected. In fact, they 
are correlated owing to the presence of the non-zero off-diagonal blocks in the discrete matrix (fT2l i. 
It is therefore desirable to develop a continuum representation that automatically takes into account 
the proper regularization. We shall achieve it in the following sections. First, the Green's functions 
should be discussed. 
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2.2 Green's functions 



According to the basic properties of the Gaussian integrals, see Appendix lAl the correlator of the 
two bosonic fields is given by 



f 2N 



D[#] (pj(/>j'exp 



iG ! 



(19) 



Note the absence of the factor Z ' in comparison with the analogous definition in the equilibrium 
theory 1341 . Indeed, in the present construction Z = 1. This seemingly minor difference turns out to 
be the major issue in the theory of disordered systems (see further discussion in Section|6l devoted 
to fermions with the quenched disorder). Inverting the discrete matrix in (fT2l ). one finds 





1 


pe'' 


pe^'' 


pe^/i 


pe* 


P 




e-" 


1 


pe'' 


pe'' 


P 


pe-'' 


1 




e-" 


1 


P 


pe-'' 


pe-2" 


1-p 






1 


1 


pe-" 


pe-^" 




e-" 


1 




e" 


1 


pe-" 




_ 1 


e" 






e" 


1 



(20) 



where p s p(a)o), and following the discussion after ( fTSl ), we have put (1 i/i)^ e*-''" and (1 -/i^)-' 
1 . In terms of the fields ^ y± (hereafter j - I, . . .,N and therefore the 2N x 2A^ matrix above is labeled 



as 1,...,A' - l,N,N,N - !,...,!) the corresponding correlators read as 

('Pj+'Pj'-} = iG)j, = riB exp{-(y - j')h} , 
{(l>j-~^y^) = iG)y = («B + 1) exp{-(; - j')h] , 

0{i-i')iG)^,+9(i' -j)iG<^,, 



bj+4>j'+) = iGjj, = - 5jy 



{4>j-4>j'-} ^ iG\, = \ 6;y + 0U ' - j)iG)j, + 9{j - j ')iG<j, . 



(21a) 
(21b) 

(21c) 
(21d) 



where the bosonic occupation number «b stands for riBioJo) = pi {I - p) and symbols T (respectively 
T) denote time-ordering (respectively anti-time-ordering). The step-function 9{j) is defined such 
that 0(0) = 1/2, so e{j) + 0(-j) = 1. 

Obviously not all four Green's functions defined above are independent. Indeed, a direct inspec- 
tion shows that 

+ G^' - G> - G^ = -i6jj> , (22a) 



G^-G^ ^sign(j-j')iG>-G<) 



(22b) 



where sign(7) = 0(j) - 9{-j). We would like to perform a linear transformation of the fields to benefit 
exphcitly from these relations. This is achieved by the Keldysh rotation 



1 , 



(23) 



with the analogous transformation for the conjugated fields. The superscripts "cl" and "q" denote 
the classical and the quantum components of the fields, respectively. The rationale for this notation 
will become clear shortly. First, a simple algebraic manipulation with (l2Th)-(l2Ttl) shows that 



jj' 



G«., 

JJ' 



(24) 
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where hereafter a,l3 - (cl, q). The expUcit form of the {q, q) element of this matrix is a manifestation 
of identity ( I22al i. Superscripts R,A and K denote the retarded, advanced and Keldysh components 
of the Green's function, respectively. These three Green's functions are the fundamental objects of 
the Keldysh technique. They are defined as 

G« , = -/<05' -4,],) = e(j - i ') (g>., - G<.,) = -tea - i ') e-'J-j'^'^ , (25a) 
Gjj, = -i{cf>] -f;,) = e(j ' - i) (g<., - G>.,) = i0{j ' - ;>-(>-^'>" , (25b) 

Gn' = -'<'^5'^y'> = + + G<^. = 5^^, - i{2nB + l)e<'-i''' . (25c) 

Since by definition [G*] ' = -G* [cf. (ISTTil. one may notice that 

G'^ = [G*]"'", G^ = -[G'^]"'" . (26) 

The retarded (advanced) Green's function is lower (upper) triangular matrix in the time domain. 
Since a product of any number of triangular matrices is again a triangular matrix, one obtains the 
simple rule; 

Gf o Gf o . . . o Gf = G* , (27a) 

Gf o G^ o . . . o Gf = G"* , (27b) 

where the circular multiplication sign is understood as a convolution in the time domain (i.e. it 
implies integration over an intermediate time). 

One can now take the continuum limit (A^ — > 00, while N5t const) of the Green's functions. 
To this end, one defines tj - jS, and notices that exp{-(y - j')h] exp{-/(yo(f - ?'))■ A. less trivial 
observation is that the factors 6jj', see (l24l i and ( 1251 ). may be omitted in the continuum limit. The 
reason for this is twofold: (i) all observables are given by the off-diagonal elements of the Green's 
functions, e.g. the mean occupation number at the moment tj is given by {ng{tj)} = 'Gj^.^j = 

(ii) the intermediate expressions contain multiple sums (integrals) of the form Y^jj' ^H'Gj' j — > 
6^N ^ 0. As a result the proper continuum limit of the relations derived above is 

- imt) nt')) = G"\t, r') = ( ;;j ^"^i; ^'^ ) , (28) 

where 

= -i0{t - t') ^ (e - wo + '0)"' , (29a) 

= /6»(f'-f)e"'"°""'''^(e-wo-iO)"', (29b) 

G^ = -/ [2nB(wo) + 1] e-''^'>^'-'"> -2m[2nB{e) + l]6{e - loq) . (29c) 

The Fourier transforms with respect to f - f' are given for each of the three Green's functions. An 
important property of these Green's functions is [cf. ( l25b l 

G*(f, t) + G^{t, f) = . (30) 

It is useful to introduce graphic representations for the three Green's functions. To this end, let us 
denote the classical component of the field by a full line and the quantum component by a dashed 
line. Then the retarded Green's function is represented by a full arrow-dashed line, the advanced by 
a dashed arrow-full line and the Keldysh by full arrow-full line, see Figure|2] Note that the dashed 
arrow-dashed line, which would represent the (0*0'^) Green's function, is absent in the continuum 
limit. The arrow shows the direction from cfi" towards 4>^- 

Notice that the retarded and advanced components only contain information about the spectrum 
and are independent of the occupation number, whereas the Keldysh component does depend on it. 
In thermal equilibrium p - e^^"^, while «b - {e^"^ - 1)"' and therefore 

G^ie) - [g«(6) - G^(e)] coth^^ . (31) 
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<t>"(f) ^'(t) 



<P"(t') 



Figure 2: Graphic representation of G^, G'*, and G^. The full line represents the classical field 
component 0'', while the dashed line the quantum component 0*. 



The last equation constitutes the statement of the fluctuation-dissipation theorem (FDT). The FDT 
is, of course, a general property of thermal equilibrium that is not restricted to the toy example 
considered here. It implies the rigid relation between the response and correlation functions in 
equilibrium. 

In general, it is convenient to parameterize the anti-Hermitian, see ( |26] l. Keldysh Green's func- 
tion by a Hermitian matrix F - F\ as follows 



G'^ ^G'^ o F - F oG^ 



(32) 



where F - F{t, f), and the circular multiplication sign implies convolution. The Wigner transform 
(see below), f (r, e), of the matrix F is referred to as the distribution function. In thermal equilibrium 
f(e) = coth(e/2r) ( see (O)- 



2.3 Keldysh action and causality 

One would like to have a continuum action, written in terms of (p'^', (p'', that properly reproduces the 
correlators ( |28] ) and ( [29] l. To this end, one formally inverts the correlator matrix ( |28] ), and uses it in 
the Gaussian action 



[G;l,r ^ 
[ [G;}.f [G;},f 



where 



^^-ljR(A) ^ ^QRWyl ^ e _ ^„ + /O ^ 6,,,' (id, - OJQ ± iO) , 



(33) 

(34a) 
(34b) 

where we used that the Fourier transform of e is 6,j'id, and parametrization ( |32] | was employed in the 
last line. It is important to mention that the actual discrete matrix action (fTTI) - (fT2l) . being transformed 
to 0'^', (f)'' according to ( |23] ), does not have the structure of ( |33] ). The action ( |33] | should be viewed as 
a formal construction devised to reproduce the continuum limit of the correlators according to the 
rules of the Gaussian integration. It is, however, fully self-consistent in the following sense: (i) it 
does not need to appeal to the discrete representation for a regularization; (ii) its general structure is 
intact in every order of the perturbative renormalization. 

Here we summarize the main features of the action ( l33T l, which, for lack of a better terminology, 
we call the causality structure. 

• The cl — cl component is zero. It reflects the fact that for a pure classical field configuration 
(0? = 0) the action is zero. Indeed, in this case (p+ = (p^ and the action on the forward part of 
the contour is canceled by that on the backward part (safe for the boundary terms, that may be 
omitted in the continuum limit). The very general statement is, therefore, that 



S[0'',O]^O. 



(35) 



Obviously this statement should not be restricted to the Gaussian action of the form given by 
Eq. (O. 
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• The cl-q and q-cl components are mutually Hermitian conjugated upper and lower (advanced 
and retarded) triangular matrices in the time representation. This property is responsible for 
the causality of the response functions as well as for protecting the cl - cl component from a 
perturbative renormaUzation (see below). Relation ( l30b is necessary for the consistency of the 
theory. 

• The q - q component is an anti-Hermitian matrix [cf. (|26bl. In our example [G^] ' - iOF, 
where F is a Hermitian matrix, with a positive-definite spectrum. It is responsible for the 
convergence of the functional integral. It also keeps the information about the distribution 
function. 



2.4 Free bosonic fields 



It is a straightforward matter to generalize the entire construction to bosonic systems with more than 
one degree of freedom. Suppose the states are labeled by an index k, that may be, e.g., a momentum 
vector. Their energies are given by a function a>k, for example - k-/(2m), where m is the mass 
of bosonic atoms. One introduces next a doublet of complex fields (classical and quantum) for every 
state k , (0'^'(k, t), <p^'(k, ?)), and writes down the action in the form of ( |33] ) including a summation 
over the index k. Away from equilibrium, the Keldysh component may be non-diagonal in the index 
k: F - F{k, k'; f, t'). The retarded (advanced) component, on the other hand, has the simple form 

If k is momentum, it is instructive to perform the Fourier transform to the real space and to deal 
with (0'^'(r, t), 0*(r, f))- Introducing a combined time-space index x - (r, t), one may write down for 
the action of the free complex bosonic field (atoms) 



Jj dxdx'(^j',^?) 







'Pi 



(36) 



where in the continuum notations 



[G^f'^^J^^jc, x'-)^6{x-x ') [id, + -^-5; + ^] 

\ 2m / 



(37) 



while in the discrete form it is a lower (upper) triangular matrix in time (not in space). The [G"']'^ 
component for the free field is only the regularization factor, originating from the (time) boundary 
terms. It is, in general, non-local in x and x', however, being a pure boundary term it is frequently 
omitted. It is kept here as a reminder that the inversion, G, of the correlator matrix must posses the 
causality structure ( |28] l. We have introduced the chemical potential yu into ( |37] i. understanding that 
one may want to consider an effective Hamiltonian H - [iSl , where .^V is the total particle number 
operator The new term may be considered as a mean to enforce a ceratin particle number with the 
help of the Lagrange multiplier /i. For discussion of real bosonic fields see AppendixlB] 



3 Collisions and kinetic equation for bosons 
3.1 Interactions 

The short range two-body collisions of bosonic atoms are described by the local four-boson Hamil- 
tonian i/int = /1 2r b'J}\htht, where index r "enumerates" spatial locations. The interaction constant, 
X, is related to a commonly used ^-wave scattering length, a,, as /l = Anaslni (see ||35l ). The 
corresponding term in the continuum Keldysh action takes the form 

5,„t[0+,<^-] fdr r df(#)' = -/i fdr r df [(0+0+)^ - (^_0_)2] . (38) 

J JC J J-oo 
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Figure 3: Graphic representation of the two interaction vertices of the |0|^ theory. There are also two 
complex conjugated vertices with a reversed direction of all arrows. 



It is important to remember that there are no interactions in the distant past, t — -oo (while they are 
present in the future, t - +00). The interactions are supposed to be adiabatically switched on and off 
on the forward and backward branches, respectively. This guarantees that the off-diagonal blocks 
of the matrix ( fT2] i remain intact. Interactions modify only those matrix elements of the evolution 
operator ( fTOl i that are away from f = -co. It is also worth remembering that in the discrete time form 
the (/) fields are taken one time step 5, after the <p fields along the Keldysh contour C. Performing the 
Keldysh rotation ( |23] l. one finds 



dr J dt[4>''0i[if')- + (^''f]+c.c.].. 



(39) 



where c.c. denotes the complex conjugate of the first term. The collision action ( [39] l obviously sat- 
isfies the causality condition ( [35l l. Diagrammatically the action (39[ generates two types of vertices 
depicted in Figures [3](as well as two complex conjugated vertices, obtained by reversing the direc- 
tion of the arrows): one with three classical fields (full lines) and one quantum field (dashed line) 
and the other with one classical field and three quantum fields. 

Let us demonstrate that an addition of the collision term to the action does not violate the funda- 
mental normalization, Z = 1 . To this end, one may expand exp(/5' int) in powers of A and then average 
term by term with the Gaussian action ( [36] l. To show that the normalization, Z - 1, is not altered 
by the collisions, one needs to show that (5 int) = ("^int) - • ■ • - 0. Applying the Wick theorem, one 
finds for the terms that are linear ordering: {^'^ f\f'f + c .c .) ~ {G'^(t,t) + G^(t,t)\G^(t,t) = 0, and 
(^^^'^'(0*)^ + c.c) - 0. The first term vanishes owing to the identity ( [3Q] l, while the second vanishes 
because {(^^'0*) = (even if one appeals to the discrete version ( l24l i. where (0^0*,) = -iSjj' /2 0, 
this term is still identically zero, since it is given by 2j7' ^Jj'^^j'j + ~ ^' 'EOJ)- There are 
two families of terms that are second order in A. The first one is (^^^j'(0j')^<^202 (^2 )^) ~ G^ih, h) 
G\t2,h)[G''{tut2)]\ while the second is {^\^'=l{cl>'=lf (Pl4,'^{^lf) ~ [G\tut2)?G\t2,h)G\t2,h), 
where = df. . Both of these terms are zero, because G*(f2,?i) ~ 0{t2 — h), while G^(f2,fi) ~ 

G^{ti,t2T ~ S(h — h) and thus their product has no support^ It is easy to see that, for exactly the 
same reasons, all higher-order terms vanish and thus the normalization is unmodified (at least in the 
perturbative expansion). 

As another example, consider a real boson field, see AppendixiBl with the cubic non-linearity 



5int= ^ Jdr^dr^^ = ^ Jdr J df [0^ - <^^] = ^ Jdr J df 



(40) 



The causality condition ( [35] l is again satisfied. Diagrammatically the cubic non-linearity generates 
two types of vertices. Figure [U one with two classical fields (full lines) and one quantum field 



^Strictly speaking, G'*(f2, fi) and G^U^Ji) are both simultaneously non-zero at the diagonal: (i = «2. The contribution of 
the diagonal to the integrals, however, is ~ S^N — > 0, when W — > 00. 
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Figure 4: Graphic representation of the two interaction vertices of the 0-' theory. Note the relative 
factor of one-third between them. 



(dashed hne), and the other with three quantum fields. The former vertex carries the factor k, while 
the latter has weight kI3. Note that for the real field the direction of fines is not specified by arrows. 

3.2 Saddle point equations 

Before developing the perturbation theory further, one has to discuss the saddle points of the action. 
According to ( [35] ). there are no terms in the action that have zero power of both (j)'' and (j)''. The same 
is obviously true regarding 6S /6(f>'^' and therefore one of the saddle point equations 

^ = (41) 

may always be solved by 

(D« = , (42) 

irrespectively of what the classical component, cD^', is. By capital letter cD'"^** we denote solutions 
of the saddle point equations. One may check that this is indeed the case for the action given by 
e.g. ( l36b plus ( [39l ). Under condition ( |42] | the second saddle point equation takes the form 

^ . ([G«]-' - A l^^'f) O^' = i^id, + ±.dl+f,-A IcD^f j O^' = . (43) 

This is the non-linear time-dependent Gross-Pitaevskii equation, which determines the classical 
field configuration, provided some initial and boundary conditions are specified. 

The message is that among the possible solutions of the saddle point equations for the Keldysh 
action, there is always one with zero quantum component and with classical component that obeys 
the classical (non-linear) equation of motion. We shall call such a saddle point "classical" . Thanks 
to (l35T l and (l42l i. the action on the classical saddle point field configurations is identically zero. As 
was argued above, the perturbative expansion in small fluctuations around the classical saddle point 
leads to a properly normalized partition function, Z - I. This seemingly excludes the possibility 
of having any other saddle points. However, this conclusion is premature. The system may posses 
"non-classical" saddle points, such that <l>'^ 9^ 0. Such saddle points do not contribute to the partition 
function (and thus do not alter the fundamental normalization, Z - I), however, they may contribute 
to observables and correlation functions. In general, the action on a non-classical saddle point is 
non-zero. Its contribution is thus associated with exponentially small (or oscillatory) terms. Exam- 
ples include tunneling, thermal activation (considered in the next section), oscillatory contributions 
to the level statistics, etc. 

Let us develop now a systematic perturbative expansion in deviations from the classical saddle 
point: 0'' - O'' + S<p'^' and (p"^ - + 6(p''. As was discussed above, it does not bring any new 
information about the partition function. It does, however, provide information about the Green's 
functions (and thus various observables). Most notably, it generates the kinetic equation for the 
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distribution function. To simplify the further consideration, we restrict ourselves to situations where 
no Bose condensate is present, i.e. <!>''' = is the proper solution of the classical saddle point 
equation ( |43] ). In this case = S(f)'' and thus the (5-symbol may be omitted. 



3.3 Dyson equation 

The next goal is to calculate the dressed Green's function, defined as 



^,"P{t, f) = J D[#] r (0 4>''(t') exp (iSo + /5i„t) , (44) 



here a,fS — (cl,q) and the action is given by ( |36] | and ( |39] l. To this end, one may expand the 
exponent in powers of 5 int. The functional integration with the remaining Gaussian action is then 
performed using the Wick theorem, see AppendixlAl This leads to the standard diagrammatic series. 
Combining all one-particle irreducible diagrams into the self-energy matrix £, one obtains 

G = G + Go|;oG + Go|;oGo|;oG + ... = Go(i+|;oG), (45) 

where G is given by ( |28] ) and the circular multiplication sign implies convolution in times and space 
domain as well as a 2 x 2 matrix multiplication. The only difference compared with the textbook 
diagrammatic expansion lfT2l [181 [341 is the presence of the 2x2 Keldysh matrix structure. The fact 
that the series is arranged as a sequence of matrix products is of no surprise. Indeed, the Keldysh 
index, a - (cl, q), is just one more index in addition to time, space, spin, etc. Therefore, as with any 
other index, there is a summation over all of its intermediate values, hence the matrix multiplication. 
The concrete form of the self-energy matrix, £, is specific to the Keldysh technique and is discussed 
below in some details. 

Multiplying both sides of ([45] ) by G ' from the left, one obtains the Dyson equation for the exact 
dressed Green's function, G, in the form 

(G"' - s) o G = i , (46) 

where 1 is a unit matrix. The very non-trivial feature of the Keldysh technique is that the self-energy 
matrix, 2, possesses the same causality structure as G"' (see Eq. ([33])), namely 

£ = ( °« ^^), (47) 

where l/^^^^ are lower (upper) triangular matrices in the time domain, while S*^ is an anti-Hermitian 
matrix. This fact is demonstrated below. Since both G ' and £ have the same structure, one con- 
cludes that the dressed Green's function, G, also possesses the causahty structure, like ([28] l. As a 
result, the Dyson equation acquires the form 



[G-^]"' - Z-^ \ / G^ G^ 



= 1, (48) 



where we have taken into account that [G ^^'^ is a pure regularization (~ iOF) and thus may be 
omitted in the presence of a non-zero 1.^ . Employing the specific form of [G*^'*'] (see ([37b). one 
obtains for the retarded (advanced) components 

[id, + —dl + yK - if^^A o G*''^) = 6{t - t')6(r - r') . (49) 
\ 2m j 

Provided the self-energy component Z**'*' is known (in some approximation), equation ( [49] l consti- 
tutes a closed equation for the retarded (advanced) component of the dressed Green's function. The 
latter carries the information about the spectrum of the interacting system. 
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Figure 5: Self-energy diagrams for the theory. 



To write down equation for the Keldysh component we parameterize it as - o F - F o G'^, 
cf. ( |32] |. where F is a Hermitian matrix in the time domain. The equation for the Keldysh component 
then takes the form ([G*]"' - S^) o (G* o F - F o G"^) = S'^ o G"^. Multiplying it from the right by 
{[G"^] - ^'^) and employing (EH), one finally finds 



(50) 



where [ , ] denotes for the commutator This equation is the quantum kinetic equation for the 
distribution matrix F. Its left-hand side is called the kinetic term, while the right-hand side is 
the collision integral (up to a factor). As is shown below, 2^ has the meaning of an "incoming" 
term, while 2* o F - F o E'^ is an "outgoing" term. In equilibrium these two channels cancel each 
other (the kinetic term vanishes) and the self-energy has the same structure as the Green's function: 
= o F - F o S"^. This is not the case, however, away from the equilibrium. 



3.4 Self-energy 

Let us demonstrate that the self-energy matrix, S, indeed possesses the causality structure ( |47] ). To 
this end, we consider the real boson field with the K(f)^ non-linearity ( |40l ), and perform calculations 
up to the second order in the parameter k. Employing the two vertices of figure |4] one finds the 
following. 

(i) The cl-cl component is given by the single diagram, depicted in FigurelS^. The corresponding 
analytic expression is 2'''"'''(f, f') - AiK^G'Xt, t')G\t, t') - 0. Indeed, the product G*(f, t')G^{t, t') has 
no support (see footnote in section [TTT l. 

(ii) The cl- q (advanced) component is given by the single diagram, FigureH)?. The correspond- 
ing expression is 

I.Ht, f) = AiK^G\t, t')G^{t, t') . (51) 

Since £'*(?, t') ~ G^{t, f) ~ 0{t' - t), it is, indeed, an advanced (upper triangular) matrix. There is a 
combinatoric factor of four, associated with the diagram (four ways of choosing external legs x two 
internal permutations x 1/(2!) for having two identical vertices). 

(iii) The q- cl (retarded) component is given by the diagram of Figure |5}; 

lf{t, t') = WG'^it, t')G'^{t, t') , (52) 

that could be obtained, of course, by the Hermitian conjugation of (ISTT i with the help of (|26] |: = 
[S'*]^ Since S*(f, t') ~ G'^(t, t') ~ 6{t - t'), it is indeed a retarded (lower triangular) matrix. 

(iv) The q - q (Keldysh) component is given by the three diagrams. Figure |5}l-l5f . The corre- 
sponding expression (sum of these diagrams) is 

2^(f,f') = 2iK\G'^{tj')f + 6i(^^K[G^{t,t')f +6i(^^K[G'^(tj')f 

= 2iK^ ([G^(f, t')f + [G^(f, t') - G^{t, t')f) . (53) 
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The combinatoric factors are two for diagram Figure |5}J, and six for Figure and f. In the last 
equality the fact that G^(f, f)G^{t, t') - 0, owing to the absence of support in the time domain, has 
been used again. Employing (|26] |. one finds "L^ = -[S^]^. 

This demonstrates that the self-energy £ possesses the same structure as G One may check 
that the statement holds in higher orders as well. In (|5T1|-(|53]) one has omitted the spatial coordinates, 
which may be restored in an obvious way. 

3.5 Kinetic equation 

To make further progress in the discussion of the kinetic equation it is convenient to perform the 
Wigner transformation (WT)0. The WT of a distribution function matrix, F(r, r'; f, f'). is a func- 
tion f (R, k; T, e), where r and R are central time and coordinate, respectively. According to the 
definition ( [32] i. the f function appears in a product with - G'^. The latter is a sharply peaked 
function at e = Wk for free particles, while for the interacting systems this is still the case as long 
as quasi-particles are well-defined. One therefore frequently writes f(R, k, t), understanding that 
e = Wk- 

To rewrite the kinetic term (the left-hand side of dSOl )) in the Wigner representation, one notices 
that the WT of id, is e, while the WT of is -k^. Then, e.g., [F,^^]^ ^ [k^f]_ + /Vkk^VRf = 
2;kVRf, where the commutator vanishes, since WTs commute. In a similar way: [F, id,]- — > -idrf- 
If there is a scalar potential V{r)bl.br in the Hamiltonian, it translates into the term -y(0'''0^ + tp^cp'^') 
in the action and thus -V{r) is added to [G'^*'*'] This, in turn, brings the term -[F, y]_ to the 
left-hand side of the Dyson equation ( fSOl l or, after the WT, /EVkf, where E = -VrV is the electric 
field. As a result, the WT of the Dyson equation ( l50t takes the form 

{dr - t^k Vr - E Vk) f (R, k, r) = /coll [f ] , (54) 

where Wk = ^/m and /coii[f] is the WT of the right-hand side of ( l50b (times /). This is the kinetic 
equation for the distribution function. 

For real bosons with the dispersion relation e = Wk, see Appendix iBl the kinetic term takes the 
form [e^ - w^,F]_ — > 2;(ec?T - Wk(Vkii>k)VR) f = 2ie{dj - L'kVR)f, where Uk = VkWk is the group 
velocity. As a result, the kinetic equation takes the form: (dj - Wk Vr) f(R, k, r) = /coii[f], where the 
collision integral /coii[f] is the WT of the right-hand side of ( fSOl l. divided by -lie. 

Let us discuss the collision integral now, using the cp^ theory calculations of Section 13.41 as an 
example. To shorten the algebra, let us consider a system that is spatially uniform and isotropic in 
the momentum space. One thus focuses on the energy relaxation only. In this case the distribution 
function is f (R, k, t) = f(T, Wk) = f (^, e), where the dependence on the modulus of the momentum is 
substituted by the Wk = e argument. Employing (BTTi-dSSll. one finds for the WT of the right-hand 
side ofldlOli: 

I* o F - F o Z-* -> -2/ f (r, e) J dco M(t, e, co) [f (r, e - co) + f (t, co)\ , (55a) 
i:^ -2i J dco M{t, e, co) [f(T, e - w)f (r, co) + l], (55b) 

^ The Wigner transform of a matrix A(r, r') is defined as a{R, k) = J dri A + ^ ~ ^) exp{/kri ). One may show 
that the Wigner transform of the matrix C = A o B, which means C(r, r') = / dr"A(r, r")B(r", r'), is equal to 

c(R,k) = jj dndrz jj fl(R+ y,k + ki)fe(R+ ^,k + k2)exp{/(kir2-k2ri)). 

Expanding the functions under the integrals in k, and r,, one finds: c(R, k) = a(R, k) b(R, k)+(2()"' (VRfiVki'-Vii«VRfo)+. . . . 

''Only products of WTs are retained, while all the gradient terms ai'e neglected, in particular G'' — > f (g" - g''). The 
energy-momentum representation is used, instead of the time-space representation as in j5T)-(53), and in the equation for 
E'' o F - F o Z'* one performs a symmetrization between the w and e - w arguments. 
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where the transition rate is given by 




r, e - w; k - q) \{t, uj\ q) . 



(56) 
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Here Ag = /(g^ - g'*)/(27r) and g'*''^'(T, e; k) are the WT of the retarded (advanced) Green functions 
Q«(A) Qjjg jj^g substituted the dressed Green's functions into (l5Tli-(l53]l instead of the bare ones 
to perform a partial resummation of the diagrammatic series. (This trick is sometimes called the 
self-consistent Born approximation. It still neglects the vertex corrections.) Assuming the existence 
of well-defined quasi-particles at all times, one may regard Ag(T, e, k) as a sharply peaked function 
at 6 = Wk- In this case (|56] | simply reflects the fact that an initial particle with e = decays into 
two real (on mass-shell) particles with energies w = and e- oj - Wk-q- As a result, one finally 
obtains for the kinetic equation 



where the time arguments are suppressed for brevity. Due to the identity: coth(fl - b) coth(i') + 1 = 
coth(a)( coth(fl - h) + coth(^)), the collision integral is identically nullified by f(e) = coth(e/2r) 
where T is a temperature. This is the thermal equilibrium distribution function. According to the 
kinetic equation (|57] |. it is stable for any temperature (the latter is determined either by an external 
reservoir, or, for a closed system, from the conservation of total energy). Since the equilibrium 
distribution obviously nullifies the kinetic term, according to ( fSOl l the exact self-energy satisfies 
Y.^ - coth(e/2r)[Z* - S'*]. Since also the bare Green's functions obey the same relation dSTT i. one 
concludes that in thermal equilibrium the exact dressed Green's function satisfies 



This is the statement of Has. fluctuation-dissipation theorem (FDT). Its consequence is that in equilib- 
rium the Keldysh component does not contain any additional information with respect to the retarded 
component. Therefore, the Keldysh technique may be, in principle, substituted by a more compact 
construction: the Matsubara method. The latter does not work, of course, away from equilibrium. 

Returning to the kinetic equation ( fSTT i. one may identify "in" and "out" terms in the collision 
integral. It may be done by writing the collision integral in terms of the occupation numbers n^, 
defined as = 1 + 2 . The expression in the curly brackets on the right-hand side of (|57] i takes the 
form 4 [iif-^n^j - wj^xv^^i^ + n^^ + 1)]. The first term nE_„n„ gives a probability that a particle with 
energy e- o) absorbs a particle with energy <y to populate a state with energy e: this is the "in" term 
of the collision integral. It may be traced back to the Y.^ part of the self-energy. The second term 
-nf(nf_jj + + 1) says that a state with energy e may be depopulated either by stimulated emission 
of particles with energies e- uj and w, or by spontaneous emission. This is the "out" term, that may 
be traced back to the E'^*'*' contributions. 

Finally, let us discuss the approximations involved in the Wigner transformations. Although ( fSOl l 
is formally exact, it is very difficult to extract any useful information from it. Therefore, passing to 
an approximate, but much more tractable, form such as ( |54] i or d57] i is highly desirable. In doing 
it, one has to employ the approximate form of the WT. Indeed, a formally infinite series in VkVR 
operators is truncated, usually by the first non-vanishing term. This is a justified procedure as long as 
(5k » 1, where (5k is a characteristic microscopic scale of the momentum dependence of f, while 
(5R is a characteristic scale of its spatial variations. One may ask if there is a similar requirement 
in the time domain: (5e5T » 1, with de and 5t being the characteristic energy and the time scales 
of f, correspondingly. Such a requirement is very demanding, since typically 6e T and at low 
temperature it would only allow very slow processes to be treated: with 5t » 1 /T. Fortunately, this 
is not the case. Because of the peaked structure of Ag(e, k), the energy argument e is locked to Wk 
and does not have its own dynamics as long as the peak is sharp. The actual criterion is therefore that 




(57) 



= (G* - G-^) coth — . 



(58) 
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6e is much larger than the width of the peak in Ag(e, k). The latter is, by definition, the quasi-particle 
life-time, Tqp , and therefore the condition is Tqi, » l/T. This condition is indeed satisfied by many 
systems where the interactions are not too strong. 



4 Particle in contact with an environment 
4.1 Quantum dissipative action 

Consider a particle with the coordinate 0(f), placed in a potential f/(<l>) and attached to a harmonic 
string if{x, t). The particle may represent a collective degree of freedom, such as the phase of a 
Josephson junction or the charge on a quantum dot. On the other hand, the string serves to model 
a dissipative environment. The advantage of the one-dimensional string is that it is the simplest 
continuum system, having a constant density of states at small energies. Owing to this property it 
mimics, for example, interactions with a Fermi sea. A continuum reservoir with a constant density 
of states at small energies is sometimes called an "Ohmic" environment (or bath). The environment 
is supposed to be in thermal equilibrium. 

The Keldysh action of such a system is given by the three terms S - Ssa + where (see 
Appendix IbI) 



(59a) 



X d O 
df -2 (f — - ^ {^'' + + ^e^'' - ^'"> 
OO [ 

X+oo p 
df J ±x0''{xj)b-'0{x,t), (59b) 

X+oo 
At6\t)&,d,0{x,t)\ . (59c) 

Here we have introduced vectors of classical and quantum components, e.g., = (O'^', O*) and the 
string correlator, D"', that has typical bosonic form ( |36] |. with [D*'^'^']"' = -d^ + v].d\, which follows 
from (142 11 1. The S p represents a particle (see corresponding discussion in (1417b in Appendix|B]i. The 
5str is the action of the string ( I421l i. The interaction term between the particle and the string is taken 
to be the local product of the particle coordinate and the string stress at jc = (so the force acting 
on the particle is proportional to the local stress of the string). In the time domain the interaction is 
instantaneous, <^{t)dxf{x, f)|,v=o — > ^+dxf+ — 0_5.v</'- on the Keldysh contour Transforming to the 
classical-quantum notation leads to 2{<i>'^'dxf^' + '^>''dx<p^'), that satisfies the causality condition i35[ . 
In the matrix notation it takes the form of (|59b). The interaction constant is denoted sjy. 

One may now integrate out the degrees of freedom of the harmonic string to reduce the prob- 
lem to the particle coordinate only. According to the standard rules of Gaussian integration (see 
AppendixlAl). this leads to the so-called dissipative action for the particle 

dtdt'6^(t)f)-\t-t'mt'), (60a) 

(f - f ') = - d, d,,b{x-x'-t-t <T, . (60b) 

The straightforward matrix multiplication shows that the dissipative correlator X) ' possesses the 
standard causality structure. Fourier transforming its retarded (advanced) components, one finds 



k 



where we put y , = 1 for brevity. The e-independent constant (same for R and A components) may 
be absorbed into the redefinition of the harmonic part of the potential {/(O) = const + . . . and. 
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thus, may be omitted. In equilibrium the Keldysh component of the correlator is set by the FDT 

[l^-'f(€) = ([^«]-' - [1^^]-') coth ^ = ie coth ^ . (62) 

It is an anti-Hermitian operator with a positive-definite imaginary part, rendering convergence of 
the functional integral over <i>. 

In the time representation the retarded (advanced) component of the correlator takes a time-local 
form [I)^*'^^] ' = + j6(t - f) d,. On the other hand, at low temperatures the Keldysh component is a 
non-local function, that may be found by the inverse Fourier transform of ( |62] i: 

ttT^ 



{IT + C)5(t - t') ■ 



sinh-[;rr(f-f')] 



(63) 



where the infinite constant C serves to satisfy the condition J dr[D '(f)]*^ - '(f = 0)]*^ = 2iT. 
Finally, for the Keldysh action of the particle connected to a string, one obtains 



X+oo 
( 
DO 



2 df 



T{mt)r + 



2 X. 



df 



, (0«(f) - 0«(f')) 



(64) 



sinh^[;rr(f - f')] 

This action satisfies all the causality criterions listed in Sec. 12.31 Note that, in the present case, 
the Keldysh {q - q) component is not just a regularization factor, but rather a quantum fluctuations 
damping term, originating from the coupling to the string. The other manifestation of the string is 
the presence of the friction term, ~ ydt in the R and the A components. In equilibrium the friction 
coefficient and fluctuations amplitude are connected rigidly by the FDT. The quantum dissipative 
action ( l64l i is a convenient playground to demonstrate various approximations and connections to 
other approaches. 



4.2 Classical limit 

The classical saddle point equation (the one that takes <f (f) - 0) has the form; 
1 6S [O] 



2 6^1 



= — ~ 1 1 ; — = . (65) 

<i,.=o df2 2 df 5<D^' ^ 



This is the deterministic classical equation of motion. In the present case it happens to be Newton 
equation with the viscous force -(y/2)<i)'''. This approximation neglects both quantum and thermal 
fluctuations. 

One may keep the thermal fluctuations, while completely neglecting the quantum ones. To this 
end, it is convenient to restore the Planck constant in the action ( |64| ) and then take the limit Ti ^ Q. 
For dimensional reasons, the factor ft ' should stand in front of the action. To keep the part of 
the action responsible for the classical equation of motion (|65l l free from the Planck constant it 
is convenient to rescale the variables as — > M)"^'. Finally, to keep proper units, one needs to 
substitute T T jti in the last term of ( l64b . The limit ft ^ is now straightforward: (i) one has to 
expand [/(<!>''' + ftO*) to the first order in ftO* and neglect all higher order terms; (ii) in the ft — » 
limit, which is equivalent to the T — > oo, the non-local part of the action ( |64l i drops out, while the 
local term oc ((l)'?(f))^ survives. As a result, the classical limit of the dissipative action is 



df2 ■ 2 df ■ d^^' I ^ V- / • (66) 



Physically the limit ft — > means that ftQ «: T, where Q is a characteristic classical frequency of the 
particle. This condition is necessary for the last term of ( l64l i to take the time-local form. The con- 
dition for neglecting the higher-order derivatives of ?7 is ft «: y(<i)''')"', where is a characteristic 
classical amplitude of the particle motion. 
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4.3 Langevin equation 



One way to proceed with the classical action ( |66] | is to note that the exponent of its last term (times 
/) may be identically rewritten in the following way 



exp j-2rr j dt = j m] exp J- j df 



fit) 



2yT 



(67) 



This identity is called the Hubbard-Stratonovich transformation, while ^(/) is an auxiliary Hubbard- 
Stratonovich field. The identity is proved by completing the square in the exponent on the right- 
hand side and performing the Gaussian integration at every instance of time. There is a constant 
multiplicative factor hidden in the integration measure, D[^]. 

Exchanging the order of the functional integration over ^ and <1), one finds for the partition 
function: 



X ^ \y[<^''] J D[(^''] exp^-2i J df<l)«(f) 



d^O'' 



df2 



^ y dO' 



2 df 



(68) 



Since the exponent depends linearly on <l>'?(f), the integration over D[<[)^'] results in the i5-function 
of the expression in the square brackets. This functional <5-function enforces its argument to be zero 
at every instant of time. Therefore, among all possible trajectories <l)^'(f), only those that satisfy the 
following equation contribute to the partition function: 



2 df 



(69) 



This is Newton equation with a time-dependent external force ^(f). Since, the same arguments are 
applicable to any correlation function of the classical fields, e.g. {<l)'^'(f)0^'(f')), a solution strategy is 
as follows: (i) choose some realization of ^(f); (ii) solve ( |69] ) (e.g. numerically); (iii) having its solu- 
tion, 0'''(f), calculate the correlation function; (iv) average the result over an ensemble of realizations 
of the force ^(f). The statistics of the latter is dictated by the weight factor in the D[^] functional 
integral in (|68] |. It states that ^(f) is a Gaussian short-range (white) noise with the correlators 



<f(0) = 0. 



mW)) = yT6{t - t') . 



(70) 



Equation ( 1691 ) with the white noise on the right-hand side is called the Langevin equation. It de- 
scribes classical Newtonian dynamics in presence of stochastic thermal fluctuations. The fact that 
the noise amplitude is proportional to the friction coefficient, j, and the temperature is a manifesta- 
tion of the FDT The latter holds as long as the environment (string) is at thermal equilibrium. 



4.4 Martin-Siggia-Rose method 

In the previous section we derived the Langevin equation for a classical coordinate, (J'', from the 
action written in terms of O'^' and another field, O*. An inverse procedure of deriving the effective 
action from the Langevin equation is known as Martin-Siggia-Rose (MSR) technique [7). It is 
sketched here in the form suggested by DeDominicis [SJ. 
Consider a Langevin equation 

dW]^m^ (71) 

where (9[<[>'^'] is a non-linear differential operator acting on the coordinate (fV), and ^(f) is a white 
noise force, specified by (iTOl i. Define the "partition function" as 

Z[^] = r D[0^'] J[d] 6{dW] - m) = 1 ■ (72) 
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It is identically equal to unity by virtue of the integration of the (5-function, provided J'[d] is the 
Jacobian of the operator <9[<1>' ']. The way to interpret ( l72b is to discretize the time axis, introducing 
A^-dimensional vectors <t>^.' = 0''(fj) and = ^(tj). The operator takes the form: (9, - OijO'j' + 
^rijk<i>'j''S>l' + . . ., where summations are taken over repeated indexes. The Jacobian, J', is given 
by the absolute value of the determinant of the following N x N matrix: J/j = dOijd'^'^! - O/j + 

^ijk'^k ^^^^ possible to choose a proper (retarded) regularization where the Jij matrix is a lower 

triangular matrix with a unit main diagonal (coming entirely from the On = 1 term). One then finds 
that, in this case, J' - I. Indeed, consider, for example, (9[<I)'^'] = 5,0'' - U{<S>'''). The retarded 
regularized version of the Langevin equation is Oj' = O^^j + 6,{U(<S>'j'_^) + ^,-i). Clearly in this case 
Jjj = 1 and Jji-i = -1 - 6,U'{<S>'j'_^), while all other components are zero; as a result^ = 1. 

Although the partition function (fTZt is trivial, it is clear that all meaningful observables and the 
correlation functions may be obtained by inserting a set of factors: (i>'^'{t)<^'^'{t') . . . into the functional 
integral (|72T i. Having this in mind, let us proceed with the partition function. Employing the integral 
representation of the 5-function with the help of an auxiliary field ^^'(t), one obtains 

Z[^] = j D[0^',(D*] expJ-2/ J df (D«(0[O''[cD^'(f)] - ^(f)]j , (73) 

where ff^ stands for the retarded regularization of the O operator and thus one takes ^ = 1 . One 
may average now over the white noise dTOt . by performing the Gaussian integration over ^ 

Z = J D[^]expj-^ J dtfit)^Z[^] 

= J D[(I)"', <!)«] exp^- J" df [2/0'?(r)(9'*[c[)"'(r)] +2yr[0'?(f)]^]j . (74) 

The exponent in ( |74l i is exactly the classical limit of the Keldysh action, cf. ( |66l ), including the 
retarded regularization of the differential operator The message is that MSR action is nothing else 
but the classical (high-temperature) limit of the Keldysh action. The MSR technique provides a 
simple way to transform from a classical stochastic problem to its proper functional representation. 
The latter is useful for analytical calculations. One example is given below. 

4.5 Thermal activation 

Consider a particle in a meta-stable potential well, plotted in Figure |6^. The potential has a meta- 
stable minimum at <1> = 0, and a maximum at O = 1 with the relative hight Uo- Let us also assume 
that the particle's motion is over-damped, i.e. y » VC/". In this case one may disregard the inertia 
term, leaving only viscous relaxation dynamics. The classical dissipative action ( l66b takes the form 



(75) 



The corresponding saddle point equations are 



2 8^'=' 



k^'-^'-^ji:^- (76b) 



These equations possess the classical solution: <^'^{t) = whereas <b'^'{t) satisfies the classical equa- 
tion of motion: \ 6''' - -dU{<i>''')/d<i>''' . For the initial condition O'''(0) < 1 the latter equation 
predicts the viscous relaxation towards the minimum at O'^' = 0. According to this equation, there 
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a) 



b) 



Figure 6: a) A potential with a meta-stable minimum, b) The phase portrait of the Hamiltonian 
system (|77] |. Thick Unes correspond to zero energy, arrows indicate evolution direction. 

is no possibility to escape from this minimum. Therefore, the classical solution of ( |76] | does not 
describe thermal activation. Thus, one has to look for another possible solution of ( |76] l, the one with 
0. 

To this end, let us perform a linear change of variables: <!)'''(?) = q{t) and <[)^'(f) = p(t)/{iy). Then 
the dissipative action (iTSl l acquires the form of a Hamiltonian action iS - - J dt (p q - H{p,q)) 
where the effective Hamiltonian 



is introduced. It is straightforward to see that in terms of the new variables the equations of motion 
( |76] | take the form of the Hamilton equations: q — dH/dp and p - -dH/dq. One needs, thus, to 
investigate the Hamiltonian system with the Hamiltonian ( ITTI i. To visualize it, one may plot its phase 
portrait, consisting of lines of constant energy E = H{p(f),q{t)) on the {p,q) plane, see Figure|6j3. 
The topology is determined by the two lines of zero energy, p - Q and Tp - dU(q)/dq, that intersect 
at the two stationary points of the potential, q = and q = I. The p - Q line corresponds to the 
classical (without Langevin noise) dynamics (note that the action is identically zero for motion along 
this line) and thus q - Q is the stable point, while q — \ is the unstable one. Owing to Liouville 
theorem, every fixed point must have one stable and one unstable direction. Therefore, along the 
"non-classical" line p - T^^dU{q)/dq, the situation is reversed: q = is unstable, while ^ = 1 is 
stable. It is clear now that to escape from the bottom of the potential well, ^ = 0, the system must 
evolve along the non-classical Une of zero energy until it reaches the top of the barrier, q - I, and 
then continue to move according to the classical equation of motion (i.e. moving along the classical 
fine p - 0). There is a non-zero action associated with the motion along the non-classical line 



where one has used that H - along the trajectory. As a result, the thermal escape probability is 
proportional to e'^ = e"^°^^, which is nothing but the thermal activation exponent. 

Amazingly, this trick of rewriting viscous (or diffusive) dynamics in a Hamiltonian form works 
in a wide class of problems, see e.g. |[37l . The price one has to pay is the doubling of the number of 
degrees of freedom: q and p in the Hamiltonian language, or "classical" and "quantum" components 
in the Keldysh language. 

4.6 Fokker-Planck equation 

Another way to approach the action dTSl l is to notice that it is quadratic in O'' and therefore the 
D[0*] integration may be performed explicitly. To shorten notation and emphasize the relation to the 



H{p,q) = p 

7 [ 



dU(q) 
dq 



(77) 




23 



classical coordinate, we follow the previous section and use the notation O^'(f) = q{t). Performing 
the Gaussian integration over of exp {iS [<£]), with S [<!>'', Q>''] given by ( fTSl ), one finds the action, 
depending on <J>'^' = q only 



(78) 



One may now employ the same trick, which allows to pass from the Feynman path integral to the 
Schrodinger equation If38l . Namely, let us introduce the "wave function", Viq, t), that is a result of 
the functional integration of exp(/5 {q\) over all trajectories that at time t + 5t pass through the point 
qN = q- Considering explicitly the last time-step, 6,, integration, one may write TiqN, t + S,) as an 
integral of P{qN-\, t) = P{q - 6q, t) over 5q = q - qN-\'- 



d[6q\ exp 

CO 

A[6q\ exp 

oo 



S, 
'ijT 



2\ 



"sr 6, 



exp 



[~U'q{q-6,)-^{u'qfyiq-S„t) 



(79) 



where the integration measure C is determined by the condition: C J d[6q] exp( - ySq/iSTd,)) - 1. 
Expanding the expression in the square brackets on the right-hand side of the last equation to the 
second order in 6q and the first order in Si, one finds 



\ + —^U + —{u \ — —iu] 



1 i^l) 



6, 



2T 
16, 



2 4T^ 



IjT 



IT 2 



= r{t) + ^{u';qr + u'qVq + T'P';^), 



(80) 



'ST6, 



4Td,/y. 



where 

d[6q]6loxp 

CO 

Finally, rewriting the last expression in the dififerential form, one obtains 

2 d 



'dt 



d dU d'- 

+ T 

dq dq dq~ 



7 dq 



dq dq 



(81) 



This is the Fokker-Planck (FP) equation for the evolution of the probability distribution function, 
"Piq, t). The latter describes the probability to find the particle at a point q - <J>'^' at time t. If 
one starts from an initially sharp (deterministic) distribution, Viq, 0) - 6{q - q{Q)), then the first 
term on the right-hand side of the FP equation describes the viscous drift of the particle in the 
potential U(q). Indeed, in the absence of the second term {T - 0), the equation is solved by Viq, t) - 
6(q - q(t)), where q{t) satisfies the deterministic equation of motiorH (y/2)q(t) — -dU (q(t))/dq . The 
second term on the right-hand side describes the diffusion spreading of the probability distribution 
owing to the thermal stochastic noise ^(f)- For a confining potential U{q) (such that C/(+oo) — > 
oo) the stationary solution of the FP equation is the equilibrium Boltzmann distribution: P{q) ~ 
exp{-U(q)m. 

The FP equation may be considered as the (imaginary time) Schrodinger equation: P - HV, 
where the Hamiltonian, H, is nothing but the "quantized" version of the classical Hamiltonian ( T/TT i. 
introduced in the previous section. The "quantization" rule is /9 — » p 5 -d/dq, so the canonical 



^To check this statement one may substitute P(q,t) = S(q — q(t)) into the T 



FP equation: d'^{_q - qm-qd)) 



• q(t)) + U'^S'^iq - Then multiplying both parts of this equation by q and integrating over dq (by per- 

-(2/y)U'^{q(t)). 



(2/r)[[/,y(g 

forming integration by parts), one finds cj(t) ■ 
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commutation relation [q,p] - 1, holds. Notice that before applying this quantization rule, the 
corresponding classical Hamiltonian must be normally ordered. Namely, the momentum p should 
be to the left of the coordinate q, cf. ( ITTI i. Using the commutation relation, one may rewrite the 
quantized Hamiltonian as 

H^Tp^- pU:, ^t[p- U:J{2T)) [p - U'J(2T)) - (U'f/(4T) + U'^Jl 

(we put y/2 - 1) and perform the canonical transformation Q - q and P - p - U[J{2T). In 
terms of these new variables the Hamiltonian takes the familiar form H - TP^ + V(Q), where 
V{Q) = -(U'q)^/{4T) + C/gg/2, while the "wave function" transforms as P(Q, t) = ef^(e)/(2r)^_ 

4.7 From Matsubara to Keldysh 

In some applications it may be convenient to derive an action in the equiUbrium Matsubara tech- 
nique ifTTlfTSl and change to the Keldysh representation at a later stage to tackle out-of-equilibrium 
problems. This section intends to illustrate how such a transformation may be carried out. To this 
end, consider the following bosonic Matsubara action: 

°° 1 

5[<D„,] =rr^ -|e,„||<l)„,|\ (82) 

ni=-oo 

where e„, = InTm and (S>,„ - — ^^TO(T)e"^"''" are the Matsubara components of a real bosonic 
field, <1)(t), with the periodic boundary conditions 0(0) = <5(y6). Note that, owing to the absolute 
value sign, |e„,| ^ idj- In fact, in the imaginary time representation the kernel K,„ = |e„,| of the 

action (l82T i acquires the form K(t) - Yjm kmk"^'"'^ = C5{t) - nT- 1 sirrinTT), where the infinite 

rP 

constant C is chosen to satisfy dTK(T) = Ko = 0. As a result, in the imaginary time representation 
the action (182b has the following non-local form 

= 2 jj drdr' <D(t)K(t - t')<1>(t') 

= ^ AtAt' (<1)(t) - <5(T')f . (83) 

4;rJJo sin'[7rr(T-T')] ^ ' 

This action is frequently named after Caldeira and Leggett 1361 . who used it to investigate the influ- 
ence of dissipation on quantum tunneling. 

To transform to the Keldysh representation one needs to double the number of degrees of free- 
dom: <[) — > <5 = (O'^', <^''Y . Then according to the causality structure. Section |24l the general form 
of the time translationally invariant Keldysh action is 

where [T)''^'^'(e)]"' is the analytic continuation of the Matsubara correlator |e„,|/2 from the up- 
per (lower) half-plane of the complex variable e,„ to the real axis: -;e„, — > e, see ifTSl . As 
a result, [D*^'*'(e)] ' = +ie/2. In equilibrium the Keldysh component follows from the FDT: 
[T)-\e)f = ([2D*]-i - [D'^]-i)coth(e/2r) = /e coth (6/2 T), cf. (gB and Therefore, the 

Keldysh counterpart of the Matsubara action ( |82] | or ( |83] |. is the already familiar dissipative ac- 
tion ( l64l i. (without the potential and inertial terms, of course). One may now include external fields 
and allow the system to deviate from the equilibrium. 
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4.8 Dissipative chains and membranes 

Instead of dealing with a single particle connected to a bath, let us now consider a chain or a lattice 
of coupled particles, with each one connected to a bath. To this end, we (i) supply a spatial index, 
r, to the field 0(f) 0(r, f), and (ii) adds the harmonic interaction potential between neighboring 
particles: ~ D{^(r,t) - <l)(r + 1,?))^ ^ D(dr^)^ in the continuum limit, where D is the rigidity 
of the chain or membrane. By changing to the classical-quantum components and performing the 
spatial integration by parts [cf. ( 142 lb 1. the gradient term translates to: D (O'^cJjO'^' + O^'cJj'I)'^). Thus, 
it modifies the retarded and advanced components of the correlator, but it does not aff'ect the (q - q) 
Keldysh component: 

[^RWyi = i 5(r - t') Sir - r')( + d, + Ddl). (85) 

In the Fourier representation [X)**'**(k, e)]"' - i( + /e-Dk^). In equilibrium the Keldysh component 
is not aff'ected by the gradient terms, and is given by ( |62] i (in the real space representation it acquires 
the factor 6{r - r')). In particular, its classical Umit is [D"']'^ - i2TS{t - t')6{Y - r'), cf. (|63]l. As a 
result, the action of a classical elastic membrane in contact with a bath is 



5[<D"',ct)«] =2 ^ drdf -O" {d,<^'' - Ddl^'' + ^ g^^f^ j + ilTi^^f 



(86) 



where the inertia terms have been neglected and we put y/2 = 1 for brevity. One may introduce 
now an auxiliary Hubbard-Stratonovich field ^(r, t) and write the Langevin equation according to 
Sectionl4!4l 



5,cD^' - Ddl^^' + = ^(r, t) , (87) 

where ^ is a Gaussian noise with short-range correlations (^(r, t)^{v' , t')) - 2TS(t — t')6(r - r'). 

Let us consider an elastic chain placed in the bottom of the (r-independent) meta-stable potential 
well, depicted in Figure [6^. If a sufficiently large piece of the chain thermally escapes from the 
well, it may find it favorable to slide down the potential, pulling the entire chain out of the well. 
To find the shape of such an optimally large critical domain and its action, let us change to the 
Hamiltonian variables of Section |4Jl qir,t) s (^'(r, f) and p(r,t) = 2/<l)'?(r, f)- The action ( [86] l 
takes the Hamiltonian form iS - - ff drdr (pq - H(p, q)) with 

H^pDdlq-p^^ + Tp\ (88) 
oq 

and the corresponding equations of motion are 

? = ^ = Dd^q - U'^(q) + 2Tp , (89a) 
SH 

P^-—^-Ddlp + pu:;g(q). (89b) 

These are complicated partial differential equations, that cannot be solved in general. Fortunately, 
the shape of the optimal critical domain can be found. As was discussed in Section |4!7l the minimal 
action trajectory corresponds to a motion with zero energy, H -Q. According to Eq. dSSl l, this is the 
case if either p - (classical zero-action trajectory), or Tp - U'^{q) - Dd^q (finite-action escape 
trajectory). In the latter case the equation of motion for ^(r, f) takes the form of the classical equation 
in the reversed time: q - -Dd^q + U'^jiq) = Tp . Thanks to the last equality the equation of motion 
for ^(r, f) is automatically satisfied 0. In the reversed time dynamics the q{r, t) - configuration is 
unstable and therefore the chain develops a "tongue" that grows until it reaches the stationary shape: 

-Ddlq + U:,(q)^0. (90) 



''Indeed, Tp = d,q = -Ddjq + qU^'^ = Ti-Ddjp + pU^'^). This non-trivial fact reflects the existence of an accidental 
conservation law: H{p(r, t), q{r, f)) = - locally\ While from the general principles only the total global energy has to be 
conserved. 
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The solution of this equation with the boundary conditions q(±oo) - gives the shape of the critical 
domain. Once it is formed, it may grow further according to the classical equation q = Dd^q - U'^iq) 
and p - with zero action. The action along the non-classical escape trajectory, paid to form the 
"tongue" is (H(p, q) = 0): 

iS ^- JJ drdtpq^ -1 JJ drdr {-Dd^^q + U'g(q))q = J dr{^{drqf + U{q)), (91) 

where in the last equality an explicit integration over time was performed. The escape action is 
given therefore by the static activation expression that includes both the elastic and the potential 
energies. The optimal domain ( |90l ). is found by the minimization of this static action ( |9T] ). One 
arrives, thus, at a thermodynamic Landau-type description of the first-order phase transitions. Note 
that the effective thermodynamic description appears owing to the assumption that H{p, q) - 0, 
when all the processes take infinitely long time. 



5 Fermions 

5.1 Partition function 

Consider a single quantum state with energy eo- This state is populated by spin-less fermions (parti- 
cles obeying the Pauli exclusion principle). In fact, one may have either zero, or one particle in this 
state. The secondary quantized Hamiltonian of such a system has the form 

// = eoc'e, (92) 

where and c are creation and annihilation operators of fermions on the state eo- They obey standard 
anti-commutation relations: {c,c') - 1 and {c ,c] — {c'^,c') - 0, where { , ) stands for the anti- 
commutator 

One can now consider the evolution operator along the Keldysh contour, C and the corresponding 
partition function, Z = 1, defined in exactly the same way as for bosonic systems (|6|l. The trace of 
the equilibrium density matrix is Trjpo) = 1 + p(eo)> where the two terms stand for the empty and 
the singly occupied states. One divides the Keldysh contour onto {2N — 2) time intervals of length 
5, ~ 1 /A^ — » and introduces resolutions of unity in 2N points along the Keldysh contour, C; see 
Figure[T] The only difference from the bosonic case of Section lZTl is that now one uses the resolution 
of unity in the fermionic coherent state basi^ 



i = J J A^jAilfje-'l'''^'\ilfj){ilfj\, (93) 

where ipj and are mutually independent Grassmann variables. The rest of the algebra goes through 
exactly as in the bosonic case, see Section lZTI As a result, one arrives at 



2N 



Tr{po)_ 



2N 



V JJ'=l 



(94) 



'The fermionic coherent state = (1 - i^c' )|0>, parameterized by a Grassmann number i// (such that {t//, i//'] = [if/, c] = 
0), is an eigenstate of the annihilation operator: cli/r) = Similarly, {il/\c^ = where i// is another Grassmann 

number, unrelated to if/. The matrix elements of a normally ordered operator, such as e.g. the Hamiltonian, take the form 
(iIj\H(c^ ,c)\iIj') = H(ijj, The overlap between any two coherent states is = 1 + ij/i//' = e\p{i//i//']. The trace 

of an operator, 6, is calculated as Tr|(5| = <0|(?|0> + <1|(9|1> = <0|(9|0> + <0|c(9ct|0> = JJ dii/dtl/e-'t"^{-t/i\d\i//), where the 
Grassmann integrals are defined as Jdi/< 1=0 and jd>fi ifi = 1. 
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where the 2N x 2N matrix G . ., is 



-1 


-p(e()) 


\-h -1 




1-/2 -1 




1 


-1 




l+/j -1 




1+/7 -1 



(95) 



and h = ieo6,. The only difference from the bosonic case is the negative sign in front of p(eo) 
matrix element, originating from the minus sign in the {-4'2n\ coherent state in the expression for 
the fermionic trace. To check the normalization, let us evaluate the determinant of such a matrix 



Det[iG"'] = 1 +p(eo)(l -/! ) 



2\N-\ 



1 +p(eo)e' 



1 + p(eo) . 



(96) 



Employing the fact that the fermionic Gaussian integral is given by the determinant (unlike the 
inverse determinant for bosons) of the correlation matrix, (see AppendixlAlfor details), one finds 



Det[;G- 
Tr{po) 



= 1 , 



(97) 



as it should be. Once again, the upper-right element of the discrete matrix ( |95i l is crucial to maintain 
the correct normalization. Taking the limit N oo and introducing the continuum notation, ijjj — > 
if/it), one obtains 



Z = ^D[\pil/] exp(i5[iA,(A]) = ^ D[i/^(/r] exp^/J" df [iA(f) G" V(0] 
where according to ( |94l i and ( |95] ) the action is given by 



2N r 



j=2 



lilfj — eo^Hl^i-\ 



5tj + /iAi[(Ai +p(eo)iA2A'] , 



(98) 



(99) 



withftj = tj-tj-i - +6,. Thus the continuum form of the operator G' is the same as for bosons (fTTji: 
G~' = id, - £{). Again the upper-right element of the discrete matrix (the last term in (|99]l), which 
contains information about the distribution function, is seemingly absent in the continuum notation. 

Splitting the Grassmann field i/r(0 into the two components ip+it) and ij/-{t) that reside on the 
forward and the backward parts of the time contour, respectively, one may rewrite the action as 



f 



(100) 



where the dynamics of and are actually not independent from each other, owing to the presence 
of non-zero off-diagonal blocks in the discrete matrix ( |95] ). 



5.2 Green's functions and Keldysh rotation 

The four fermionic Green's functions: G^^^' and G^*^' are defined in the same way as their bosonic 
counterparts, see ( |2TI ). 

<tA+(OiA-(f')) = iG^(t,t') = -«^exp{-/eo(f-f')), (101a) 



<(A-(f)<A+(f')>= iG^itj') = (l-nF)exp{-/eo(r-f')), 
{tfr+(t)tfr+(t')} = iO'^it, t') = e{t - t')iG^{t, t') + 6(t' - t)iG^{t, t') ., 



(101b) 
(101c) 
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{ifr-(t)if^-(t ')) = iG^(t, t') = 0(t' - tyCit, t') + 0(t - t')iG^{t, f) . (lOld) 

The difference, however, is in the minus sign in the expression for G^, due to the anti-commutation 
relations, and Bose occupation number is exchanged for the Fermi occupation number: iig — » nf = 
p(eo)/(l + p(fo))- Equations ( I22al ) and ( I22bl ) hold for the fermionic Green's functions as well. 

It is customary to perform the Keldysh rotation in the fermionic case in a different manner from 
the bosonic one. Define the new fields as 

<Ai(0 = + -A-CO) , «A2(0 = -^(M) - ^-(t)) ■ (102) 

V2 V2 



Following Larkin and Ovchinnikov 11391 . it is agreed that the bar-fislds transform in a different way: 

hit) = -^{h(t) - hit)) , hit) = -^ihit) + hit)) . (103) 
V2 V2 

The point is that the Grassmann fields iff are not conjugated to tf), but rather are completely indepen- 
dent fields, that may be transformed in an arbitrary manner (as long as the transformation matrix has 
a non-zero determinant). Note that there is no issue regarding the convergence of the integrals, since 
the Grassmann integrals are always convergent. We also avoid the subscripts cl and q, because the 
Grassmann variables never have a classical meaning. Indeed, one can never write a saddle point or 
any other equation in terms of if/, if/, rather they must always be integrated out in some stage of the 
calculations. 

Employing (1102b . ( I103l l along with Eq. ( llOll l. one finds 

- i{^,mbit')) = GAt, t') = iy ^""^q' ^'^ ) > (104) 

where hereafter a,b - (1,2). The fact that the (2, 1) element of this matrix is zero is a manifestation 
of identity ( I22ab . The retarded, advanced and Keldysh components of the Green's function ( 11041) 
are expressed in terms of G^*'^' and G^'^' in exactly the same way as their bosonic analogs ( [25T l, and 
therefore posses the same symmetry properties: (|26]l-(l30ll. An important consequence of ( [27T i and 
([33 is 

Tr {g12 o Gf; o . . . o Gg} {t, = 0, (105) 

where the circular multiplication sign involves convolution in the time domain along with the 2x2 
matrix multiplication. The argument (f, f) states that the first time argument of G*'* and the last 
argument of G*" are the same. 

Note that the fermionic Green's function has a different structure compared to its bosonic coun- 
terpart ( |28] i: the positions of the R,A and K components in the matrix are exchanged. The reason, 
of course, is the different convention for transformation of the bar fields. One could choose the 
fermionic convention to be the same as the bosonic (but not the other way around), thus having the 
same structure (l28T l for the fermions as for the bosons. The rationale for the Larkin-Ovchinnikov 
choice ( |104b is that the inverse Green's function, G ' and fermionic self-energy £f have the same 
appearance as G, namely 

^ [gV )' "^"^l s-J' (106^ 

whereas in the case of bosons G"' (see (|33T l) and £ (see ( |47] i) look differently from G (see (|28]l). 
This fact gives the form (1104b and (1106b a certain technical advantage. 

For the single fermionic state, after the Keldysh rotation, the correlation functions (llOlb allow 
us to find components of the matrix (1104b 



G*(f, t ') = -ie{t - t ')e-'*('-''> ^ (e - eo H- , (107a) 
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G^(t, t ') = i6(t' - t)e-''"'^'-'"> -> (e - eo - '0)"' , (107b) 
G^{t,t') = -/(I - 2«f)e"'*^'"''* ^ -27r/(l - 2«f )5(e - eo) , (107c) 
where the right-hand side provides also the Fourier transforms. In thermal equilibrium, one obtains 

G'^ie) = [G\e) - G\e)\ tanh ^ . (108) 

This is the FDT for fermions. As in the case of bosons, FDT is a generic feature of an equilibrium 
system, not restricted to the toy model. In general, it is convenient to parameterize the anti-Hermitian 
Keldysh Green's function by a Hermitian matrix F = as 

G'^ ^G'^ oF -F oG^ . (109) 

The Wigner transform of F{t, t') plays the role of the fermionic distribution function. 



5.3 Free fermionic fields and their action 

One may proceed now to a system with many degrees of freedom, labeled by an index k. To this 
end, one changes eo — > fk and performs summations over k. If k is a momentum and = k^/(2m), 
it is instructive to transform to the coordinate space representation i/'(k, t) — > i^(r, f), while ek = 
k^/(2m) -dl/{2m). Finally, the Keldysh action for a non-interacting gas of fermions takes the 
form 

5o[<A,<A] = rr dxdx' ^ 4r„{x)[G-\x,x')l,ifri,(x'), (110) 

a.h=l 

where x - (r, f) and the matrix correlator [G"']ni has the structure of ( 1106b with 

[G'^^'^\x, x')r^ = 6{x - x') (id, + —dz +iA. (Ill) 
\ 2m I 

Although in continuum notation the R and the A components look seemingly the same, one has to 
remember that in the discrete time representation, they are matrices with the structure below and 
above the main diagonal, respectively. The Keldysh component is a pure regularization, in the sense 
that it does not have a continuum limit (the self-energy Keldysh component does have a non-zero 
continuum representation). All of this information is already properly taken into account, however, 
in the structure of the Green's function (11041 ). 



5.4 External fields and sources 

According to the basic idea of the Keldysh technique, the partition function Z = 1 is normalized by 
construction, see ( |97] |. To make the entire theory meaningful one should introduce auxiliary source 
fields, which enable one to compute various observable quantities: density of particles, currents, etc. 
For example, one may introduce an external time-dependent scalar potential V(r, f) defined along 
the contour C. It interacts with the fermions &5 Sv - / dr df y(r, t)ij/{r, t)>f/{r, t). Expressing it via 
the field components residing on the forward and backward contour branches, one finds 

//-•+O0 
dr I df - y_iA-iA-] 

%J —OQ 

dr df[r'(0-+iA+-'A-'A-) + V9((A+<A+ + 0--<A-)] 

%J—oo 

dr J dt [r''(iAi<Ai + h^2) + v"(iAi<A2 + h^i)] , (112) 
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where the y"?* components are defined in the standard way for real boson fields, y '*'') = (V+ + 
V-)/2, way. We also performed rotation from to i//i(2) according to (I102l i and ( I103l l. Note that 
the physical fermionic density (symmetrized over the two branches of the Keldysh contour) g = 
+ lA-^A-) is coupled to the quantum component of the source field, V. On the other hand, 
the classical source component, V', is nothing but an external physical scalar potential, the same at 
the two branches. 

Notation may be substantially compactified by introducing two vertex y-matrices: 











(J 






i) 



f'^[, \ ], r^[ \ oj. (113) 

With the help of these definitions, the source action (II 12b may be written as 

dr df 2 [r''hY!t<lJt + V^^arli/^h] = Tt[^V% (1 14) 

a.b=l 

where we have introduced Keldysh doublet *f and matrix V, defined as 
where a - (cl, q). 

In a similar way one may introduce external vector potential into the formalism. The corre- 
sponding part of the actiorlfls'A = /^r j[l df A(r, f)j(r, f) represents the coupling between A(r, t) and 
the fermion current j(r, t) - -^.[i^ir, r)5ri^(r, t) - drif/ix, t)iff{r, t)]. By splitting dt into forward 
and backward parts, performing Keldysh rotation, one finds by analogy with the scalar potential 
case (11121 ) that 

= Tr{#Avf f } , A = A^r = ( a'' ) " ^^^^^ 

We have linearized the fermionic dispersion relation near the Fermi energy and employed that -idr ~ 
Pf and \f = Pf /m. 

Let us now define the generating function as 

Z[V'',r'] = {expiiSv)) , (117) 

where the angular brackets denote the functional integration over the Grassmann fields and i/f 
with the weight exp(;5'o), specified by the fermionic action dllOl ). In the absence of the quantum 
component, V = 0, the source field is the same at both branches of the time contour Therefore, the 
evolution along the contour brings the system back to its exact original state. Thus, one expects that 
the classical component alone does not change the fundamental normalization, Z - 1 . As a result, 

Z[V''',0] = l, (118) 

as we already discussed in Section |2] see (|35] |. Indeed, one may verify this statement explicitly 
by expanding the partition function ( 11171 ) in powers of V ' and employing the Wick theorem. For 
example, in the first order, one finds Z[y'^', 0] = 1+J df Tr[7'^'G(f , f)] = 1 , where one uses that y'^' - 1 
along with ( 1105b . It is straightforward to see that for exactly the same reason all higher-order terms 
in V^' vanish as well. 

A lesson from ( 11181 ) is that one necessarily has to introduce quantum sources (which change sign 
between the forward and the backward branches of the contour). The presence of such source fields 
explicitly violates causality, and thus changes the generating function. On the other hand, these 



The vector source A(r, ;) that we are using here differs from the actual vector potential by the factor of e/c. However, we 
refer to it as the vector potential and restore electron charge in the final expressions. 
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fields usually do not have a physical meaning and play only an auxiliary role. In most cases one 
uses them only to generate observables by an appropriate differentiation. Indeed, as was mentioned 
above, the physical density is coupled to the quantum component of the source. In the end, one 
takes the quantum sources to be zero, restoring the causality of the action. Note that the classical 
component, V', does not have to be taken to zero. 

Let us see how it works. Suppose we are interested in the average fermion density g at time t in 
the presence of a certain physical scalar potential V^'it). According to ( 1112b and ( 11171 ) it is given by 

g(x; V'') = -- — ^ Z[V'', r']\ , (119) 

where x - (r, t). The problem is simplified if the external field, V"', is weak in some sense. One may 
then restrict oneself to the linear response, by defining the susceptibiUty 



n*(jc,jc') = ; g{x; V') 



■ ^2 2^ycl^ yq-^ 



(120) 

V'i=Vi=0 

We add the subscript R anticipating on the physical ground that the response function must be re- 
tarded (causality). We demonstrate it momentarily. First, let us introduce the polarization matrix 
as 

ti'^f^u x') = -i ^''"^['^ - 1 n^(^'-^') ] (121) 
II (x,x)- 2 svHx')sv^'{x) y^,~\n>'(x,x') n^(x,x')j- 

Owing to the fundamental normalization, dllSb . the logarithm is redundant for the R and the A 
components and therefore the two definitions (1120b and (1121b are not in contradiction. The fact 
that n''''' = is obvious from dl 18b . To evaluate the polarization matrix, 11, consider the Gaussian 
action dUOb . Adding the source term (II 141) . one finds; 5o + 5v - J dx *1'[G"' + Vy"]^. Integrating 
out the fermion fields if/, if/ according to the rules of fermionic Gaussian integration. Appendix lAl 
one obtains 

Z[y"', y«] = Det [;G"' + iV"j"] = Det \i + G VV] = exp (Trln[i + G Vy"]] , (122) 

Tr{po) ^ J L J I I 

where ( |97] i has been used. Since Z[0] = 1, the normalization is exactly right. One may now expand 
ln[l + G y"y"] to the second order in V". As a result, one finds for the polarization matrix 

tP^ix,x') = TT{y''G(x,x')fG(x',x)] , (123) 

which has a transparent diagrammatic representation, see Figure [T] 

Substituting the explicit form of the gamma matrices, (II 13b . and the Green's functions ( 1104b . 
one obtains for the response and the correlation components 

n«™,„', = -i |C».,„OC'(,',„ .C'„,,',C««„',,)] , (124a, 

nV.', = 4|c'(,.,')G'(,',„ + GV.')G'(,'.„ + GV..',GV.,)|. (124b) 

From the first line it is obvious that n^^'^^(x, x') is indeed a lower (upper) triangular matrix in the 
time domain, justifying their superscripts. Moreover, from the symmetry properties of the fermionic 
Green's functions one finds: 11* = [H'^] ' and Ft'*^ - -[Ft^] ' . As a result, the polarization matrix, H, 
possesses all the symmetry properties of the bosonic self-energy E, see ( |47] ). 

Equation ( 11241 ) for 11* constitutes the Kubo formula lfT2l l40l for the density-density response 
function. In equilibrium it may be derived using the Matsubara technique. The Matsubara routine 
involves the analytical continuation from discrete imaginary frequency (d,„ - InimT to the real fre- 
quency o). This procedure may prove to be cumbersome in specific applications. The purpose of 
the above discussion is to demonstrate how the linear response problems may be compactly formu- 
lated in the Keldysh language. The latter allows to circumvent the analytical continuation and yields 
results directly in the real frequency domain. 



32 




Figure 7: Polarization operator fY'^(x,x'): each solid line stands for the fermion matrix Green 
function ( 1104b . wavy lines represent external classical or quantum potentials y^'**', and x = (r, t). 
The loop diagram is a graphic representation of the trace in ( 1123b . 



5.5 Applications I: Quantum transport 
5.5.1 Landauer formula 

Let us illustrate how Keldysh technique can be applied to calculate Landauer conductance pTI of a 
quantum point contact (QPC). For that purpose consider quasi-one-dimensional adiabatic constric- 
tion connected to two reservoirs, to be referred to as left (L) and right (R). The distribution functions 
of electrons in the reservoirs are Fermi distributions ni(i{){ek) - [exp[(ejr - iihk))IT] + 1] with 
electrochemical potentials shifted by the voltage /i^, - /ig - eV. Within QPC electron motion is 
separable into transverse and longitudinal components. Owing to the confinement transverse motion 
is quantized and we assign quantum number « to label transverse conduction channels with (p„{v^^ 
being corresponding transversal wave functions. The longitudinal motion is described in terms of 
the extended scattering states, i.e. normalized electron plane waves incident from the left 

and the right 

onto mesoscopic scattering region (Figure |8]i. Here k{e) is the electron wave vector and t„(^) and 
r„{k) are channel specific transmission and reflection amplitudes. Second quantized electron field 
operator is introduced in the standard way 

❖(r, t)^Yj [^«(^' ^)««(^' + ')««(^' ' (127) 

nk 

where if/„ (k, f) are fermion destruction operators in the left and right reservoirs, respectively. For 
the future use we define also current operator 

/(X, M,t,-A;"(^, o-a' (^', , (128) 

nk,n'k' 

with the matrix elements 

Mf„,ix-k,k') ^^Jdr^ [u:"(k,r)dy„,(k',r) - [d,u:^{k,r)]4.{k' ,r)\ , a^L,R, (129) 
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Figure 8: Two terminal scattering problem from the quantum point contact. 

which are constructed from the scattering states ( I125l l-( ll26l l. Based on the orthogonality condition 
J drj^0„(rj^)0*,(rx) - 6nn', direct calculation of M„„>{x\ k, k') for x > gives 

^ ' '~ ' "" \ K{k)t„{k') vl{k)v„{k') 

where Vf — kf/m is Fermi velocity. For x < the expression for M is similar and different from ( 11301) 
by an overall sign and complex conjugation. The second approximate relation on the right-hand side 
is written for the case when the transmission amplitudes depend weakly on the wavenumber k on 
the scale dictated by temperature or the applied bias, and thus their momentum dependence may be 
disregarded. 

One can set up now the partition function for this transport problem as 

Z[A] = J D[#]exp|i#[G-' +AM]>f} , (131) 

here *F = (i^^, i^*), G - diag{G£, Gr) is 4 x 4 Green's function matrix, whereas is 2 x 2 matrix in 
the Keldysh space, and A is auxiliary vector potential, c.f. ( 1116b . Since the functional integral over 
fermionic fields in (I131I I is quadratic, one finds upon Gaussian integration 

In Z[A] = Tr In [i + GAM] . (132) 

In analogy with (11 19t the average current is generated from Z[A] via its functional differentiation 
with respect to the quantum component of the vector potential (/) - -{i /2)6lnZ[A] / 6A''(t)\Ai=Q- By 
expanding trace of the logarithm to the linear order in A, as Trln[l + GAM] x Tr[GAM], one finds 
for the current 

(133) 

'Equation )130) is obtained as a result of certain approximations. The exact expression for the cuiTent matrix explicitly 
depends on coordinate x. There are two types of terms. The first depends on x as exp(±i(<: + k')x) k exp(±2i<:;rA:), where kf 
is Fermi momentum, it represents Friedel oscillations. Their contribution to the cun'ent is small as (k - k')/kf « 1, and thus 
neglected. The second type of terms contains exp(±((<: - k')x) a 1, since \k - k'\ ~ L^^ <s .v"' , where Lj = Vf IT is ballistic 
thermal length, and the coordinate x is confined by the sample size L -i: Lj. See corresponding discussions in 1421 . 
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where we used Keldysh trace 

Tr{G„f = G,f (f, f , fc) = J^Ga > k) , 

and introduced QPC transmission probability T„{ek) - \i„{ky^. The last step is to take Keldysh 
component ofthe Green's function G,f(e, A:) = -27r;(5(e-eA+A'fl)[l-2«/r(e)], with - Vfk [see dlOTI il. 
and to perform momentum integration which is straightforward owing to the delta function in G^. 
The result is 

(/> = 1^ X / r„(e)[«L(e) - «R(e)] . (134) 
For a small temperature and applied voltage ( 11341 ) gives a conductance (/) - gV, where 

and all transmissions are taken at the Fermi energy T,, = T„(£f) (note that we restored Planck 
constant ti in the final expression for the conductance). Equation (1135b is known as a multi-channel 
Landauer formula (see Il43ll44ll for detailed reviews on this subject). 

5.5.2 Shot noise 

Based on the previous example we can make one step forward and calculate the second moment of 
the current fluctuations, so-called noise power, defined as the Fourier transform of current correla- 
tions 

S{aj,V) = J dte'"\6I{t)6i(0) + 6I{0)6i(t)), 61{t) ^ I{t) - (I) . (136) 



Within Keldysh technique this correlator may be deduced from Z[A] (see ( |132t ). Indeed, one needs 
now to expand trace of the logarithm in ( |132t to the second order in auxiliary vector potential A and 
diff'erentiate In Z[A] oc Tr[GAMGAM] twice over the quantum component, A^: 



1 6HnZ[A] 

V) 



2 6Ai{aj)5A'i{-aj) 



(137) 

A'i=Q 



This expression automatically gives properly symmetrized noise power ( 1136b . As a result of the 
diff'erentiation one finds 



nkk' ^ 

+T„R„Tr{Gde^)rGR{e-)r} + r„7?„Tr{Gfi(e+)y«Gi(e_)r«) + r2Tr{G«(e+)r«Gfl(e_)r«)l , (138) 

where we already calculated partial trace over the left/right subspace, assuming that transmissions 
are energy independent, and used notations e± - e + cojl and R„ - \ - T^. Calculation of Keldysh 
traces requires ( |104b and dl 131 ) and gives 

Tr{GarGbr] = G^Gf + G^„G^, + G^Gf . (139) 

The remaining step is the momentum integration. One uses Gf ^\e, k) - (e - Vfk + + ;0)"' and 
Gf,{e,k) = -Inidie - vpk + //„)[! - Inpie)] from ([l03, and finds that Y.kk' J deTriG afGhy^} = 
vf /de [1 - (1 - 2«fl)(l - Irib)]. As a result, the final expression for the noise power obtained by 
Lesovik ||45]| reads as 



S(oj, y)^-^Yj / '^^ [T^nBLLie) + T„R„Blr{€) + T„R„BRL(e) + TlBRR(e)\ , (140) 
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where statistical factors are Bab(e) - ««(£+)[! - «/>(£-)] + «i(e-)[l - na(e+)] and we again restored H 
in the end. Despite its complicated appearance, e integration in (11401 ) can be performed in the closed 
form0 



2 

S{cj, V)^^Yj [^«'^''°*(^) + ^"^^ " ^"^^"^^ + w)coth( 



eV + a>\ 
IT ] 



(141) 



There are two limiting cases of interest, which can be easily extracted from (11411 ). The first one 
corresponds to the thermally equilibrium current fluctuations, V — » 0. In this case 

>S(w,0) = 2g^^coth(^) , (142) 

where we used ( I135I I for conductance g. This result is nothing but familiar fluctuation-dissipation 
relation for the current fluctuations. Note, that despite the complicated dependence on transmission 
amplitudes in ( |14Q| i the equilibrium noise power (I142l i is written in terms of conductance ( 11351 ) only. 
The other limiting case is fully non-equilibrium noise at zero temperature T ^ Q and a finite bias 
V. For such a case one finds from (1141b for the excess part of the noise 

S{u, V) - S(u, 0) = ^{\eV + + \eV 2|M) ^ T„{1 - T„) , (143) 

n 

which is called the shot noise. An important observation here is that in contrast to equilibrium 
noise (1142b . shot noise can not be written solely in terms of the conductance g. Only for the case 
of tunnel junction, where all transmissions are small, T„ «; 1, Equation ( 11431 ) reduces to S(0, V) = 
2eVg - 2e{I), which is known as Schottky formula (for a review of shot noise in various systems 
see e.g. Il46l|l7l|48l). 



5.5.3 Coulomb drag 

Drag effect proposed by Pogrebinskii 1491 and Price ifSOl by now is one of the standard ways to 
access and measure electron-electron scattering. In bulk two-dimensional systems (two parallel 
two-dimensional electron gases, separated by an insulator) the drag eff'ect is well established ex- 
perimentafly l5l]|52l|53l|5l|55l and studied theoretically l56l 1221 ES |59l . Recendy a number of 
experiments were performed to study Coulomb drag in quantum confined geometries such as quan- 
tum wires 1601 [61] |62] |63l, quantum dots l64l|65l or QPCs l66l . In these systems a source-drain 
voltage V is applied to generate current in the drive circuit while an induced current (or voltage) is 
measured in the drag circuit. Such a drag current is a function of the drive voltage V as well as gate 
voltages, Vg, which control transmission of one or both circuits. Figure|9h shows an example of such 
a setup, where both drive and drag circuits are represented by two QPCs. 

The Keldsyh technique is an efficient way to tackle the drag problem both in linear response 
regime and away from the equilibrium, when a relatively large bias is applied to the drive circuit. 
Within each QPC electrons are assumed to be non-interacting and their motion is separated into 
quantized-transversal, and extended-longitudinal, see Section 15.5.11 The action describing non- 
interacting point contacts is 

!5qpc = !Tr{#G-''P), (144) 

where ^ - (i/^^^^, i/'*^) and G = (Jj/diaglG/,, Gr}. Index j - 1,2 labels QPCi(2) respectively, n is the 
transverse channel index within each QPC, and Gl(r) is a 2 x 2 Keldysh matrix ( 11041) . 
The interaction term between the two QPC is 

i5in, = 2 dtdt' Il{t)Kf^(t-t')I^Jt'), (145) 



'"Deriving 11411 one writes statistical factors as Bab(^) = 5 [1 -tanh[(e+ -fia)/2T] tanh[(e_-/ii,)/2r]] and uses the integral 
J^" d.v [1 - tanh(.v + i/)tanh(.v - y)] = 4i/coth(2y). 
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Gate voltage - Vg 



Figure 9: a) Two coupled QPCs and surrounding electric circuitry. The Coulomb coupling is due 
to mutual capacitances Cc- Gate voltage Vg control transmission of, e.g., drive QPC. b) Schematic 
representation of conductance of the drive QPC along with the drag current as a function of the gate 
voltage. 



where IjR(L)it) are current operators, on the right (left) of QPCy, coupled by the kernel Ko/,(f - t'), 
which encodes electromagnetic environment of the circuit. The retarded and advanced components 
of the interaction kernel are related to the trans-impedance matrix K*^'^*(a;) = Z^^^\(jj) I {a> ± iO). 

The latter is defined as Z**'^'(w) = d<b „(+(!>) I dlh{+ui), where the corresponding local fluctuating 
currents !„ and voltages (Dq are indicated in Figure |9^. The Keldysh component of the interaction 
kernel is dictated by the fluctuation-dissipation theorem: K^^(w) = [K*^(w) - K^^(w)] coth(w/2r), 
i.e. we assume that the surrounding electric environment is close to equilibrium. Finally the current 
operators are given by ( 11281 ) and ( I130l l. 

The drag current is found by averaging h over the fermionic degrees of freedom 

/o = J D[#] Tr [0-2^(^2] exp (/5qpc[0"A] + iSU^^\) ■ (146) 
Expanding the exponent to the second order in the interaction term S int, one obtains 

/d = ^ J D[#]Tr[(A2MiA2]Tr[/iK/2]Tr[/iK/2]exp(/5QPc[#]). (147) 

The remaining Gaussian integral over the fermionic fields is calculated using the Wick's theorem. 
One employs expression (1128b for the current operators with the M-matrix given by ( 1130b and takes 
into the account all possible Wick's contractions between the i/r-fields. The latter are given by the 
Green's functions ( 1104b . This way one finds for the drag current 

/J 
^Tr [t{co)Sx (w, y)Z(-w)f2(w)] . (148) 

The drive circuit is characterized by the excess part Sf'{a), V) = Sah{(^, V)-Sab{oJ, 0) of the current- 
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current correlation matrix Sahico, V) - j dt e'"'{{SIa{t)6Ib(0) + 6Ib(0)6Ia{t))), given by, e.g., 
Srr(oj, V) = ^ y r de[BLL(e)\t^{e^f\t^(e-f + Bz.«(e)|t^(e+)|V*(e-)|' 

+Bfiz,(6)|r*(e+)ht^(6-)|2 + B««(6)[l - rf (e+)r«(e_)][l - r«(e+)rf (e_)]] , (149) 

where e± = e + w/2 , t^***(e±) = t^***(e± + eV^R)) and r^'^''\e^) = r,^'''\e± + eV^R)), while Rq = 
2nh/e- is quantum resistance, and statistical occupation form-factors Babie) are given by ( I140l i. 
Here SRR{a>, V) generahzes ( I140l l to the case of energy dependent transmissions B2l . Expressions 
for other components of the noise matrix Sll, Slr, and >Sr£ are similar, see Refs. Il42ll68l . 

The drag circuit in (1148b is characterized by the rectification coefficient r2(w) - r2{a>)g~ of ac 
voltage fluctuations applied to the (near-equilibrium) dra g Q PC? , where is the third Pauli matrix 
acting in the left-right subspace. Rectification is given bvr'l 

riico) = ;^ J ["^(^-^ - nF(e+)][|t„(e+)|' - \t„{e-)\^] . (150) 

Characteristics of the QPC2 enter through its energy-dependent transmission probabilities |t„(e)P. 
This expression admits a transparent interpretation: potential fluctuations with frequency at, say 
on the left of the QPC, create electron-hole pairs with energies e± on the branch of right moving 
particles. Consequently the electrons can pass through the QPC with the probability |t„(e+)p, while 
the holes with the probability |t„(e-)|~. The difference between the two gives the dc current flowing 
across the QPC. Note that the energy dependence of the transmission probabilities in the drag QPC is 
crucial in order to have the asymmetry between electrons and holes, and thus non-zero rectification 
r2(w). At the diagrammatic level ( 1148b has transparent representation shown in Figure [TOl 

Focusing on a single partially open channel in a smooth QPC, one may think of the potential 
barrier across it as being practically parabolic. In such a case its transmission probability is given by 

mf = [expiieVg - 6)/Aj} + l)"' , (151) 

where A; is an energy scale associated with the curvature of the parabolic barrier in QPC; and gate 
voltage Vg shifts the top of the barrier relative to the Fermi energy. This form of transmission was 
used to explain QPC conductance quantization [69] and it turns out to be useful in application to the 
Coulomb drag problem. Inserting ( 11511 ) into Eq. ( 11501 ) and carrying out the energy integration, one 
finds 



Ytfcj) In 

Ro 



sinh^(w/2A2) 



1 + y-^ — (152) 

cosh2(ey,/2A2)j 

for r « A2. In the other limit, T » A2, one should replace A2 — > T in ( 1152b . Note that for small 
frequency w <*: A2 one has r2 ~ lt, thus making the integral in ( 1148b convergent in w — > region. 

Linear drag regime. For small applied voltages V one expects the response current lo to be 
linear in V. Expanding S\{(jl>, V) to the linear order in V, one finds that only diagonal components of 
the current-current correlation matrix contribute to the linear response and as a result. 



oto 



coth — 
2T 



ri(w)^, + 0(n, (153) 



"in teiTns of the Keldysh matrices the rectification coefficient is given by the following trace r2ici>) = 
TrlGf MGy^' MGf-'' M]. Finding r2((<;) in the form of fBol one uses Keldysh trace Tr [Gy«Gf ''Gf '''] = 
i;+[G'^(e)G'^(e±tj)G'^(6) + G'^(6)G'^(6±w)G'^(e) + G'^(e)G'*{6±w)G'^(6)]. To simplify this expression further one 
should decompose each Keldysh component of the Green's function using fluctuation-dissipation relation G'^(e) = [G'*(e) - 
G'^Ce)]!! - 2n{e)] and keep in the resulting expression only those terms, which have a proper causality, i.e. combinations 
having three Green's functions of the same kind, like G'*G'*G'* and G'^G'^G", do not contribute. In this way, one finds for 
the Keldysh trace Tr [Gy''G-y'^'Gy'''j oc [jif - nf {£+)]. Remaining trace in the left-right subspace over the current vertex 
matrices M reduces to the transmission probabilities at shifted energies, namely Tr[MMM] k |t„(6+)P - |t„(6_)P, leading 

to (Bo). 
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Figure 10: Drag current lo in the second order in inter-circuit interactions K - Z/w (wavy lines). 
The drag circuit is represented by triangular rectification vertex FaCw), while the drive circuit by the 
non-equilibrium current-current correlator >Si((jj, V) (loop). 

where T\{<jj) is obtained from ( |150t by substituting transmission probabilities of QPC2, by that of 
QPCi. Inserting ( 1153b into (1148b one finds 



where dimensionless interaction kernel a+{u) is expressed through the trans-impedance matrix as 



Equation (1154b has the same general structure as the one for the drag current in bulk two-dimensional 
systems ||58ll59l . Being symmetric with respect 1 <-> 2 permutation, it satisfies Onsager relation for 
the linear response coeflicient. Performing remaining frequency integration in ( |154b . it is sufficient to 
take the interaction kernel at zero frequency. Indeed, frequency scale at which a+(w) changes is set 
by inverse /?C-time of the circuit. If load impedance of the drag circuit is large compared to that of 
the drive circuit Zi <*: Z2 <sc Rq, which is the case for most experiments, and the mutual capacitance 
of the two circuits is small Cc <*: Cr^l,s, see Figure |9^, one finds t^^ = (ZiC,) ' » T . Since Iq 
in ( 11541 ) is determined by w < T, it is justified to approximat^H a+{<jS) ^ a+{Q). Substituting ( 11521 ) 
into ( 11541 ). one finds for, e.g., low-temperature regime T <sc Ai_2 




°^+(") ^ T^Tr[Z(w)^,Z(w)5^;] . 




V a +{0)71- 1 



(155) 



Rq 6 A1A2 cosh^(ey3/2Ai) ' 



■For the circuit sliown in tlie Figure|2]one finds for tlie low-frequency limit of the trans-impedance kernel 
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where we assumed that the gate vohage of QPC2 is tuned to adjust the top of its barrier with the 
Fermi energy and wrote lo as a function of the gate voltage in QPCi. The resulting expression 
exhibits a peak at V^, = similar to that depicted in Fig.|9}5. This expression describes rectification 
of near-equilibrium thermal fluctuations (hence the factor T-), which is due to the electron-hole 
asymmetry (hence, non-monotonous dependence on Vg). 

Nonlinear regime. At larger drive voltages drag current ceases to be linear in V. Furthermore, 
contrary to the linear response case, S\{(jj, V) does not require energy dependence of the transmis- 
sion probabilities and could be evaluated for energy independent |t„|~ (this is a fare assumption for 
T,eV <K Ai). Assuming in addition T «: eV, one finds Sf{u,V) = [^^^(w, V) - .S„,,(w, 0)]s^o, 
where S\{<jd, V) is given by ( 1143b (recall that T„ = |t„|^). Inserting it into ( 1148b . after the frequency 
integration bounded by the voltage, one finds for the drag current 

Id = t^-(O) y r„(l - T„) . (156) 

Here again we assumed that the detector QPC2 is tuned to the transition between the plateaus. We 
also assumed eV ^ (ZiCj)"' to substitute 

^-('^^ ^ T7;2"Tr[Z(w)^oZ(w)?;J 

by its dc value, a^iQ). One should notice that while a+ > 0, the sign of is arbitrary, since 
a- oc Cj^ - C^, see Figure and Note 12. For a completely symmetric circuit Qf_ = 0, while for 
extremely asymmetric one Iff-I x a +12. Although we presented derivation of ( 1156b for T <«c eV, one 
may show that it remains valid at any temperature as long as T <«; min{Ai, (ZiC.,)"')- 

Equation ( 1156b shows that the drag current is the result of the rectification of the quantum shot 
noise and is hence proportional to the Fano factor pSi of the drive circuit. It exhibits the generic 
behavior depicted in Figure |9}5, but the reason is rather different from the similar behavior in the 
hnear regime. The direction of the non-linear drag current is determined by the inversion asymmetry 
of the circuit (through the sign of a J) rather than the direction of the drive current. As a result, for a 
certain polarity of the drive voltage, the drag current appears to be negative. Finally, assuming that 
for a generic circuit a+ ~ and comparing ( 1155b and ( 11561 ) one concludes that the transition from 
the linear to the non-linear regime takes place at V a: V* with eV* - T^/Ai <K T, for T <K Ai . In the 
opposite limit, T > Ai, the crossover voltage is given by the temperature eV* = T. Further details 
and discussions can be found in 1681 . 



6 Disordered fermionic systems 

One is often interested in calculating, say, density-density or current-current response functions, in 
the presence of static (quenched) space-dependent disorder potential f/dis(r). Moreover, one wants 
to know their averages taken over an ensemble of realizations of t/dis(r), since the exact form of 
the disorder potential is, in general, not known. The response function in the Keldysh formula- 
tion, may be defined as variation of the generating function and not the logarithm of the generating 
function. More precisely, the two definitions with, and without the logarithm coincide owing to the 
fundamental normalization Z - \. This is not the case in the equilibrium formalism, where the 
presence of the logarithm (leading to the factor Z ' after differentiation) is unavoidable in order to 
have the correct normalization. Such a disorder-dependent factor Z ' - Z '[t/dis] formidably com- 
plicates the averaging over f/dis- Two techniques were invented to perform the averaging: the replica 
trick II25I l26l l27l l28l I29I and the supersymmetry 13111321 . The first one utilizes the observation that 
InZ - lim„^o(Z" - to perform calculations for an integer number, «, of replicas of the same 
system and take n — > in the end of the calculations. The second one is based on the fact that 
Z"' of the non-interacting fermionic system equals to Z of a bosonic system in the same random 
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potential. One thus introduces an additional bosonic replica of the fermionic system at hand. The 
Keldysh formalism provides an alternative to these two methods ensuring that Z = 1 by construc- 
tion ll2n I22I l23l . The purpose of this section is to show how the effective field theory of disordered 
electron gas, known as the non-linear cr-model (NLSM), is constructed within Keldysh formalism. 



6.1 Disorder averaging 

We add the disorder dependent term to the fermionic action S dis [if', tj/j-J^dtJ dr t/dis(r)i^(r, t)if/(r, t), 
where t/dis(r) is a static scalar potential, created by a random configuration of impurities. It is usu- 
ally reasonable to assume that impurities are short-ranged and distributed uniformly over the system, 
thus having the correlation function of the form (t/dis(r)t/dis(r')) ~ <5(r - r')- Assuming, in addition, 
Gaussian distribution of the impurity potential, one ends up with the disorder averaging performed 
with the help of the following functional integral: 



<. . .>di,s 



J D[t/dis]...exp|-7ryTei J dr C/2,(r)| , (157) 



where the disorder strength is characterized by the elastic mean free time Tei, and v is the electronic 
density of states at the Fermi energy. Since the disorder potential possesses only the classical com- 
ponent, it is exactly the same on both branches of the Keldysh contour. Thus, it is coupled only to 
- 1 vertex matrix. Next, we perform the Gaussian integration over f/dis of the disorder-dependent 
term of the partition function (at this step we crucially use the absence of the normalization factor) 
and find 

J D[C/diJ exp J- J dr nvr^.Ul^ir) - /f/di,(r) J dr^«(r, t)f^,ij,\r, t) j 

= exp(-^^ J'^'" IT ^^Jt'^'^'"' ^'^j) ' ^^^^^ 

where a,b - 1,2, and summations over all repeated indices are assumed. One can rearrange 
[(/'"(r, f)iA°(r, t)][\l/^(r, t')ip\r, t')] = -[if>"{r, Oi/z^r, t')][il/\r, t')\p"{r, t)] in the exponent on the right- 
hand side of the last equation (the minus sign originates from anti-commuting property of the Grass- 
mann numbers) and then use Hubbard-Stratonovich matrix-valued field, Q - Q"l^{r) to decouple 
(time non-local) four-fermion term asF^ 

exp(^^ Jdr ^ dtdt'[;fr\r, f)/(r, t')M\r, t')if,yr, f)]J 

= J D[Q] exp {-^Tr{Q~}-^^j dr jj dfdr'0;(r)iAV, O-AV, ?)) • (159) 

Here we have introduced that the trace of the implies summation over the matrix indices as well 
as time and spatial integrations 

dr dfdr' Q'lfimnir). (160) 

a,b=\ 

Now the averaged action is quadratic in the Grassmann variables S [*iF, Q] — Tr{*iF[G"' + ^^Q\^], 
and they may be integrated out explicitly, leading to the determinant of the corresponding quadratic 



"since we do not keep track of the time-reversal symmetry, i.e. the fact that the Hamiltonian is a real operator, the 
following considerations are restricted to the case, where the time-reversal invariance is broken by, e.g., an external magnetic 
field (complex Hermitian Hamiltonian). This is the so-called unitary NLSM. The orthogonal NLSM, i.e. the one where the 
time-reversal symmetry is restored is considered in Section[8] devoted to disordered superconductors. 
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form: G ' + 577 2- of the matrices here should be understood as having 2x2 Keldysh struc- 
ture along with the N x N structure in the discrete time. One thus finds for the disorder averaged 
generating function Z. = (Z)dis: 



z 



j T)[Q] exp(/5[e]), 



iS[Q] = -—Tr{Q'-} + Tr\n 

4Tel 



2Tel 



(161) 



As a result, one has traded the initial functional integral over the static field C/dis(r) for the functional 
integral over the dynamic matrix field Q„'{r). At a first glance, it does not strike as a terribly bright 
idea. Nevertheless, there is a great simplification hidden in this procedure. The point is that the 
disorder potential, being ^-correlated, is a rapidly oscillating function. On the other hand, as shown 
below, the Q matrix field is a slow (both in space and time) function. Thus, it represents true 
macroscopic (or hydrodynamic) degrees of freedom of the system, which are diffusively propagating 
modes. 



6.2 Non-linear cr-model 

To proceed we look for stationary configurations of the action S[Q] in ( |161t . Taking the variation 
over Qi,'{r), one obtains the saddle point equation 



e(r)=-G-' + — e , (162) 
nv \ 2Tel 

where 2„,(r) denotes a stationary configuration of the fluctuating field Qn'ir)- The strategy is to find 

first a spatially uniform and time-translationally invariant solution Q of ( 1162b and then consider 
space- and time-dependent deviations from such a solution. This strategy is adopted from the theory 
of magnetic systems, where one first finds a uniform static magnetized configurations and then treats 
spin-waves as smooth perturbations on top of such a static uniform solution. From the structure 
of ( II62I 1 one expects that the stationary configuration Q possesses the same form as the fermionic 
self-energy ( 1106b (more accurately, one expects that among possible stationary configurations there 
is a classical configuration that admits the causality structure ( 11061 )). One looks, therefore, for a 
solution of ( 1162b in the form of the matrix 

/ A* , A^ , \ 

G,_„=A,-.=( - ^r). (163) 

Substituting this expression into (1162b . which in the energy/momentum representation reads as A^ = 

^ 2/7 " f/j + T77^e) with €p = 12m - ep, one finds 

= - E - ±1 ' (164) 

where one adopts the convention Yjp - ■ ■ ^ ^/ de^. The signs on the righ-hand side are chosen 
so as to respect causality: the retarded (advanced) Green's function is analytic in the entire upper 
(lower) half-plane of complex energy e. One has also assumed that I/tbi <k ep to extend the energy 
integration to minus infinity, while using constant density of states v. The Keldysh component, as 
always, may be parameterized through a Hermitian distribution function: A^ = Af^oF-FoK^ , where 
the distribution function F is not fixed by the saddle point equation (1162b and must be determined 
through the boundary conditions. In equilibrium, however, F is nothing but the thermal fermionic 
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A. = | A], (165) 



distribution function F^'^ = tanh thus Af = (Af - A^)F^'' = IFf. Finally we have for the 
stationary Q matrix configuration 

If 2F, 
-1^ 

where we have introduced the retarded and advanced unit matrices to remind about causality struc- 
ture and the superscript "e^" in the distribution F was suppressed for brevity. Transforming back to 
the time representation, one finds aJ^*^' - ±6(t - f' + 0), where +0 indicates that 5-function is shifted 
below (above) the main diagonal, t - f . As a result, Tr{A) - and 5 [A] = 0, as it should be, of 
course, for any purely classical field configuration ( 11631 ). One should note, however, that this partic- 
ular form of the saddle point solution ( 1165b is a result of the approximation that the single-particle 
density of states v is independent of energy. In general, it does depend on e and thus retarded (ad- 
vanced) components of A^ are analytic functions of energy in the upper (lower) half-plane, which 
do depend on energy on the scale of order of the Fermi energy €f. Therefore, the infinitesimally 
shifted 5-functions in aJ^'^' = +5(t - f' + 0) should be understood as 5,+o - f±(t)6(±t), where 9{±t) 
is the Heaviside step function, and f±(f) are functions that are highly peaked for \t\ < e^' and satisfy 
the normalization j^^"' dtf+it) = +1. This high-energy regularization is important to remember in 
calculations to avoid spurious unphysical constants. In particular, for this reasons \f ^,Mf,, - 0, and 
l'^_,,Mf,j = 0, where Mf^'^ is an arbitrary retarded (advanced) matrix in the time space. 

Now we are on a position to examine the fluctuations around the saddle point ( 1165b . The fluc- 
tuations of Q fall into two general classes; (i) massive, with the mass oc v/rgi and (ii) massless, i.e. 
such that the action depends only on gradients or time derivatives of these degrees of freedom. The 
fluctuations along the massive modes can be integrated out in the Gaussian approximation and lead 
to insignificant renormalization of the parameters in the action. The massless, or Goldstone, modes 
describe diffusive motion of the electrons. The fluctuations of Q matrix along these massless modes 
are not small and should be parameterized by the matrices satisfying a certain non-linear constraint. 
To identify the relevant Goldstone modes consider the first term in the action S[Q\ of (1161b . The 
stationary configuration given by (1165b satisfies 

Q'-['l = (166) 

Note that Tr{2^} = Tr{l*) + Trjl'*) = 0, owing to the definition of the retarded/advanced unit 
matrices. The fluctuations of Q which do not satisfy (1166b are massive. The class of Q matrix 
configurations, that obeys the constraint (1166b . is generated by rotations of the stationary matrix A^ 
and may be parameterized as follows 

Q^n-^ okoii. (167) 

The specific form of 'R is not important at the moment and will be chosen later. The massless modes, 
or spin waves, if one adopts magnetic analogy, which are associated with '^(f'(r) are slow functions 
of f -H f' and r and their gradients are small. Our goal now is to derive an action for soft-mode Q 
field configurations given by ( 11661 ) and ( 1167b . 

To this end, one substitutes (1167b into (1161b and cyclically permutes matrices under the trace. 
This way one arrives al'RoG-^ o-Rr^ = (J-i +'ko[G-^ ° ■^-'] = G'^ +i'kd,'k-^ +m\Fdr'Rr\ where 
one has linearized the dispersion relation near the Fermi surface = p^/2m — ef ~ Vf p — > -iypdr- 
As a result, the desired action has the form 

iS [Q\ = Trln [l + iQnd/R-^ + ignypdA'^] , (168) 

where we omit circular multiplication sign for brevity. Here Q is the impurity dressed Green's 
function matrix, defined through the Dyson equation (G ' + ■;^^K)0 - 1. For practical calculations 
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it is convenient to write Q in the form 



^"(? |a ) = ^^''[1+A] + ^^^[1-A], (169) 
with retarded, advanced and Keldysh components given by 

g^'^Kv, e)^[€-ep± //2Tei]-' , §"(9, e) = ^«(p, e)F, - F,g\p, e) . (170) 



One may now expand the logarithm in (1168b in gradients of the rotation matrices "R to the linear 
order in di'Rr^ and to the quadratic order in cJr'^ ' terms (contribution, linear in the spatial gradient, 
vanishes owing to the angular integration). As a result 

iS [Q] « iTrigm''] + ^Tr{^(^Vf5r^-')^(!Rvf5r5R-')l . (171) 

Since 2^ ^(p, e) = -invA^, which directly follows from the saddle point Equation ( 1162b . one finds 
for the df term in the action /Tr{^'^5,!R ' ) = 7TvTr{d,Q}. For the dr term, one finds -^nvOTiKdrQ)^}, 
where D - u^Tei/c/ is the diffusion constant and d is the spatial dimensionality. Indeed, for the prod- 
uct of the Green's functions one uses Yjp S'^iP^ f)vF^'^(p, e)\F = 2nvT^\v~pld - InvD, while the cor- 
responding R — R and A — A terms vanish upon performing integration. Employing then ( |169b . 
one arrives at Tr{[i+A](:R5r:^"i)[i-A]C^5rlR"')) = -\Tx{{dr{1i-^ kli)f] = -\Tx[{drQf]. Finally, 
one finds for the action of the soft-mode configurations lISTl l22l l23l 

iS [Q] = ~TT[D{dM - Ad,Q] . (172) 

Despite of its simple appearance, the action ( 11721 ) is highly non-linear owing to the constraint = 
1 . The theory specified by ( 1166b and ( 1172b is called the matrix non-linear cr-model. The name came 
from the theory of magnetism, where the unit-length vector cf{r), represents a local (classical) spin, 
that may rotate over the sphere cP- = I. 

One may now incorporate source terms Sy and Sa (see (II 12b and ( II 16b ) into the fermionic part 
of the action: 



' + 7r^e + V + V/.A 

2rei 



After Gaussian integration over and one finds for the source-fields-dependent partition func- 
tion, compare with (cf. ( 11611 )) 

Z[A,y] = jD[e]exp(/5[e,A,y]), 

(173) 



iS[Q,A, V] = -i^Trjg^) +Trln 



4Te 



G"' + —Q + V + \fA 



2Ter 



Expanding trace of the logarithm in gradients of Q with the help of ( 1167b . one assumes that source 
fields V and A are small in some sense and do not disturb the stationary configuration ( 1165b (see 
Section|7]for discussions of this point). Then, similarly to ( 11721 ). one finds from (11731 ) 

iS [Q, A, y] = y Tr{ y,T,y} - ^ Tr{D(4e)' -4d,Q + 4iVQ} , (174) 
where o", is the PauU matrix acting in the Keldysh space, and we have introduced covariant derivative 

drQ^drQ-i[kQ]. (175) 

A few comments are in order regarding ( 1174b . First, it is still restricted to the manifold of Q matrices 
satisfying Q~ = I. The second trace on the right-hand side of ( 1174b . containing Q, originates from 
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Yip ^fQ^^fG^ and Yp combinations in the expansion of the logarithm. On the other hand, 
the first term on the right-hand side of ( 1174b originates from Yp Q^Q^ and Yjp @^Q^ combinations. 
These terms should be retained since the matrix Y"{e - e')T" is not restricted to the 1 /xei shell near 
the Fermi energy. This is so, because the scalar potential shifts the entire electronic band and not 
only energy strip |e|, |e'| < 1 /rei. Thus, it is essential to follow the variations of the electron spectrum 
all the way down to the bottom of the band to respect the charge neutrality. To derive Triyo-vV) one 
has to employ the fact that for any physical fermionic distribution function F^^^oa — » +1. Equations 
(1174b and (1175b generalize an effective cr-model action given by (1172b . Additional technical details 
needed to derive (1174b from ( 1173b are provided in Appendix ICl 



6.3 Tunneling action 



Consider two metallic leads separated by a tunneling barrier, such that upon applying external volt- 
age a current may flow between them. In this case one has to add corresponding tunneling term to 
the Hamiltonian of the system 

r dr r dr'[r„,0'l(r)'AR(r') + 7';r''AR(r')'AL(r)], 

JrGl Jr'e« 

where i^l(r) is the electron annihilation operator to the left(right) from the tunneling barrier. The 
i^lj^^j is corresponding creation operator The T^^, and r*^, are tunneling matrix elements whose 
range is restricted to the vicinity of the junction, since the overlap of electron wave functions decay 
exponentially away from it. Tunneling Hamiltonian translates into the fermionic tunneling action 

/5 r = £ df jj drdr'[r„,0'L(r, f)iA«(r', + T*^^R{r\ f)<AL(r, 0] . 

Since 5 r is still quadratic in fermion fields, the Gaussian integration over them is straightforward, 
leading to the disorder averaged action in the form 



Deriving (1176b one has to introduce two Q matrices to decouple disorder mediated four-fermion 
term ( 1159b in each of the two leads independently. In doing so it was assumed for simplicity that 
both disordered samples are characterized by equal mean free times and bare electronic densities of 
states. Equation (1176b contains an additional 2x2 matrix structure in the space of left-right electronic 
subsystems, described by 2l(R), respectively. Note also that the tunneling matrix elements entering 
S \Ql, 2r] are unit matrices in the Keldysh subspace f„' - T„'0-q. 

Introducing the notation G^' = G^' + 277 one identically rewrites the last term of the action 
S [Ql, Qr] in dnlll as 



Trln 



ft 



= Trln 









G«' 



+ Trln 



1 





G«fi- 



GlT 




(177) 



Expanding now TrlnG^' in gradients of Qa matrix around the saddle point A,,, one obtains sigma 
model action ( 1172b . for each of the two leads independently. The coupling between them is described 
by the second term on the right-hand side of ( 1177b . which defines tunneling action 5 ^[Ql, Qr\. For a 
small transparency tunneling junction, one may expand trace of the logarithm to the leading (second) 
order in T and obtain 



/5r[eL,eR] =Trln 



1 





Gfif- 



GlT 




-TrlG^f G«f '"} + 



(178) 
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Employing the local nature of matrix elements Trr' and the fact that at the soft-mode manifold 
Qa - ;^Ga(r, r), see ( 1162b . one finds for the tunneling part of the action 

iSAQl, Qr] = ^Tx{QlQr] = -^Tx[{Ql - Qr?] . (179) 

4ge 8gQ 

Here we approximated the tunneling matrices as Trr' - T()6{r - r') and introduced the tunneling 
conductance gj - Ajre'^\T(!^v'^ , and the quantum conductance gg = e~ I{2jiK). The tunneling ac- 
tion ( |179l l is a generalization of the Tr{D(5r2)~) term of the NLSM action ( 11721 ) for the tunneling 
geometry. 

If the tunneling amplitudes T^' are not small one needs to keep higher orders in the expansion 
of the logarithm in ( |178l l. It is convenient to express products of the even number of the tunneling 
amplitudes Trr' through the transmission probabilities of individual transverse channels T„ (see, for 
example. Appendix C of lf73l ). With the help of ( 1162b . one may show that expansion of the logarithm 
in ( 11781 ) is order by order equivalent to the expansion of the following action llTOl ITTl l72l 



iSrlQL, Qr]^\Yj^^^A^-^{Ql- Qr) 



(180) 



If all transmissions are small, T„ <^ 1, one may expand dlSOb to the leading order in T,, and re- 
cover (1179b . identifying the tunneling conductance as gj - gg 2„ T„, c.f. ( 1135b . Equation dlSOl) 
goes beyond this limit and allows mesoscopic transport to be treat in arbitrary two-terminal geome- 
tries. Its generalization for multi-terminal case was also developed by Nazarov et. al. Il70ll74ll75l . 



6.4 Usadel equation 

Let us return to the action specified by (1172b . Our goal is to investigate the physical consequences 
of NLSM. As a first step, one needs to determine the most probable (stationary) configuration, 
2^ ,(r), on the soft-mode manifold (1166b . To this end, one parameterizes deviations from G^^,(r) as 
Q^'Rr^ oQo'R and chooses 'R - exp(^/2), where ^,/'(r) is the generator of rotations. Expanding 
to the first order in HV, one finds Q - Q— ['W ° Q]/2. One may now substitute such a Q matrix 

into the action ( 11721 ) and require that the terms linear in "W vanish. This leads to the saddle point 
equation for Q. For the first term in the curly brackets on the right-hand side of ( 11721 ) one obtains 

^Tri'WdrD [(drQ)Q- QdrQ]} = -Tri'WdrDiQdrQ)}, where one has employed OrQoQ+QodrQ = 

0, since 2=1. For the second term one finds Tr{'Wtt'idi + - Tr{^{(9,, Q}]. Demanding 

that the linear term in "W vanishes, one obtains 

driDQodrQ)-{d„Q}^0. (181) 

This is the Usadel equation 1761 for the stationary Q-matrix. If one looks for the solution of the 
Usadel equation in the subspace of "classical", having causality structure, configurations, then one 
takes Q - A, with as-yet unspecified distribution function F,i'{r). Therefore, in this case the Usadel 
equation is reduced to the single equation for the distribution function F„'{r). Substituting A from 
( 11651 ) into ( |181b and performing the Wigner transformation 

^'.'(r) -f^F, (r, T) e-^^^'-''\ r = ^ , (182) 

one obtains 

dr[D(r)drF,(r, r)] - 5,F,(r, r) = , (183) 

where we allowed for a (smooth) spatial dependence of the diffusion constant. This is the kinetic 
equation for the fermionic distribution function of the disordered system in the non-interacting limit. 
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which happens to be the diffusion equation. Note that it is the same equation for any energy e and 
different energies do not "talk" to each other, which is natural for the non-interacting system. In 
the presence of interactions, the equation acquires the collision integral on the right-hand side that 
mixes different energies between themselves. It is worth mentioning that elastic scattering does not 
show up in the collision integral. It was already fully taken into account in the derivation of the 
Usadel equation and went into the diffusion term. 

As an example, let us consider a disordered quasi-one-dimensional wire of length L, attached to 
two leads, kept at different voltages ITTTI . We look for the space dependent, stationary function Fdx) 
with X being coordinate along the wire, that satisfies D d^Fdx) - 0, supplemented by the boundary 
conditions F^{x = 0) = Fiie) and F^{x - L) - Fr{€), where Fk(£)(6) are the distribution functions 
of the left and right leads. The proper solution is 

F,{x) = Fde) + [FR(e) - F^e)] | . (184) 

The distribution function inside the wire interpolates between the two distribution linearly. At low 
temperatures it looks like a two-step function, where the energy separation between the steps is 
the applied voltage, eV, while the relative height depends on the position x. Comparing ( |183t with 
the continuity equation, one notes that the current density (at a given energy e) is given by j{€) = 
D dxF^x) - D[Fs(e) - Fi(e)]/L. The total electric current, is thus 

r evD r ,vD 

where the Drude conductivity of the diffusive wke is given by (Td - e~vD. 

6.5 Fluctuations 

Following the discussions in previous sections we consider fluctuations near the stationary solution 
Q ,(r) - A,-,', see (11651 ). We restrict ourselves to the soft-mode fluctuations that satisfy - I 
and neglect all massive modes that stay outside of this manifold. The massless fluctuations of the 
2-matrix may be parameterized as 

Q^i{oe-'^l^o&,oe''''^-oqA-\ (185) 

where rotation generators are given by 

' <u^ii-^4 \ ' \- (186) 



^ £/ / ' I -1 , 

Here dtfir) and dtfir) are two independent Hermitian matrices in the time space. One, thus, under- 
stands the functional integration over Qtf (r) in ( |173l l as an integration over two mutually independent 
Hermitian matrices in the time domain, dti'{r) and J„'(r). The physical meaning of t/„'(r) is a devia- 
tion of the fermionic distribution function F„'{r) from its stationary value. At the same time, d,t'(r) 
has no classical interpretation. To a large extent, it plays the role of the quantum counterpart of 
t/„'(r), that appears only as the internal line in the diagrams. The reason for choosing Q in the form 
of ( 11851 ) can be justified as follows. First, one notes that 2 = A = 'ZYiXj'i/ Second, one should 

realize that the part of 'W that commutes with Q does not generate any fluctuations, therefore, one 

restricts 'W to satisfy 'W o"; + &^ 'W = 0. Thus, 'W has to be off-diagonal and most generally 
parameterized by two independent fields, d and d, see ( 1186b . 

One may expand now the action ( 1172b in powers of J„'(r) and di,'{r). Since Q , was chosen to 
be a stationary point, the expansion starts from the second order If stationary F, ,>{r) is spatially 
uniform, one obtains 

iS ["Wj^-Y J dr JJ dtdt' d„. (r) [-D dl + d, + d,.\ d,,,{r) . ( 1 87) 
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The quadratic form may be diagonalized by transforming to the energy/momentum representation 
^^^,(q) = J dr JJ dfdf'^„'(r)exp(/ef- /e'f')exp(-iqr). 

As a result, the propagator of small Q matrix fluctuations is 

2 (5^ (5 2 

<flf,,,,(q)4^(-q))^ = Trf^ = ^.,.3^...4 ^(q> ^) > (188) 

TTV Dq-^ + loj nv 

where cj = ei - €2 ^ €3 - £4 and object lf^'^\q, co) = lf^'^\q, ei - 62) = [Dq'^ + /(ei - 62)]"' 
is called the dijfuson. The higher-order terms of the action ( 11721 ) expansion over £/,/'(r) and dn'iv) 
describe non-linear interactions of the diffusive modes with the vertices called Hikami boxes ll78l 
l79l . These non-linear terms are responsible for weak-localization corrections ll79l [8OI [STl l82l . If 
the distribution function F,i'(r) is spatially non-uniform, there is an additional term in the quadratic 
action -(7rvD/2)Tr{ii((9r/^)^(5r/^)}. This term generates non-zero correlations of the type {dd)'w, 
which are important for some applications. 



6.6 Applications II: Mesoscopic effects 
6.6.1 Kubo formula and linear response 



It was demonstrated in Section 15.41 how the linear response theory is formulated in the Keldysh 
technique. Let us see now how the polarization operator of the disordered electron gas may be 
obtained from NLSM action. To this end, one uses general definition of the density response function 
n''(x, X ') given by ( 1121b along with the disorder averaged action ( |174t , which gives 



n {x,x ) 



2 6V'''ix')6Viix) 



= y6(r - r')6{t - f) + Unvf {Tr{rQ„ir)]Tr{r'Q,','(r')}) 



(189) 

where x - (r, t) and angular brackets stand for the averaging over the action ( |172| i. The first term 
on the right-hand side of ( |189t originates from the differentiation of Trfyo- j V) part of the action 
( 11741 ). while the second term comes from diff'erentiation of TrfyQ}. Equation ( 11891 ) represents the 
cr-model equivalent of the Kubo formula for the linear density response. 
In the Fourier representation the last equation takes the form 

n^(q, o.) = y + Ujivf jj ^ (Tr{f''e,+„,,(q)}Tr{r'(2.',.'+„(-q))>^ . (190) 

Employing ( |185t and ( |186t . one finds in the liner order in the diffusive fluctuations (the only contri- 
bution in the zeroth order is v, indeed Tr{y'''A) = 0) 

Tr{r'e,,,,,+^(-q)) = d,,,,,U-q){F,'^ca - F,,) , 

Tr{f*e,+„,,(q)} - d.^^Am - PeF (q). (191) 

Since {dd}^ = only the last term of the last expression contributes to the average in (|190t . The 
result is 



1 + 



vDq^ 



n(q,w) = v+— — d£{F, - F,+^){d,+^,,{q)d,_,+U-q)) - , . . , 

4 \ [ Dq--ico\ Dq- - iio 

(192) 

where we have used the propagator of diffusons ( 11881 ) and the integral j de (Fc - F^+u) — -2a>. The 
fact that n^(0, w) = is a consequence of the particle number conservation. One has obtained the 
diffusion form of the density-density response function. Also note that this function is indeed re- 
tarded (analytic in the upper half-plane of complex ai), as it should be. The current-current response 
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function, K {q, co), may be obtained in the similar manner. However, more straightforward way is 
to use continuity equation q ■ j + <yp = 0, which implies the following relation between density and 
current response functions A'^(q, oj) = w^n*(q, cj)/q'. As a result the conductivity is given by 

cr(q, w) = — K'^iq, oS) ^ — n^(q, oS) = e'vD — , (193) 

iuj q~ Dq-^ - icj 

which in the uniform limit q — > reduces to the Drude result ctd = cr(0, u)) - e^vD. 



6.6.2 Spectral statistics 

Consider a piece of disordered metal of size L such that L » I, where / s UfTf.\ is the elastic mean 
free path. The spectrum of the Schrodinger equation consists of a discrete set of levels, e„, that 
may be characterized by the sample-specific density of states (DOS), v(e) = 2„ d{e - e„). This 
quantity fluctuates strongly and usually cannot (and does not need to) be calculated analytically. 
One may average it over realizations of disorder to obtain a mean DOS: (v(e))dis. The latter is a 
smooth function of energy on the scale of the Fermi energy and thus may be taken as a constant 
{v(e/r))dis = V. This is exactly the DOS that was used in the previous sections. 

One may wonder how to sense fluctuations of the sample-specific DOS v(e) and, in particular, 
how a given spectrum at one energy e is correlated with itself at another energy e'. To answer this 
question one may calculate the spectral correlation function 

/?(e,e') = <v(e)v(e'))dis-v'. (194) 

This function was calculated in the seminal paper of Altshuler and Shklovskii Il83l . Here we derive 
it using the Keldysh NLSM. 
The DOS is defined as 

y(e) = /^(G«(k, e) - G^(k, e))/{2n) = «(AnAi> - <-A2«A2»/(2;r) = -<#<T,^)/(27r) , 

k 

where the angular brackets denote quantum (as opposed to disorder) averaging and the indices are 
in Keldysh space. To generate the DOS at any given energy one adds a source term 

iSdos^- j ^JeJ dr%,r)cr,^(€,r) ^ - Jj drdf' J dr#(r, f)y,_,,(r,»#(r, f') , 

to the fermionic action (1172b . After averaging over disorder and changing to the Q matrix represen- 
tation the DOS source term is translated to 



iS DOS ^^v J J drTi{Q^^{r)d--} . 

Then the DOS is generated by v(e) - 6Z[J]/SJir. It is now clear that {v{e)}^i^ = jv{T^i{Q€eO-z)}q- 
Substituting Q^^ - one finds (v(e))dis = v, as it should be, of course. It is also easy to check that 
the fluctuations around A do not change the result (all the fluctuation diagrams cancel owing to the 
causality constraints). We are now in a position to calculate the correlation function ( 11941 ), 



R{£, £) = V -V 

6Je6Je' 



\ (Tr{4e<^^)Tr{4v'<^z))o - 1 



(195) 



Employing the parametrization of ( |185t , one finds, up to the second order in the difiiisive fluctuations 
Tr{ea-,} = i[4-2FoJ-2joF + c/oJ + Joc/]. (196) 
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Since {dd)'w - 0, the only non-vanishing terms contributing to ( 11951 ) are those with no d and d at 
all (they cancel term) and those of the type {dddd)^/. Collecting the latter terms one finds 

R{e, e') = ^ Jdr ^ {{d,,^., + d,,,d,,,) [d,,,M,, + d,,^_d,,,,))^ . (197) 

Now one has to perform Wick's contractions, using correlation function {dee'de'e)'W °^ I^{e — e'), 
which follows from ( 11881 ). and also take into account J dei [^^'^'(q, e - ei)]^ — 0, owing to the 
integration of a function which is analytic in the entire upper (lower) half-plane of ei . As a result. 



(198) 



where the momentum summation stands for a summation over the discrete modes of the diffusion 
operator Dd^ with the zero current (zero derivative) at the boundary of the metal. This is the result 
of Altshuler and Shklovskii 1831 for the unitary symmetry class. Note that the correlation function 
R{e, e') depends only on the energy difference o) - e - e' . Diagrammatic representation of R(e, e') 
function is shown in Figure [TT] Adopting an explicit form of the diffusion propagator, we find 
spectral correlation function in the form 



7?(e-e') = -:7^Rey 



1 



(199) 



where ql - Y^ii^^n^^il with ^ - x,y,z', - 0,1,2... and are spatial dimensions of the 
mesoscopic sample. 



(q,e-e') 




Figure 11: Diagram for calculation of mesoscopic fluctuations of the density of states, 
R(e,e'), see ( 1195b . It is generated from the Wick contraction {d^f^deiEd^'^^_d^^e'}'w — » 
{decide' e2)'w{deiede2e')'\v <^ [^''^(q, f - e')]^Sete' 6e2e, scc ( 1197b . There is also a similar diagram with 
the advanced diffusons. 

For a small energy difference co Ej-i, = D/L^ only the lowest homogenous mode, q„ = 0, 
of the diffusion operator (the so-called zero mode) may be retained and, thus, R{oj) - -I l{2jr<J^). 
This is the universal random matrix result. The fact that the correlation function is negative means 
that the energy levels are less likely to be found at a small distance a) from each other This is a 
manifestation of the energy levels repulsion. Note that the correlations decay very slowly — as the 
inverse square of the energy distance. One may note that the random matrix result 1841 

1 - co^ilnojIS) 

RrmAcS) = r-r^ , (200) 
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where 6 is the mean level spacing, contains also an oscillatory function of the energy difference. 
These oscillations reflect discreteness of the underlying energy spectrum. They cannot be found in 
the perturbation theory in small fluctuations near the A "point". However, they may be recovered 
once additional stationary points (not possessing the causality structure) are taken into account ll85l . 
The saddle point method and perturbation theory work as long as w » 5. Currently it is not known 
how to treat the Keldysh NLSM at w < 5. 



6.6.3 Universal conductance fluctuations 

Similarly to the discussions of the previous section consider an ensemble of small metallic samples 
with the size L comparable to the electron phase coherence length, L ~ L^. Their conductances 
exhibit sample-to-sample fluctuations owing to diff'erences in their specific realizations of disorder 
potential. These reproducible fluctuations are known as universal conductance fluctuations (UCFs). 
Theoretical studies of UCFs were initiated by Altshuler ll86l . and Lee and Stone ISTl . More detailed 
study of UCFs was given later in ll83l[88l . The technique developed in this section for treating UCFs 
in the framework of the Keldysh non-linear cr-model is closely parallel to the theory developed by 
Altshuler, Kravtsov and Lerner [[89l in the framework of the zero temperature replica non-linear 
sigma model. 

Our starting point is the expression for the dc conductivity within the linear response given by 



„2 



1 



2 n^o Q 



(201) 



,M;;'(n)A^(-n)j^,,^^,^„ 

where indices jj., v stand for the spatial Cartesian coordinates. Expanding action ( 11741 ) to the quadratic 
order in the vector potential with the help of ( 1175b one finds that corresponding term in the generat- 
ing function reads as J3[A''', A*] = ^{Tr{AgAg))g. At the saddle point 2 = A, after consecutive 
diff'erentiation over the vector potential in (12011 ) one finds for the average conductivity 



<cr^v>dis = Sfiv lim ^^Tr iy''' A,+nV' K-n] = Sfiv^^^ lim ^ f ^{F.+n - F.^n) = croSf^y (202) 

where cto = e^vD, as it should be of course. At this level, retaining fluctuations 'W of the Q matrix 
around the saddle point A, one can calculate weak-localization corrections 17811791 [80l[8Tll82ll to the 
average conductivity. In what follows we are interested in calculation of the irreducible correlation 
function for the conductivity fluctuations which is defined in the following way 

In view of ilQli this irreducible correlator can be expressed through the Q matrix as 



<^^-^'^-^>- - i^J n i 6A^,ia^4,i-J (T'-iAeAeiTriAeAei) 



cTn^Mivi^M^v, , (203) 



where we have used ( 1174b and expanded exp(;5' [2, A]) up to the forth order in the vector potential. 
Now one has to account for fluctuations of the Q matrix up to the second order in generators 'W. 
There are two possibilities here: within each trace on the right-hand side of (1203b one may expand 
each Q matrix either to the linear order in "W resulting in 7~i ['W] - Ti{A&~'WA&-'W}; or alterna- 
tively set one of Q matrices to be A, while expanding the other one to the second order, resulting in 
T2['W] = TT{Ad-.Ad-,'W^], where A = tl'^Atl. As a result, ( i203b takes the form 

6^ 

6A'Jmi)dAU-Q.d 



<5cr^,v,(5cr^,v2>dis = (^) f] him ^ 



[{Ti V\V]Ti VW])^ + {T2VW]T2VW])^] - o-lS,,y,6,,y, . (204) 
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a) ^ b) 

Figure 12: Diagrams for the variance of conductance fluctuations. 

It is convenient to represent each average here diagrammatically, see Figure [12] A rhombus in 
Figure [T2h correspond to the term with T\Y\V], where the opposite vertices represent matrices A, 
while rectangles with adjacent vertices in Figure [T2b correspond to the term with T2['W]. The 
vertices are connected by the diffuson propagators of the field 'W. Equation ( I204l i should also 
contain the cross-contribution 2{Ti['W]T2['W]}<w, which vanishes, however, upon "W averaging. 

Differentiating each term of the ( |204t individually, multiplying matrices and using diffuson prop- 
agators from (II88I 1. one finds for (1204b 



deide2 



ny]JL^ .^.12 



where 



[2T cosh(ei /2T) cosh(e2/2r)] 

)+Rs[lfiq, 61-62)]% . (205) 

9 

The first term in the square brackets of ( 12051 ) corresponds to Figure [T2h and the second one to 
Figure fT2b. Introducing ei - 62 = w and ei +62 = 2e, and integrating over e. Equation (12051 ) may be 
cast into the form 

{5cr^,v,5cr^,v,)dis - 0-j(5^,^,5v,v, + ^/Jiv,5v,;j,) + CTj 5^i,v,5/j2V2> (206) 

M' r-r(-)y^^, (207a) 

and dimensionless function is given by Tix) - [xcot\\{x) - 1]/ sinh~(x). Here a\ may be regarded 
as contribution from the mesoscopic fluctuations of the diffusion coefficient. Figure fT2h . while (x\ 
as the corresponding contribution from the fluctuations of the density of states. Figure [T2b . The 
fact that (7~i ['VV']7~2['VV])w - implies that mesoscopic fluctuations of the diffusion coefficient and 
density of states are statistically independent. In general, <j\ and cr^ contributions are distinct. At 
zero temperature w — » 0, however, they are equal, resulting in 

<5cr^,y;5cr^,vj) - Cd ^^^j ('^AiiAi2^viV2 + ^//iV2^vi;/2 + ^lnviS^^y^) , (208) 

where - {4-/7t)~ Ym^i^^ti^ti) '^ dimensionality- and geometry-dependent coefficient (note that 
in the final answer we have restored Planck's constant). This expression reflects the universality 
of conductance fluctuations and, of course, coincides with the result obtained originally from the 
impurity diagram technique Il83l[87l . for review see 
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6.6.4 Full counting statistics 

When current /(f) flows in a conductor it generally fluctuates around its average value {/). One is 
often interested in calculation of the second, or even higher moments of current fluctuations. An 
example of this sort was already considered in the Section l5.5.2l Remarkably, in certain cases one 
may calculate not only a given moment of the fluctuating current, but rather restore full distribu- 
tion function of current fluctuations. Theoretical approach, utilizing Keldysh technique, to the full 
counting statistics (FCS) of electron transport was pioneered by Levitov and Lesovik and cowork- 
ers II9TI |92l |93l . In the following we consider its application to the diffusive electronic transport 
developed by Nazarov 19411 . 

Consider two reservoirs, with the chemical potentials shifted by externally applied voltage Y . 
It is assumed that reservoirs are connected to each other by difiiisive quasi-one-dimensional wire 
of length L. The wire conductance is g^ = o-qJIIL, with J[. being wire cross section. Describing 
difiiisive electron transport across the wire one starts from the disorder averaged generating function 
Z\x\ = / D[2] exp(/5' [2, A^]). The action is given by ( I174l i. while the auxiliary vector potential 

enters the problem through the covariant derivative ( 1175b . We choose to be purely quantum, 
without classical component, as 



f^i X 0<t<tQ 

2L 1 otherwise 



Here the quantum Keldysh matrix is given by ( |113t andx is called counting field. The action 
S [2, A^] is accompanied by the boundary conditions on Q{x) matrix at the ends of the wire; 

2(0) - ( I ) , Q(L) = ( I j . (210) 

Knowing Z,\x] one can find then any moment {q") of the number of electrons transferred between 
reservoirs during the time of measurement to via differentiation of Z,[x\ with respect to the counting 
field;if. The irreducible correlators are defined as Ci - {q) - qo andC„ = {{q-qo)") with « = 2, 3, . . ., 
where q - ^ J^^" I(t)dt and ^0 - fogoW^ - h){I}/e, where go is the average diffusive conductance. 
They may be found through the expansion of the logarithm of Z,lx] in powers of the counting field 

lnZM = y^C„. (211) 

One calculates Xlx]^ by taking the action at the saddle point Q-Kj^ which extremizes S [Q, A^]. 
The difficulty is that the action S [Q, A^] depends explicitly on the counting field and solution of the 
corresponding saddle point equation is not know for an arbitrary A^. This obstacle can be overcame 
by realizing that vector potential (12091 ) is a pure gauge and it can be gauged away from the action 
S [2' A;^] — » S [Qx] by the transformation 

Q{x ; f, f') - exp [ixA^it)] Q^ix ; f, f') exp { - ixA^it')] . (212) 

It comes with the price though: the boundary conditions ( 121 Oi l change accordingly 

G^(0) = e(0), e^(L) = exp ( - ixr/2)Q(L) exp {ixf /2) . (213) 

The advantage of this transformation is that the saddle point equation for Qj^, which is nothing else 
but the Usadel equation ( 11811 ) 



4 

ox 



dx 



Qx° — 



0, (214) 



can be solved explicitly now. To this end, notice that ° dxQx - ~dxQx ° Qx = is a constant, 
i.e. x-independent, matrix. Since - I, J anti-commutes with Qx, i.e. o / + y o = 0. As a 
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result one finds = Qj^{Q)exp{Jx). Putting x - L and multiplying by Q^{0) from the left, one 

expresses as-yet unknown matrix J through the boundary conditions (I213l i: J = j^ln [Q^{Q)Q^{L)]. 

Having determined the saddle point configuration of the Qj^ matrix, for an arbitrary choice of the 
counting field one substitutes it back into the action S [Qj^] to find the generating function 

InZM = iS[Q;,] = TriOMxf} = ^ TrlT^) , 

where we used anti-commutativity relation (2^(0), J} = 0. Calculating time integrals one passes to 
the Wigner transform ^dfdf' — > fo / , where fo emerges from the integral over the central time, 
and finds 

InZM = Iff / deTTln^[Q(0)sxp{-ixr/2)Q(L)exp{ixr/2)] • (215) 

In the following we analyze ( |215t in the zero-temperature limit, T = 0, where - tanh(e/2r) 
sign(e). The algebra can be further significantly shortened by performing the rotation Q = O^QO 
with the help of the matrix 

1/1-1 



One should note also that expi+ixj'^ /2)0 - exp(+/;\^(T /2). It is not difficult to show that for 
r = the only energy interval that contributes to the trace in ( I215l l is that where < e < eV. 
Furthermore, at such energies rotated (3-matrices are energy independent and given by 

QiO)^[-l Q(L)^[_l (217) 

As a result, the e integration in Eq. ( 1215b gives a factor eV and inserting Q into InZDif] the latter 
reduces to 

InZM^-^Trln^ _^ j . (218) 

Since the trace is invariant with respect to the choice of the basis, it is convenient to evaluate it in 
the basis where matrix under the logarithm in ( 12181 ) is diagonal. Solving the eigenvalue problem and 
calculating the trace, as the final result one finds 



InZM = — — In- 

4e 



.Jpf^], Pj,^2e'^-l. (219) 



Knowing In ZDf ] one can extract now all the cummulants of interest by expanding in powers of x 
and employing (121 lb . for example, C\ - qo, C2 = ^o/3, C3 - qo/lS, etc. For a review devoted to 
FCS see ||95J. 



7 Interactions and kinetic equation for fermions 
7.1 Interactions 

Consider a liquid of electrons that interact through the instantaneous density-density interactions 

i/mt = ^ drdr' : e(r)f/o(r - r')e(r') : , 

where ^(r) = {fi'\r)^{r) is the local density operator, Uoir - r') is the bare Coulomb interaction 
potential and : . . . : stands for normal ordering. The corresponding Keldysh contour action has the 
form 

S ,nt [lA, -A] = - ^ J dr jj drdr' 0-(r , ^r' , t) Uo(r - r')if,{r' , t)tfr(r, t) . (220) 
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One may now perform the Hubbard-Stratonovich transformation with the help of a real boson field 
(f>{r, t), defined along the contour, to decouple the interaction term 



exp (i5i„t[iA, "A]) = J D[0] exp ^ df ^ drdr>(r, f) C/„ i(r - r') 0(r', f)j 

X exp {a ^ df ^ dr 0(r, f)0-(r, f)i/r(r, f)j , (221) 

where f/^' is an inverse interaction kernel, i.e. Jdr"f/o(r - r")UQ^{r" - r') = 6{r - r'). One 
may notice that the auxiliary bosonic field, 0(r, f), enters the fermionic action in exactly the same 
manner as a scalar source field (1112b . Following dl 14b . one introduces (p^'"-''"' = (0+ ± (f>-)l2 and 
rewrites the fermion-boson interaction term as \iia(lfY^iy4'h , where summations over a,b = (1,2) and 
a = (cl, q) are assumed and gamma matrices y are defined by dl 13b . The free bosonic term takes 
the form ^(pU^^cj) — » <p"lJ^^lTf ^ . Following ( 1221b one may integrate fermions explicitly to obtain 
the partition function for the interacting disordered electron liquid 

= J D[0] exp(/Tr{<5^t/o'(T,-(5)) J D[e] exp(i5[e, <D]) , 



z 

{222) 



4Tel 



G"' + — g + d + VfA 

2Tel 



where we introduced doublet <P^ - {(f)'^', (p'^) and matrix <!) = (f>"y" . This should be compared to 
the non-interacting version of the action given by (1174b . An extra complication, which stems from 
interactions, is an additional functional integral over the dynamic bosonic field O entering (1222b . 

Varying the action in ( 1222b over the Q matrix 6S [ 2, <1>] / 6Q - 0, at zero external vector potential 
A = 0, one obtains the following equation for the saddle point matrix Q - Q[^\. 



,(r) = — G-' + -^2 + 6) , (223) 



which is a generalization of (|162b for the interacting case. Our strategy will be to find a stationary 
solution of ( 12231 ) for a given realization of the fluctuating bosonic field 6, and then consider space- 
and time-dependent deviations from such a solution. 

The conceptual problem here is that the saddle point Equation ( 12231 ) can not be solved exactly for 
an arbitrary <l)(r, t). Note, however, that ( 1223b can be solved for a particular case of spatially uniform 
realization of the boson field, Q> - 0(f). This is achieved with the help of the gauge transformation 
of the non-interacting saddle point 

g^, mm = exp J/ J df 4)(f) j A,_,, exp (-/ J df 6(f) j . (224) 

The validity of this solution can be verified by acting with the operator G"' + il{2T^\)Q + 6 on 

both sides of (1223b . and utilizing the fact that A,_,' solves ( |223b with 6 = 0. We also rely on the 
commutativity of the vertex matrices [y^', y*] = 0, in writing the solution in the form of ( 1224b . 
This example shows that a properly chosen gauge may considerably simplify the task of finding the 
saddle point and performing perturbative expansion around it. We show in the following that there 
is a particularly convenient gauge (the % gauge) suited for calculations of interaction effects. 

7.2 % Gauge 

Let us perform a gauge transformation from the old Q matrix to a new one, which we call Qjc matrix. 
It is defined as 

e3c(r; f, f') = exp [-i%{r, t)) Q„,(r) exp (/D<:(r, f')) , (225) 
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where the matrix 3C(r, t) - OCir, f)-y" is defined through two scalar fields OC^ir, t) with a - (cl, q), 
which are specified below. Substituting Q - e'^Qxe^'^ into the action (1222b and using the invari- 
ance of the trace under a cyclic permutations, we can rewrite the action as[ 



;5[e3c,4'] = --piTr{e^) + Trln 

4Tel 



/Tel 2ot 



(226) 



where we have introduced the notation C{r, t) - OjcCr, t) + v^- Ag<;(r, t) along with the gauge trans- 
formed electromagnetic potentials 

a)3c(r,0 = <i>(r,f)-5,a(:(r,f), A3c(r,f) = A(r,f)- Wr,f). (227) 

We assume now that the saddle point of the new field Q'x is close to the non-interacting sad- 
dle point A, ( 1165b . and use the freedom of choosing two fields %" to enforce it. To this end, we 
substitute Q-x = A + dQ-x into ( |226b and expand it in powers of the deviation SQ-x as well as the 
electromagnetic potentials, encapsulated in C. The first non-trivial term of such an expansion is 

iS [6Qx, <I)] = -^Tr{gcgSQx] + ■■■, (228) 

where we have employed the fact that A is the saddle point of the non-interacting model and, thus, 
in the absence of the electromagnetic potentials, there are no linear terms in deviations 6Qx- We 
have also neglected the diamagnetic (drOQ^/lm term, since it is quadratic in %, and hence (as shown 
below) in <1). 

We now demand that this linear in SQx,tt'(i^) term vanishes. Performing the Fourier transform, 
one notices that this takes place for an arbitrary 523c,e_e+(q), if the following matrix identity holds 
for any e, u and q 

^(P+, f+)C(q, ^m9-, f-) = , (229) 

p 

where p-t = p ± q/2 and €± - e + ull. Condition ( |229b represents matrix equation, which ex- 
presses as-yet unspecified gauge fields %" through <1>" and A". Employing ( |169b . and the following 
identities 

Yj ^''(P=t' e±)^^(p^ , e^) * 27TVT,x , (230a) 
p 

^*(P±' f±)^^(P^> f^) ^ +27rivreiZ)q , (230b) 

one may transform (|229b into 

y ^(p+,e+)C(q,w)^(P-,e-) = {T - KT ^^^-Wx ' iKT - T K_)D diw =0. (231) 

TTVTel ^ 

It is in general impossible to satisfy this condition for any e and w by a choice of two fields %"{r, ai). 
In thermal equihbrium, however, there is a "magic" fact that 

\" ^coth— s jj^, (232) 

Fe^ - Fe_ 2T 



'■^Deriving {226\ one uses obvious equality between the traces Tr{Q^| = Tr{Q-|. As to tlie logarithm term, one writes 
Tr je"'* In [G"' + i + Vf A + j^^e'^ Qxe"'*] <?'* ) = Trln [e"'*G"'e'* + <! + Vf A + Qx], where familiar algebraic 
identity Tr{Lf{A)L] = Tr|/(LAL-')| was used, which holds for any analytic function / of matrix A. Finally, one rewrites 
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which depends on u only, but not on e. This allows for the condition ( 12311 ) to be satisfied if the 
following vector relation between the gauge transformed potentials ( 1227b holds; 

<53c(r, ^) = ( J ^^"^ j D divXocCr, . (233) 

This equation specifies the 3C-gauge for both classical and quantum components of the electromag- 
netic potentials. 

The advantage of the 3C gauge is that the action does not contain terms linear in the deviations 
of the Q'x matrix from its saddle point A and linear in the electromagnetic potentials. Note that 
there are still terms which are linear in bQ'x and quadratic in electromagnetic potentials. This means 
that, strictly speaking, A is not the exact saddle point on the Q-x manifold for any realization of the 
electromagnetic potentials. However, since the deviations from the true saddle point are pushed to 
the second order in potentials, the % gauge substantially simplifies the structure of the perturbation 
theory. Moreover, this state of affairs holds only in equilibrium. For out-of-equilibrium situations 
condition ( 12311 ) cannot be identically satisfied and terms linear in bQ'x and electromagnetic fields 
appear in the action. As we explain below, it is precisely these terms which are responsible for 
the collision integral in the kinetic equation. Still the 3C gauge is a useful concept in the out-of- 
equilibrium context as well. In such a case one should define the bosonic distribution function 
in (12331) as 



1 r+°° 

B^{X,T) = — de [1 - Ff+„/2(r, T)Ff_„/2(r, t)] 

^^^> J-oo 



(234) 



where Ff(r, r) is WT of the fermionic matrix F, ,'(r). 

With the help of (12271 ) the definition of the 3C gauge ( I233I I may be viewed as an explicit relation 
determining the gauge fields 3C" through the electromagnetic potentials and A". Taking A = 
for simplicity, one finds for the quantum and classical components of the gauge field 

(Ddl - iaj)3CHr, co) = (D''(r, w) , (235a) 

(M? + iaj)X'='(r, oj) + 2B^Ddl%Hr, w) = -0'^'(r, w) . (235b) 
In general case it is convenient to cast these relations into the matrix form 

aC(q, cj) = £)-i(q, w)(S„''5(q, cS}-D&,q- X(q, w)) , (236) 
with the vector OC^ = %''). Here we have introduced diffiison bosonic matrix propagator 

having matrix components 

2/<'*'(q, OJ) = {Dq^ + ia)y\ £)^(q, oj) = B„[2/(q, oj) - £)^(q, w)] , (238) 

and 



2B,., 1 



R 



Equation ( 12361 ) provides an explicit linear relation between the gauge fields %" and the electro- 
magnetic potentials. It thus gives an explicit definition of the gauge transformed field Q-x, cf. ( |226l l. 
The latter has the saddle point which is rather close to the non-interacting saddle point A (with de- 
viations being quadratic in electromagnetic fields). Returning to the original gauge, one realizes that 
the following Q matrix 

(r) = exp {i%"{r, f)f " ) A,_,, exp (-;3C^(r, t')f) , (240) 
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provides a good approximation for the solution of the generic saddle point Equation (12231 ) for any 
given realization of the fluctuating potentials. This statement holds only for the equilibrium condi- 
tions. Away from equilibrium, ^6Qx terms reappear and have to be taken into the account to obtain 
the proper form of the kinetic equation (see further discussions in Section ITSl l. In addition, terms 
~ ^^6Qjc exist even in equilibrium. They lead to interaction corrections to the transport coefficients 
(details are given in Section lT!6] l. 



7.3 Non-linear cr-model for interacting systems 

Performing gradient expansion for the trace of the logarithm term in ( I226l l (this procedure is closely 
analogous to that presented in Section |672] i. one obtains an effective action written in terms of Qx 
matrix field and electromagnetic potentials in the % gauge 

iS[Qx,<b] = yTr{<l)3c<T,4)3c} - ^Tv[D(drQoc)- - 4d,Qoc + 4i^xQ3c} , (241) 

where 

drQx = drQx - ilAjc, Qx] ■ (242) 

Equation ( 12411 ). together with the saddle point condition ( |236l )-( |239l ), generalizes the effective cr- 
model action ( |172t to include Coulomb interaction effects. Employing the explicit form of the long 
covariant derivative ( I242l i. and the relation between the 3C and <1) fields at A = (see ( |235t ), one 
finds for the partition function 

Z = J D[0]exp(/Tr{0^;7«;,^(5)) J mQx] {'S o[Qx] + iS dQx, drT] + iS.lQx, drOC]) , 

(243) 

where Si, with / = 0, 1, 2 contain the l-th power of the electromagnetic potentials and are given by 
iSolQx] = ~TT{D(drQxf - MdtQx] , (244a) 

iSdQx, drOC] = -i7TvTT[D{dS)QxidrQx) + ^xQx} , (244b) 
iS2[Qx,drX] = ^TT{(dr±)Qx(dS)Qx - (dr±)A(dS)A] . (244c) 



The effective interaction matrix Urpa is nothing but the screened interaction in the random-phase 
approximation (RPA) 

URPAiq, cj) = [C/o '<^.v + n(q, w)]"' , (245) 

where n(q, oS) is the density-density correlator. According to (I121l i and (|192| i it has a typical form 
of a bosonic propagator in the Keldysh space 



"('I'^^^^l n«(q,^) n-(q,J) I' (246) 



with the components 



n«(^>(q, oj) = ^''f!" , n^(q, CO) = B^[n«(q, co) - n\q, co)] . (247) 
Dq~ + 10) 

To derive ( I243l l-( I247| | one has to add and subtract the term Tr{(5r^K!)A(5r^K!)A}, and employ the 
equation 

deTr [rf - rKfK] = ^i^iK'f > (248) 
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where e± - e+oj/l, and matrices A and S are defined by ( 1165b and (|239t correspondingly. Equation 
(12481 ) is a consequence of the following integral relations between bosonic and fermionic distribution 
functions 

X + OO p + CX) 

d£(F,^-F,_) = 2aj, Ae{\-F,^F,_)^2uB,,. (249) 

Equations ( I243b - (l247b constitute an effective non-linear cr-model for interacting disordered 
Fermi liquid. The model consists of two interacting fields: the matrix field Qx, obeying non-linear 
constraint - 1, and the bosonic longitudinal field drOC (or equivalently <1)). As will be apparent 
later, Qjc field describes fluctuations of the quasi-particle distribution function, whereas (or 3C) 
represents propagation of electromagnetic modes through the media. 

7.4 Interaction propagators 

For future applications we introduce correlation function 

^"/^{r - r', f - f') = -2/<3C''(r, f)ac'^(r', f')> = -2iJ D[<1)] X^ir, f)3C''(r', f') exp (!Tr{<5^t/^^^<5)) , 

(250) 

where factor ~2i is used for convenience. Since O and % are linearly related through (12361 ), one may 
evaluate this Gaussian integral and find for the gauge field correlation function 

^^(q,^) = D(q,cj)^-JURPA(q,oj)CB-_'jD^(-q,-aj). (251) 

The bosonic correlation matrix ^{q, to) has the standard Keldysh structure 

with the elements 

n/«(->(q, a>) = - \. ' + , (253a) 

'V'^iq, CO) = B^l-V^iq, to) - ^\q, to)] . (253b) 

This propagator corresponds to the screened dynamic Coulomb interaction, dressed by the two dif- 
fusons at the vertices. Figure [T3b . Thus, the role of the gauge field % is to take into account au- 
tomatically both the RPA-screened interactions. Figure [T3b . and its vertex renormalization by the 
diffusons. Owing to the linear dependence between <1) and 3C, (see ( 1236b ). we use averaging over <1) or 
DC fields interchangeably. The essence is that the correlator of two "X" fields is given by ( I250l )-( l253b . 



7.5 Kinetic equation 

The aim of this section is to show how the kinetic equation for the distribution function F appears 
naturally in the framework of the Keldysh formulation. In Section|63]it was demonstrated that the ki- 
netic equation for non-interacting fermions is nothing but the saddle point equation for the effective 
action of the Q matrix. In the case of interacting electrons it is obtained from the action S [Qoc, <1>], 
(see ( 1241b ). by first integrating out fast degrees of freedom: diffusive, 'W, and electromagnetic, % 
(or, equivalently, O). 

Let us outline the logic of the entire procedure, which leads from the partition function ( |243b 
and ( 1244b to the kinetic equation. As the first step we separate slow and fast degrees of freedom in 
the action S i[Qx,dr'K], where / = 0, 1,2 (see ( 1244b ). The former are encoded in the distribution 
function Fi,'{r), while the latter are carried by diffusons "Wn'ir) and electromagnetic modes !}C(r, t). 
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Figure 13: a) Diagrammatic representation of the gauge field propagator 'V{q, co): wavy line rep- 
resents Coulomb interaction. Vertices dressed by the diffusons are shown by the ladders of dashed 
lines, b) Screened Coulomb interaction in RPA, URPAi(l,CL>)- Bold and thin wavy lines represent 
screened and bare interactions correspondingly, the loop represents polarization operator dressed by 
the diffusion ladder. 

This separation is achieved by an appropriate parametrization of the Qac -matrix. One convenient 
choice is Qx = 'Z/- ° Qfast ° where rotation matrices 

(254) 



Z -1 / ' ^ \Z -\+ZoF ^ 

with Ao B - J dt'A„>B,>,» carry information about slow degrees of freedom, and the fast part of Qjc 
matrix is parameterized by the diffuson fields Qf^st - exp{-^/2) o &~ o exp{'iV/2) (compare this 
parametrization with that given by ( 1185b ). In the last equation Z„'(r) (not to be confused with the 
partition function) may be thought of as the quantum component of the distribution function F„'{r). 
Although Z„'(r) is put to zero in the end of the calculations, it was emphasized in the 1961 that Z„'(r) 
must be kept explicitly in Q parametrization to obtain the proper form of the collision integral in the 
kinetic equation. 

As the second step, one performs integrations over <1) (or equivalently 3C, since the relation 
between them is fixed by ( I236l l). and over 'W fields in the partition function ( I243l l, to arrive at the 
effective action 

Z = jD[eK,a>] exp(/5[^,5r3C]) = J T>[F,Z] expiiS ,dF,Z]) . (255) 

Note that after the decomposition given by (I185I I. with the 'Z/. and 'Z/r' matrices in the form of ( I254l i. 
one understands the functional integral over Qjc matrix in the (I255l l as taken over the independent 
matrix fields F, Z and 'W. As a result, the effective action S^ff will depend on F and its quantum 
component Z, and possibly the classical external fields, such as, e.g., scalar or vector potentials. One 
then looks for the saddle point equation for the distribution function F: 

SS,s[F,Z] 



6Z 



= 0, (256) 

z=o 



which is a desired kinetic equation. 

Proceeding along these lines, one expands the action (1244b in terms of F, Z, 'W, and electromag- 
netic potentials C) and %. For the slow part of the action one finds from (1244k ') that Tr{(5rQ3c)^} = 
8tr{c)ri^f/'5rZ,',) + 0(Z~) and Tr{dtQx} = 2tr{c)fZ„'F,', - dtFu'Zft], where tr{. . .) denotes the spatial 
and time integrations only and Keldysh structure was traced out explicitly. Passing to the Wigner 
transform representation ( 1182b . one obtains 

iSQ[F, Z] = 27rvtr {[M?F,(r, r) - d.F^r, T)]Z^(r, r)} , (257) 
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where t — (t + f)/2. Already at this stage, differentiating So[F,Z] with respect to Z one recovers 
from ( 1256b the non-interacting kinetic Equation ( 1183b . In a similar fashion, one finds dynamic part 
of the action for the fast degrees of freedom, 

iSol^] = -y tr[d,{r,T)[Ddl - dr]d,(r,T)] , (258) 

which is nothing else but Wigner representation of ( 1187b . 

We continue now with the coupling terms between the 'W and <i> modes. For S i ['W, F, Z] part 
of the action, which follows from (1244b ) upon expansion, one obtains 

iS 1 [OV, F, Z] = -invti [{[F, X'J] + Xl - FXlF)d + ([Z, Xl'] -Xl+ ZFXl + XlFZ)d] , (259) 

where 

XI = (D« - d,X" ± DdlX". (260) 

Deriving the functional relation between Q> and % fields, our logic was to nullify S \ part of the 
action (recall ( 12281 )). This step turns out to be impossible to implement for the non-equilibrium 
situation. However, we may still satisfy ( |235b ) by imposing a condition - 0. Although the 
Keldysh component of ( |231b cannot be satisfied identically, it still makes sense to demand that JC^' 
obeys the following non-equilibrium generalization of equation ( 1235b ) 

{Ddl + ;w)3C"'(r, w) + 2B„(r, r)M^aC«(r, w) = -cE)'''(r, w) , (261) 

where non-equilibrium bosonic distribution function is defined by (1234b . Note, however, that this 
generalization does not imply that linear in 'W (i.e. in d and d) terms vanish in ( 1259b . Indeed, us- 
ing (I235h ) which relates quantum components of O and %, and ( 12611 ), performing Wigner transform, 
one finds that S i ['W, F, Z] part of the action can be brought to the form 

iS 1 [OV, F, Z] = -OTvtr [l[F]Xl{r, aj)d,_(r, T)e"'"^ + Z,(r, r)X!'(r, Lo)[d,_(r, r) - d,^(r, T)]e"'"^} , 

(262) 

where e± = e + w/2 and we have introduced functional 

I[F] = B^{r, T)[F,_„(r, r) - F,{r, r)] + 1 - F.^Jr, T)F,(r, r) . (263) 

Note that, in equilibrium, I[F] = 0. In ( 1262b one keeps an explicit co dependence, thus not perform- 
ing expansion for small w as compared to e in the conventional Wigner transform sense. In addition, 
equation ( 1262b should also contain terms proportional to FZX'Ld, which will not contribute to the 
kinetic equation after % averaging, thus omitted for brevity. 

The remaining S 2 part of the action ( 1244b ) is already quadratic in the fast degrees of freedom 
Si (5r3C)-, therefore it can be taken at 'W = 0: 

iSilF Z] = 4nvD tr{(5r3C'')(5r3C«)Z - {drX'')FidrX'')FZ - {drX'')(dr%'')Z + {dr%'')F{dr'2C^)FZ 
. i6.X^')Zi8.X^')F - ii3.X^)id.X^')FZ - \i3.X^)id.X^')ZF} . (264) 

The next step is to perform the Gaussian integration over the fast degrees of freedom: diff'usons 
{d, d) and gauge fields (!}C^', %''). For S \ part of the action, employing ( 1258b and (12621) we obtain 

(exp {iSoVW] + iS 1 V\V, F, Z]))^^ = exp {iS'^^[F, Z]) , (265) 

where 

Z] = -4;7rvtr [{Dq^f[!D^{(i, oS)'V\({, u) - £)*(q, w)'V'^(q, w)]J[f ]z} . (266) 
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To derive S'"^^ in the form of ( 12661 ), one observes that upon integration the terms trlI[F]Xld} and 
t[{ZX'_'d} in ( |262] | prod uce an effective interaction vertex between F and Z, namely: { exp(/5' i))^ = 
exp {tT{I[F]Xl3y^ZX'^}). The latter has to be averaged over X, which is done observing that 

{x-^iq, co)Xl{-q, -co))^ = -4D^ {dlX^'iq, co)dl3CH-q, = -liiDg'-f^^iq, to) . (267) 

The last equation is a direct consequence of ( I260l l and ( 1261b . and correlator given by ( 12501 ). 
For S 2 part of the action, using ( |264t . one finds 

<exp(/52[F,Z])>3, = exp(/5f„'[F,Z]), (268) 

where 

iS f^[F, Z] = 2i7Tv tr [Oq^'V'^iq, w) - ^■^{q, tj)]I[F]z] . (269) 

To derive equation (|269l l one has to use interaction propagators for the gauge fields (12511 ). and adopt 
quasi-equilibrium FDT relation for the Keldysh component at coinciding arguments 

'V'^ir, r, t) = B^ir, r) [^''(q, oj) - ^^(q, oj)] , (270) 

which holds in the non-equilibrium conditions as long as F^{r, r) changes slowly on the spatial scale 
Lj = ^JD/T (this implies that gradient of F^{r, t) is small). The correction to the (1270b is of the 
form ocu j dr'lf(r - r', oj)drB^(r', T)d^D^(r' - r), see |22|. 

As the final step, one combines So[F,Z] from Eq. (1257b . together with S^^g^^\F,Z] parts of the 
action given by Eqs. (1266b and ( 1269b , and employs Eq. ( 12561 ) to arrive at the kinetic equation 

DdlF,(r, r) - 5,F,(r, t) = I,oi[F] , (271) 

where the collision integral is given by 

with the kernel 

M(q,w) = -iDq^[['y'^{q,(^)-^'^{q,aj)]-2Dq^[D^{q,Lj)'V'^{q,Lj)-D'^{q,(^)'V'^{q,u)^^ (273) 

This equation can be simplified by noticing that the gauge field propagator 'y*'^'^'(q, <y) may be writ- 
ten in terms of the diffiisons and screened RPA interactions, as 'V*(q, at) = -[£)*(q, w)]^C/^p^(q, a)) 
and similarly for the advanced component, which is direct consequence of (1245 b . Using this form of 
<y«(A)(^q^ w), after some algebra the interaction kernel Aiiq, at) reduces to 

M{q, w) = 2 Re[2>^(q, w)] Im[f/^p^(q, w)] . (274) 

For the conventional choice of the fermion distribution function n^{r, t) - {I - F^(r, r))/2, one 
can rewrite the collision integral ( 1272b in the usual form with "out" and "in" relaxation terms. Indeed, 
employing ( |249b , one identically rewrites the right hand side of ( 127 lb as ||97l |98l 

dwde' K(q, oj) [n^n^.-M - n^')( 1 - n^-^) - n^'ii^-U 1 - n^)( 1 - n^'-^)] , (275) 

oo 

where collision kernel is K(q, a)) = 2A1(q, aii/jrai. 

There are several important points which has to be discussed regarding the general structure of 
the kinetic equation, (i) The term trlZFX^id}, neglected in the ( 1262b , produces an effective vertex of 
the type tr{I[F]XlD^ZFX'i] after 'W integration, which indeed vanishes after X averaging, since 
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{X'lXl}x = 0. Thus, it indeed does not generate any additional terms into the collision integral, 
(ii) Throughout the derivation of the collision integral we persistently neglected all spatial drF^{r, t) 
and time drF^ir, r) derivatives of the distribution function, e.g. in ( I270l i. This is justified as long as 
there is a spatial scale at which F^Cr, t) changes slowly. In fact, gradients of the distribution, if kept 
explicitly, contribute to the elastic part of the colHsion integral 19611991 . (iii) We kept in the effective 
action only terms which are linear in the quantum component of the distribution function. There 
are, however, terms which are quadratic in Z^(r, t). These terms are responsible for the fluctuations 
in the distribution function and lead to the so-called stochastic kinetic equation or, equivalently, 
Boltzmann-Langevin kinetic theory 1461 IIOOIIIOIII . It was shown recently that Keldysh cr-model 
with retained Z^(r, r) terms is equivalent to the effective Boltzmann-Langevin description 1 1 021 
II03I . (iv) A collision integral similar to Eq. ( 12721 ) was derived within Keldysh cr-model formalism 
in l22ll . However, the S''^^ part of the effective action was overlooked and as a result, the obtained 
kernel of the collision integral turns out to be correct only in the universal limit f/^ ' — » 0. One 
finds from ( 12741 ) for t/^ ' — » that M(q, w) reduces to A1(q, co) - Im[£)*(q, w)], which is result 
of l22ll . (v) Finally, the present discussion can be generalized to include a spin degree of freedom. 
Corresponding kinetic equation and collision kernel were obtained in ||I04||I05l . 

7.6 Applications III: Interaction effects in disordered metals 
7.6.1 Zero-bias anomaly 

Having discussed in Section 16.61 several examples, where non-interacting version of the cr-model 
may be applied, we turn now to consideration of interaction effects. The first example of interest is 
the modification of the bare single particle density of states v of free electrons by Coulomb inter- 
actions. The question was addressed by Altshuler, Aronov and Lee 01061 II071 II08I . Although in 
their original work only leading order interaction correction was calculated, one may extent treat- 
ment of zero-bias anomaly beyond the perturbation theory l28l II 101 II III II 121 . Here we follow the 
sigma-model calculation of l22ll . 

We are interested in the single-particle Green function at coinciding spatial points 

g"\t - t') = -/«iA„(r, t)Mr, , (276) 

where {(. . .)) denotes both the quantum and disorder averaging. One may evaluate it introducing a 
corresponding source term into the action which is directly coupled to the bilinear combination of the 
fermion operators. Following the same algebra as in the Section|6| performing Keldysh rotation and 
disorder averaging, one finds that this source term enters into the logarithm in ( 1161b . Differentiating 
the latter with respect to the source and putting it to zero, one obtains for the Green's function 

-1 

sxp(iS[Q,^]). (277) 

H'.rr 

One evaluates the integral over the Q matrix in the saddle point approximation, neglecting both the 
massive and the massless fluctuations around the stationary point. Then, according to (12231 ). the 
pre-exponential factor is simply -invQ At the saddle point Q matrix is given by Eq. ( 12401 ). As a 
result, one obtains for ( 1277b the following representation 

g{t - t') = -inv J Dm exp {iTr{S'^ U^j,^6}) exp [i±(r, t)) A,_,, exp (-/±(r, f')) . (278) 
Since 3C is the linear functional of <1), given by (1235b . the remaining functional integral is Gaussian. 



g(t-t')^ jD[<5]exp(iTr{(5^f/o'M)) J ^Q] 



+ 2 + 0, 

2Tel 
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To calculate the latter, one rewrites phase factors of the gauge field aj*^ 

g±/3C"r° ^ 1 1^ g±,X3C''+3C") _^_ g±/(3C''-3C") j yl ^ 1 |^ g±/(3C''+3C') _ g±i(3C''-3C'') j ^ (279) 

Performing Gaussian integration in (1278b with the help of ( |279t . the result may be conveniently 
expressed in the form 

m = -inv 2 (r Ar /) B''^(f) , (280) 

where the auxiliary propagator ©"^(f) has the standard bosonic structure [as in, e.g., ( I252b 1 with 

B"<'*>(f) = /exp(/['y'^(f) -'V^(0)]/2)sinCy*('^'(f)/2), (281a) 

B^(f) = exp(/['V^(f) - 0^^(0)1/2) cos ([^^^(f) - ^"^(t)]/!) . (281b) 

The gauge fields propagator, 'Vir, t), defined by M52\ and M52>\ . enters ( 128 11 1 at coinciding spatial 
points 

(282) 



Mt)^ j ^exp(-/W)^4'(q,d^). 



Knowledge of the Green's function (|280t allows to determine the density of states according to the 
standard definition 



2n 



(283) 



In the thermal equilibrium, the Green's functions obey FDT (see ( |108t ), which together with the re- 
lations g'^ie) = 0^{e)+0''{£) and g^ie) = - exp(e/r)^^(e) allows to rewrite ( |283] l in the equivalent 
form 



v(e)= — ^"(e)[l+exp(-6/r)]. 
2n 



(284) 



Using ( |280t one relates greater (lesser) Green's functions to the corresponding components of 
the auxiliary propagators B^**^': 



>«)to>«) 



(0. 



The latter are found explicitly to be 



B>M(f)= -exp 



— coth — ( 1 - cos ) ± / sin J Im y (q, 



(285) 



(286) 



where we employed dSSTT i along with the bosonic FDT relations B*(0 - B'*(f) = B>(f) - B"=(f), and 
B^(f) = B>(f) + B^(f). Finally, combining ( 1284b and ( I285l l together, one finds for the density of 
states 



v(e) = 



J dtF,B'^ 



(f) exp(/ef) . 



(287) 



tanh(e/2r) 

Expanding (1286b to the first order in the interaction, 'V(q, a>), and substituting into (1287b . one recov- 
ers Altshuler and Aronov result for the zero-bias anomaly II106I . 



'^Equation )279t is based on the following property: consider an arbitrary function which is linear form in Pauli matrices 
f(a + hif), where a is some arbitrary number and b some vector The observation is that f(a + bo") = A+ Bcr, where A is some 
new number and B a new vector To see this, let us choose i-axis along the direction of the b vector Then the eigenvalues of 
the operator a + b& are a ± fo, and corresponding eigenvalues of the operator f{a + b(5") are f(a ± b). Thus one concludes that 
A=k ma + h) + fia - h)] and B = |; [/(a + fo) - f(a ~ b)]. 
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We restrict ourselves to the analysis of the non-perturbative result, (I286l l and ( I287l i. only at zero 
temperature. Noting that for T - 0, F, - (inty^ one obtains 



v(e) 



-I! 



sin \e\t 
df exp 



,co){l - cos tot) 



Jo ^ ^ 

— Im ) 'y*(q,w)sinwf 
^ V 



xcos 



In the two-dimensional case (I253I I with f/o = 2ne^/q leads to 
dw 

Jo 



(288) 



^yim[n/^q,a;)] 



1 - COS ClJt 

sin wf 



1 f ln(f/rei)ln(fTdw2) + 2Cln(too) ^289) 
87r2g \ 7rln(too) 



where gr = vD is the dimensionless conductance, wo = ^''c^, - 27re^v is the inverse Thomas-Fermi 
screening radius and C = 0.577... is the Euler constant. Since the fluctuations "W of the Q matrix 
were neglected, while calculating functional integral in (|277t , the obtained result ( |288t does not 
capture corrections, which are of the order of ~ gr"' ln(r/Tei) (in d - 2), see Section |7.6.2| Therefore, 



(I288l l can only be trusted for e not too small, such that {Sn^gy^ In(eTei)"' <*: 1, however, ln^(f/Tei) 
terms have been accounted correctly by the preceding procedure. If, in addition, gr"' ln((yoTei) 1, 
the time integral in ( 1288b may be performed by the stationary point method, resulting in 



v(e) = vexp 



1 



ln(|e|Tei) ' ln(Teiw5/|e|) 



(290) 



Thus, one achieved a non-perturbative resummation of anomalously divergent, oc In (exei), terms for 
a single-particle Green's function. The non-perturbative expression for the density of states essen- 
tially arises from the gauge non-invariance of the single-particle Green's function. The calculations 
above are in essence the Debye-Waller factor tl09J owing to the almost pure gauge fluctuations of 
electric potential, cf. ( 1278b . Gauge-invariant characteristics (such as conductivity, for example) do 
not carry phase factors, and therefore are not affected by the interactions on this level of accuracy 
(fluctuations of Q matrix should be retained, see next section). 



7.6.2 Altshuler-Aronov correction 

Here we consider yet another example where interactions are essential, namely electron-electron in- 
teractions correction (Jctaa to the Drude conductivity cro of the disordered metal II106II1071I108II . In 
contrast to the previous example, where density of states of an interacting disordered electron liquid 
was considered (Section r7.6.1l ), the correction to the conductivity is not affected by the interactions 
at the level of trial saddle point Qx = A and fluctuations "W must be retained. In what follows, 
we restrict our consideration to the lowest non-vanishing order in the expansion of the action ( 1244b 
over 'W, ( |187b and ( |188b . and identify those terms of the action which are responsible for interaction 
correction <5crAA- 

One starts from the part of the action S i[Qx,drOC] given by (1244b ). To the linear order in 
fluctuations "W one finds 

iSii'W, dr%] = -^Tr {[D5?3C«(Af"A - y") + (0)" - d,X"WA - Af')] w) , (291) 

where W ^tlo'Wo see ( fT85] l and (fT86] l. Note that in thermal equilibrium iS i ["W, drOC] = 0. 
Indeed, the expression in the square brackets on the right-hand side of ( 129 lb coincides with ( 1231b . 
which was used to determine the 3C[^] functional. In equilibrium it was possible to solve ( 12311 ) by 
an appropriate choice of 3C[<I>], see ( 1236b . This was precisely the motivation behind looking for the 
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Figure 14: Diagrammatic representation of an effective four-leg vertex Vaa, see (I296l l. which gen- 
erates Altshuler-Aronov correction to the conductivity. 

saddle point for each realization of the field O to cancel terms linear in 'W. Since it was not possible 
to find the exact saddle point, such terms do appear, however, only in the second order in drX. These 
latter terms originate from the S2[Qoc, drOC] part of the action. Expanding (1244b ) to the linear order 
in 'W one finds 

iS2['W,drX] = ^ Tr{5r3C"(e, - ej) [TKfK " KrKf] W,,,,dr%^e2 - 63)} 

= nvDTx[d,%'{ei - eo) [Mid,^,, + Ml^^ J,^,]drX{e2 - e^)] , (292) 

where we used notation Xj - (X'^', X''), and introduced coupling matrices between diffusons {d, d} 
and the gauge fields JC^'^*' 

/O \ , J2F,-F,-F, l^F,F,-2F,F, \ 

'•^^^ ^0 -2F,J' \ -l-F,,F,,+2F,,_F,, F,^+F,,-2F,,F,,_F,J- ^^^^> 

Employing now the general expression for the conductivity ( 12011 ). we show that Altshuler-Aronov 
interaction correction to the conductivity dcTAA is obtained from ( 12921 ) 

= exp(/..[nV,.,:K]))^^ , (294) 

where the averaging goes over the diffusive modes as well as over the fluctuations of the electric 
potential. Note also that as compared to ( 1201b here we perform differentiation over drX and not the 
vector potential A itself. The two definitions are the same since the vector potential and the gauge 
field enter the action ( 1241b in the gauge invariant combination (1227b . 

Having ( 1187b and ( |292b we deal with a Gaussian theory of the diffuson modes d and d fluctua- 
tions, which allows for a straightforward averaging in ( 1294b . Integrating over the diffuson modes, 
one finds 

(exp(/52[^, drX]))^ = exp (/Vaa[3C]) . (295) 
This way the (drX^ effective four-gauge-field vertex is generated 

Vaa[3^] = 47rvD-Tr { F^F^ - F,, - F,,)drX\r, ei - e2)drXHr, £2 - €3) 

X £)*(r - r', 63 - ei)dr'X''{r', 63 - £4)dr'X'='(r', 64 - ei)} . (296) 

Its diagrammatic representation is depicted in Figure[T4l This vertex originates from TT{drXM'^ddrX} 
and Tr{drXM''ddrX} parts of the action ( 12921 ) after we pair d and d by the diffuson propagator 
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Figure 15: Diagrams for the interaction correction to the conductivity 6a aa- These diagrams are 
constructed from the effective vertex Vaa[3C] by keeping one classical and one quantum leg to be 
external, while connecting the remaining two by the interaction propagator ^^^(q, w). 

{dd)<w <^ ly^- The factor F^^ originates from q- q element of the matrix M'', while the combination 

296i is the cl - cl element of the matrix M''. By 



2F, 



Ff, of the distribution functions in 



writing Vaa[3^] in the form of ( |296t we kept only contributions with the lowest possible number 
of quantum gauge fields drOC^- However, matrix M'' has all four non-zero elements, thus Vaa[3<^] 
in principle also contains contributions with four and three legs carrying the quantum gauge fields. 
The latter are to be employed in calculations of the corresponding interactions corrections to the 
shot-noise power, see jllSl for details. 

Having performed "W averaging, one brings now Yaa [3^] into (|294| i and integrates out X field. 
For the conductivity correction this gives 



2 ff ^(^-. + ^-x^-^^.-^^^-)^(q'^)(^r3c^'(q'^)^r3<:^'(-q'-^^^ 

(297) 

where new integration variables e = (63 + ei)/2 and w = 63 - ej were introduced. The 3C averaging 
produces two diagrams. Fig. [15] for (5cr aa, which follows naturally from the effective vertex shown 
in Figure [141 ^.ftei one pairs two external legs by the interaction propagator. In the universal limit of 
strong interactions {/^ ' — > the propagator ^*(q, co) takes the simple form. As a result. 



1 



2vD Dq^ - id) ' 



(298) 



which follows from (1251b and (1253b . Inserting ( |298b into Eq. ( |297b and carrying e integration one 
finds 



60-A 



ctd 



--y 



a> coth — 
2T 



1 



{Dq^ - icoY 



(299) 



In two dimensions this expression leads to the logarithmically divergent negative correction to the 
conductivity: ^ctaa - ln(l/rrei), where the elastic scattering rate t^j' enters as an upper cutoff 
in the integral over the frequency oj. A detailed review of the eff'ects of the interaction corrections 
on disordered conductors can be found in 11081 , see also 



7.6.3 Relaxation rate 

Kinetic equation discussed in Section [73] ma v be used to find energy relaxation rate Il98]|108llll4l 



II 151 . Focusing on the "out" term of the collision integral in (1275b . one may introduce the out 
relaxation rate for an electron of energy e, as 



1 



Tout(e) 



dwde' K(q, oj) n/r (e)[l - iif (e - w)]n/r(e')[l - n/r(e + w)] , 



(300) 
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where all electron distributions were substituted by Fermi functions. This is appropriate if one 
is interested in small (linear) deviations of from its equilibrium value nf{e). Equation ( 13001) 
simplifies considerably at zero temperature, T = 0. Indeed, Fermi distribution functions Umit energy 
integration to two ranges -o) < e' < and < w < e, where the product of all occupation numbers 
is just unity. In the universal limit of strong interactions, f/g' — > 0, the kernel acquires a form, see 
Eq. ( 12741 ) 

K(q,w) = --— ^. (301) 
Inserting K(q, to) into (1300b . one finds for the out relaxation rate the following expression 

1 4 r''' . r° . , 1 ki 



.,(e) 



±y fd^fd.- ^\ ,.J±, (302) 



where g - vD and momentum integral was performed for the two-dimensional case. For an arbitrary 
dimensionality d, out rate scales with energy as t „„'t(e) oc {\lvd){elDYI^, see IJOSil for further details. 



7.6.4 Third order drag effect 

Discussing Coulomb drag in Section [533] it was emphasized that the effect appears already in the 
second order in inter-circuit interactions and the particle-hole asymmetry is crucial. In the linear re- 
sponse at small temperatures the drag conductance appears to be quadratic in temperature, see ( I155l l. 
Here we discuss the third order in the inter-layer interaction contribution to the drag conductance. 
Although, being subleading in the interaction strength, it does not rely on the electron-hole asym- 
metry (in bulk systems the latter is due to the curvature of dispersion relation near the Fermi energy 
and thus very small). We show that such a third order drag is temperature independent and thus 
may be a dominant effect at small enough temperatures 11 161 . Technically the third-order contribu- 
tions originate from the four-leg vertices (see Figure [14] and corresponding ( 12961 )), which describe 
induced non-linear interactions of electromagnetic fields through excitations of electron-hole pairs 
in each of the layers. 

Following 111 161 we consider two-dimensional electron gas bilayer and apply NLSM to calculate 
the drag conductivity. From the general expression ( 12011 ) with the help of ( 12961 ) one defines drag 
conductivity as 

-2 1 / Waa[3C] Waa[3C] \ 



'"''-^ = - 2 n^l, Q U(5r3Cf (Q)) <5(5r^C^(-^2))/^ ' ^^^^^ 

where indices 1 , 2 refer to the drive and dragged layers, respectively, following notations of Sec- 
tion 15.5.31 The averaging over the fluctuating gauge field % is performed with the help of the 
correlation function 



n/«,(q, a>) = 2i{X'M, u)Xl(-<i, -oj))^ = / / /! , . , (304) 

(Dfl^-^ - iLo){Dt,q~ - icj) 



where a,b — (1,2) and [/^^(q, w) is 2 x 2 matrix of retarded screened intra- and inter-layer in- 
teractions calculated within RPA. It is a solution of the following matrix Dyson equation, f/^ - 
Uo + UqTI'^U'^, where 



v\D\q Q 



(305) 



Off-diagonal components of t/o matrix represent bare Coulomb interaction between the layers, 
where d is the inter-layer spacing. Note also that the polarization operator matrix n*(q, oS) is diag- 
onal, reflecting the absence of tunneling between the layers. 
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q/2+q',a}/2+w' 




Figure 16; Two diagrams for the drag conductivity crjiag in the third order in the inter-layer in- 
teractions, 'V^jCq, w), denoted by wavy lines. The intra-layer diffusion propagators £)*(q, w) = 
{Daq^ - iw) ' are denoted by ladders. 



We are now on the position to evaluate the third-order drag conductivity. Inserting ( |296l l into ( 1303b 
and performing averaging with the help of ( 13051 ). one finds the following expression for drag con- 
ductivity 



o"drag = 32e Tvi 



dwdw' 
4;r2 



"H'i(w,w')'H2(w,w') 



: 2 Im [Sf (q, c.)l)«(q, a;)^f,(q, (f " Q'' f " -'l^n (f + Q', f + 

9.9' 



(306) 



The two functions '7Yi(w, Q) and '7Y2(w, Q) originate from the integration over the fast electronic 
energy e. Figure [15] in the active and passive layers, respectively. In the dc limit they are given by 



"Hiico, w') = 2 - + w/2) - - w/2) + , 
<H2(w, w') = r — H- w/2) - - w/2)] , 



„ (X) 1 iJLl \ 

SM.-coth(-). 

The corresponding diagrams are constructed from the two vertices of Figure [141 one for each of the 
layers, see Figure[T6l It turns out that there are only two ways to connect them, using the propagators 
^«i(q,w), since (JC^JOC*) = 0. 

In the following we assume identical layers and consider the experimentally most relevant case 
of the long-ranged coupling, where /cc/ » 1 . Here k - 2ne~v is the Thomas-Fermi inverse screening 
radius. In this limit the effective interlayer interaction potential, Eqs. ( 13041 ). acquires a simple form 



(307a) 
(307b) 
(307c) 



g KdDq- - lico ' 



(308) 



where g = vD. Next, we substitute D'^iq, a>) along with ( 1307b and (1308b into (1306b and perform the 
energy and momentum integrations. Inspection of the integrals shows that both energies w and w' 
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are of the order of the temperature u ~ to' ~ T . On the other hand, the characteristic value of the 
transferred momenta \s q ~ q' ~ yjT HDkcI) «: y/T/D, cf. (I308l l. Therefore, we may disregard Dq^ 
in comparison with ico in the expressions for £)f (q, cl>), approximating the product T^D^ in ( 13061 ) 
by -Lo^^. Such a scale separation implies that the four-leg vertices are effectively spatially local, 
while the three inter-layer interaction Unes are long-ranged. 

Rescaling energies by T and momenta by -\/T/{DKd), one may reduce expression (1306b for the 
drag conductivity to o-±-i^g = g^^ {kcIY^x (dimensionless integral). The latter integral does not 
contain any parameters, and may be evaluated numerically I116J . In the limit CTjjiag glRg the drag 
resistance pdrag is given by pdrag - o-dmgRg/g^, resulting finally in pdrag ~ 0. 27 R gg^^ {I /kcI)^. This 
is the temperature-independent drag resistivity, which may be larger than the second order (in the 
inter-layer interactions) contribution. The latter goes to zero at small temperatures as T^. Further 
details and discussions can be found in 111 161. 



8 Superconducting correlations 

8.1 Generalization of the cr-model 

So far we have been discussing the unitary version of Keldysh cr-model, i.e. the one, where the 
time-reversal symmetry was supposed to be broken by, e.g., external magnetic field. We now switch 
to the orthogonal symmetry class, with the unbroken time-reversal invariance. The case in point is 
superconducting fluctuations in disordered metals. The Keldysh sigma-model, generalized for the 
disordered superconductors was developed by Feigel'man, Larkin and Skvortsov Ii24l|l 171 . It is also 
applicable for treating weak-localization effects in normal metals. 

We proceed to describe disordered superconductors by adding the BCS term to the Hamiltonian 
of a metal 



^Bcs = -^ j driA;"(r)iAj(r)iAi(r)iAT(r), 



which corresponds to the short-range attraction in the particle-particle (Cooper) channel mediated 
by electron-phonon interactions, where A is dimensionless coupling constant. In a standard way 
^BCS translates into the Keldysh action 



'BCS 



^^J^dtjdr ^^(r, f)<Ai(r, OiAi(r, t)tf/^{r, t) , 



where the time integral is calculated along the Keldysh contour This four-fermion interaction term 
may be decoupled via Hubbard-Stratonovich transformation, by introducing an auxiliary functional 
integral over the complex field A(r, f): 



exp(i5BCs) = J"D[A]exp^; J" dx -j\A(x)f + A(x)if/^(x)4fi(x) + A*(x)iffi(x)i//f(x) 



(309) 



here x = (r, t) and J dx - J^dt J dr. To make further notations compact it is convenient to introduce 

a bispinor fermionic vectors ^' — 1/ V2(i^|, i/'j, tf/i, -i^i)^ and ^'^ - 1/ V2(i^|, tf/i, defined 
in the four-dimensional space Q, which can be viewed as the direct product S ® T of the spin 
{^■\,4'[) and time-reversal spaces {(//,ijf)- In principle, choice of the bispinors is not unique. One 
can rearrange components of the bispinors in a different manner, separating explicitly the Gor'kov- 
Nambu ffll81|119l (A^) (i^i, i^^) and spin spaces. Finally one may equally think of ^ as acting in the 
direct product of the Nambu and time-reversal subspaces. These three representations are equivalent 
Q. = S ®T cc N®S cc N®T and the choice between them is dictated by convenience in calculations 
for a particular problem at hand. In most cases we use N ® S choice and omit the S part, since 
the theory is diagonal in spin subspace. Vectors ^' and are not independent and related to each 
other - (C^P)^, by the charge-conjugation matrix C = ify ® Sv, where f, and i,, for / - 0, x, y, z, 
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are Pauli matrices acting in the Nambu and spin subspaces, respectively; matrices, as before, act 
in the Keldysh sub-space. To avoid confusions, we shall specify, where appropriate, Keldysh and 
Nambu sub-spaces by subscripts K and correspondingly. 

After the Hubbard-Stratonovich transformation ( 13091 ). along with the standard treatment of dis- 
order and Coulomb interactions, the action appears to be quadratic in fermion operators. Performing 
thus Gaussian Grassmann integration, one obtains for the disordered averaged partition function 



Z = J D[(D,A] exp(^Tr{6f/oit<i)}-^Tr{A"'"tA)) J \i\Q\ exp(/5 [g. A, A, O]) , 

(310) 



!5[e,A,A,(D] = --^Trle^l + Trln 
4rei 



+ 7-^2 + <i) + Vf 5A + A 



2rer 



which generalizes ( 1222b . In the last equation and throughout the rest of this chapter we use the check 
symbol O to denote 4x4 matrices acting in the K®N space, while hat symbol O for the 2x2 matrices 
acting in Nambu and Keldysh subspaces. Equation ( 13101 ) contains matrices T = (Xiigifo, S = (To®t-., 
G"' = iRd, + d\j2m + yu, and matrix fields 

<E)(r, f) = W\x, t)&o + <l)«(r, t)&,] ® tq , A(r, f) = [A"'(r, f)<3-o + A^r, t)&,] ® tq , 

A(r, t) = [A^'(r, t)&o + A^r, t)&,] ® f + - [A*''(r, t)&o + A'-'^ir, t)&,] ® t_ , (311) 

with f± = (r.v + ity)/2; Q matrix also has 4x4 structure in Keldysh and Nambu spaces along with 
the matrix structure in the time domain. 

We next perform the gauge transformation in ( 13101 ) with the help of ^^^'^^'(r, f) fields, as ir0 ( 12251 ). 
and expand the logarithm under the trace operation in gradients of Qx matrix (similar to the calcu- 
lation presented in Section|6|. As a result, one obtains the action of disordered superconductors in 
the following form 

S[Q,A,A,^]^SA + S^ + Sa, (312a) 
Sa = -^Tr{A;;,tA3c) , 5o - ^TrlOjctOjcl , (312b) 

Sa ^'■^Tr{D{drQxf-4Zd,Qoc+4iijcQx+4iAjcQ3c]. (312c) 

Here gauged electromagnetic potentials and Ajc are related to the bare ones <1) and A by ( |227T i. 
while the gauged order parameter field is given by 

A3c(r,f) = exp(-iM(r,f))A(r,f)exp(/3(:(r,f):§). (313) 

As compared with ( 1242b the covariant spatial derivative in ( 1312b ) contains an extra H matrix due to 
Nambu structure, i.e. 

drQjc = drQoc - &0C, Qx] ■ (314) 

Varying the action ( 1312b with respect to Qx, under the constraint - 1, yields the saddle 
point equation 

dr(DQxodrQx)-\^d„Qx]^+i[^x+Ax,Qx]^0, (315) 



which for 3C = coincides with the dynamic Usadel equation 11391 . The classical solution of this 
equation is to be sought in the form 

with retarded, advanced and Keldysh components being matrices in Nambu subspace. 



'*In the superconducting case the gauge transformation contains phase factors exp{±iZOC), which is different from {225) 
by an extra matrix H in the exponential. 
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Varying the action with respect to the quantum component A*''(r, t) of the order parameter field, 
one finds the self-consistency equation for the classical component of the order parameter 

A^^(r, t) = 7r^Tr{(,T, ® f-)Qx} . (317) 

Finally, varying the action with respect to the quantum components and A'' of the electromag- 
netic potentials one obtains set of Maxwell equations, which together with the dynamic Usadel equa- 
tion (1315b and self-consistency condition (1317b represent the closed system of equations governing 
dynamics of the superconductor 

In the generalized cr-model action ( 1312b . and subsequent dynamical equations for Qt,'{r) and 
A(r, f), all the relevant low-energy excitations have been kept indiscriminately. The price one pays 
for this is the technical complexity of the theory. In many practical cases this exhaustive description 
is excessive and the theory may be significantly simplified. For example, one often considers a 
superconductor in the deep superconducting state T ^ Q, with well defined gap |A|, and studies 
dynamical responses when perturbing frequency at of the external field is small w <«; |A|, thus dealing 
with the quasi-stationary conditions. For this case quasiclassical kinetic equations of superconductor 
can be derived from ( 1315b . As an alternative, one may consider temperature range in the vicinity of 
the transition \T - Td <*c T^, where the order parameter is small |A| <sc T^, and develop an effective 
theory of the A(r, t) dynamics, i.e Ginzburg-Landau theory. Both approximations follow naturally 
from the general cr-model theory and will be considered in the next sections. 



8.2 Quasiclassical approximation 

In the superconducting state, choosing an optimal gauge field 3C(r, e) that is valid in the whole energy 
range is a complicated task. However, it had been shown in the 01201 that in the deep subgap limit 
(e <& |A|) the effect of the electric potential on the quasiclassical Green's function Q is small in the 
parameter e/|A| «: 1 and hence as an approximation one may set 3C(r, e) = 0. This assumption will 
be used below I '1 

In a spatially uniform, equilibrium superconductor the saddle point Usadel equation is solved by 
the the following Q-matrix 

e«'^>(e) = ±— J==( ' M , (318) 
V(e + /0)2 - |A|2 I -f In 

while Q'^ - tanh ^fiQ'^ ~ Q^)- We have suppressed superscript cl, writing the order parameter as 
A (its quantum component will not appear within this section). Substituting (13181) into the self- 
consistency condition ( 13171 ). one obtains the standard BCS gap equation 

de e 

A = AA , tanh — , (319) 

which has a non-zero solution for |A| below a critical temperature T^- 

In presence of boundaries or proximity to a normal metal one faces the problem of spatially 
non-uniform superconductivity. In this case, both A and g*^'*' acquire a coordinate dependence and 
one should look for a solution of ( 13151 ) and ( 13171 ). In doing so, we will assume that Qtf is static, i.e. 
independent of the central time and pass to the Wigner transform representation. From the retarded 
block of the 4x4 matrix Usadel equation at (D = and A = we obtain 

dr{D Q%Q'') + i€[f„ Q"] + i[A, 2^] = . (320) 

With the similar equation for the advanced block of the matrix Usadel Equation ( 1315b . The Keldysh 
sector provides another equation, which is 

driD Q^drQ'' + D Q^drO^) + ie[f„ e^] + i[A, Q"] ^ . (321) 



17 



Within this section the subscript 3C is suppressed in the notations of Qx matrix, Qx — > Q, and all other gauged fields. 



72 



The non-linear constraint ff' - I imposes the following conditions 

= = e«e^ + e^e^ = o. 022) 

They may be explicitly resolved by the angular parametrization M121I for the retarded and advanced 
blocks of the Green's function matrix: 

e«(r,e) = ( . ,7^^. . . ^i'^^^'^Pj'^) ) , (323a) 
^ ^ ^ y -sinh6exp{-ix) -cosh 6* ^ ^ 

e^(r,e) = -f-[e^]^r- = ( .,7°'''/.-, -^i"h0expOB\ (323^) 

where 0{r, e) and;if(r, e) are complex, coordinate- and energy-dependent scalar functions. As to the 
Keldysh component, it can be always chosen as 

Q'^^Q^oF-FoQ^, (324) 

where F may be thought of as a generalized matrix distribution function. Following Schmidt- 
Schon 11221 . and Larkin-Ovchinnikov 012311 we choose 

Ar, .) = ( I '> ^^^^ » ^^^^ I - Far. . Frir. , (325) 

where abbreviations Fkt) refer to the longitudinal and transverse components of the distribution 
function with respect to the order parameter Physically Ft corresponds to the charge mode of 
the system and determines the electric current density, while Fl corresponds to the energy mode, 
determining the heat (energy) current (further discussions may be found in books of Tinkham 112411 
and Kopnin |fT25l ). 

Substituting in the form of (1323b into ( 1320b . one finds from the diagonal elements of the 
corresponding matrix equation 

D dri sinh^ drx) = 2/| A| sinh 6* sin(<^ - x) , (326) 

where the order parameter is parameterized as A(r) = |A(r)| exp{!^(r)}. From the off-diagonal block 
of the matrix Equation ( 1320b . using ( 1326b . one obtains 

Z) 5^6* + 2/esinh6'-2/|A| cosh 0cos(^-;t') = ^(drxf sinhlO . (327) 

We proceed with the equation for the Keldysh component of the Green's function matrix Q'^ . Using 
decomposition (13241 ) and substituting it into Eq. ( 13211 ). one obtains 

D [dlP + Q^drQ^drP - drPQ^'drQ'' - driQ^d^F Q"")) + ie (e^[t-, F] - [f„ F]Q') 

+i (e«[A, F] - [A, F]Q^) = . (328) 

Now using ( 1325b for F and: (i) taking Nambu trace of the above matrix equation; (ii) multiplying 
the above equation by fj and then tracing it; one finds two coupled kinetic equations for the non- 
equilibrium distribution junctions Fl^t), which can be written in the form of conservation laws Ml 261 

drCDAFi - DdrFj Y) + DdrFjJs = rU ' (329a) 
drCDrdrFT + DdrFLY) + DdrF^Js = J*„ii . (329b) 
Here we have introduced energy- and coordinate-dependent diff'usion coefficients 

DL(r,e) = ^Trjfo - Q''&']^ = y [l + |cosh0|2 - | sinh^l^ cosh(2ImM)] , (330a) 
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2)Kr,e) = jTr{fo-f,e^f,e^}^ = - [l + | cosh 01^ + | sinh^l^ cosh (2ImM)l , (330b) 
density of the supercurrent carrying states 

Js(r, e) = ^TrlUQ^drQ^ - ^d.Q")]^ = -Im (sinh^ 5,;^) , (331) 
and the spectral function 

F(r,e) = ]tx[qH,Q^]^ = i| sinh^l^ sinh(21mM) ■ (332) 
Finally, the right hand side of ( |329t contains the collision integrals 

^cou = ^ Tr {f,(e«A + W)]^ = 2Fr|A|Re [sinh sin(^ - , (333a) 

-^cou = ^Tr{e«A + Ag^}^ = -2FHA|Im[sinh0cos(^-;^)] . (333b) 

Collision integrals associated with the inelastic electron-electron and electron-phonon interactions 
are not discussed here, one may find corresponding derivations in the book of Kopnin 11251 . Equa- 
tions ( I326l l, (13271 ) and ( |329t , together with the spectral quantities ( I330b - (l333b represent a complete 
set of kinetic equations for disordered superconductors applicable within quasi-classical approxi- 
mation. These equations are supplemented by the self-consistency relation, see (13171 ) 

A(r) = ^ J"de|[sinh0exp(;^) + sinh 0exp(!]j^)]F/, - [ sinh exp(i^) - sinh 0exp(!j)]F7-| , (334) 

and the boundary conditions for the Green's functions, expressing the current continuity 1170111271 
[1281 [T29l . 

(TL^LQidrQi = (TRJiRQRdrQR = griQi, Qr\ , (335) 
where cr and J[ are the bulk normal-state conductivity and the cross section of the wire next to the 
interface, LjR denote left/right from the interface, respectively, and gj is the interface tunneling 
conductance. 

An analytic solution of the system of kinetic equations ( I326l l-( l329l l is rarely possible. In general, 
one has to rely on numerical methods. To find solution for a given transport problem, one should 
proceed as follows 11211 . 

1. Start with a certain A(r). Usually one takes A - const everywhere in the superconductors and 
A = in the normal metals. 

2. Solve Usadel equations ( I326b -( [327b for the retarded Green function, thus determining spectral 
angles 6{r, e) and;t'(r, e). 

3. Use these solutions to calculate spectral kinetic quantities Dij{r, e), J's (r, e) and Y{r, e). 

4. Solve kinetic equations (1329b for /^L/r(r, e). 

5. Calculate new A(r) from Equation ( 1334b . and iterate this procedure until the self-consistency 
is achieved. 

Having solved the kinetic equations one may determine physical quantities of interest. For example, 
for the electric current one finds j - j„ + jj, where j„(r) - J deX)r(r, e)(9r/^r(r, e) is the normal 
component and j ^(r) - vD J deFiir, e)J7s (r, e) is the supercurrent density. 

The quasiclassical kinetic theory of disordered superconductors, outlined above, may be applied 
to study various phenomena. To name a few: the proximity related problems in the superconductor- 
normal metal heterostructures 11301 11311 11321 [1331 . non-equilibrium Josephson effect II1341 11351 . 
Hall effect 11361 , thermoelectric phenomena 1137111381 in superconductors, shot noise II139I , engi- 
neering of non-equilibrium distribution functions 11401 and many other problems may be success- 
fully tackled with the help of ( I326b -( [329l l. Several relatively simple (equilibrium) examples are 
considered in Section |8]4|for illustration. 
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8.3 Time dependent Ginzburg-Landau theory 



Gor'kov [141| had shown that the phenomenological Ginzburg-Landau (GL) theory HI 421 follows 
naturally from the microscopic BCS model in the limit when temperature is close to the critical 
one \T - Tc\ <*: Tc- Later Gor'kov and Eliashberg II143I extended application of the Ginzburg- 
Landau theory to include time dependent dynamical phenomena. It was revisited in a number of 
subsequent publications Qli [l45l [l46l [1471 [Mi [1491 [B^ and books IHl [l25l [BD . Within the 
cr-model terminology the static GL functional may be obtained by means of supersymmetric II152II 
or replica 11531 approaches. Here we discuss the dynamic theory in Keldysh formulation 11541 . 

The way dynamical time dependent Ginzburg-Landau (TDGL) theory is derived from (13121) 
allows to formulate it in terms of the effective action, rather than the equation for the order parameter 
only, as it is done in a traditional way. As a result, in addition to the average quantities one has an 
access to fluctuation effects, since TDGL action contains the stochastic noise term, which serves to 
satisfy the fluctuation-dissipation theorem. Moreover, one may naturally and unmistakably identify 
an anomalous Gor'kov-Eliashberg (GE) term II143L which preserves gauge invariance of the theory, 
along with the Aslamazov-Larkin (AL) II155II , Maki-Thompson (MT) II 1561 and density of states 
(DOS) terms II162I . which renormaUze the conductivity and single particle density of states owing 
to superconductive fluctuations. Although Aslamazov-Larkin term is correctly captured by most of 
the approaches to TDGL equation, Gor'kov-Eliashberg, Maki-Thompson and DOS are frequently 
lost in many works on TDGL. 

The strategy of deriving the effective TDGL theory starting from the general cr-model ac- 
tion ( 1312b is as follows, (i) Choose a parametrization of a saddle point Q matrix manifold, which 
resolves the non-linear constraint Q' - I. (ii) Integrates out Gaussian fluctuations around the saddle 
point and (iii) keeps terms up to the second order in all quantum fields (the order parameter A and 
electromagnetic potentials O and A) in the resulting action, (iv) Rely on the assumption that the 
electronic system is always in a local thermal equilibrium. This in turn implies that the external 
fields are not too large. More precisely, the electric field E is such that e|E|fo ^ T^, while the mag- 
netic field H is restricted by the condition e|H|^o ^ l/^o, where - ^DjTc is superconductive 
coherence length. The restrictions on spatial and temporal scales of the external fields along with 
the fact that electrons are in local equilibrium considerably simplify the theory. In particular, most 
of the terms in the effective action acquire a local form in space and time. Nevertheless, the effective 
theory does not take a completely local form. 

This procedure is relatively straightforward in the case of gapless superconductivity. The latter 
occurs either in the presence of magnetic impurities, or in the fluctuating regime above the critical 
temperature T > T^. In the gapped phases, T < Tc, the situation becomes more complicated. As 
noted by Gor'kov and Eliashberg |143|, the difficulty stems from the singularity of the BCS density 
of states at the gap edge. The latter leads to a slowly decaying oscillatory response at frequency 2A/^ 
in the time domain. As a result, the expansion in powers of the small parameter A/r^- <*: 1 fails. In 
principle, it may be augmented by an expansion in A/(ha)), in case of high-frequency external fields. 
To describe low-frequency responses in the gapped phase, one needs a time non-local version of 
the TDGL theory. The analysis is greatly simplified in the presence of a pair-breaking mechanism, 
such as magnetic impurities or energy relaxation. Such a mechanism may eliminate singularity in 
the density of states, leading to gapless phase in the presence of finite A. Under these conditions, an 
expansion in powers of Ar^/ft «: 1 and llit^ <*: 1 is justified and thus a time-local TDGL equation 
may be recovered (here is the pair-breaking time). Within this section only fluctuating regime, 
T > Tc, will be considered. In this case the spectrum is gapless automatically and there is no need 
in an expUcit pair-breaking mechanism. 

Proceeding along the steps (i)-(iv), outlined above, one recalls that at T > Tc energy gap self- 
consistency Equation (1317b has only trivial solution with (A^') = 0. Thus the trial saddle point of 
the action (1312b collapses back to the metallic state Qjc = A - A® f~, where A is defined by ( 1165b . 
The Gaussian integration around this Qx includes Cooper modes, which are accounted for in the 
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following parametrization of Qx-matrix: 



with the following choice of the fluctuation matrix 



c„, (r)f + - 4, (r)f _ d,,' (r)f + c/= , (r)f , 
£/„' (r)fo + d^, (r)f- c,,' (r)f + - c*,,, (r)f _ 



o^l- 



1/ = 1/ ' = 



(336) 



-1 



5 To . 

(337) 

As compared with ( 1186b . ^ contains twice as many diffusive modes, which are described by 
four Hermitian matrices in time subspace: {d,d] and {d~,d^]. It also contains the Cooper modes 
described by two independent complex matrix fields {c, c). Now substitutes the 'VV-dependent 
Qx matrix QjcV^V] into ( I312l i and expands the action up to the second order in 'W fluctuations: 
S {Q, A, A, 5 YW, A, A, <I>]. After this step the Gaussian integration over is possible (see 

details of this procedure in AppendixlDt 



(338) 



\)YW] exp [iS YW, A, A, (£]) = exp (/5efF[A, A, <D]) , 

which leads eventually to the effective TDGL action. It consists of several contributions: 

5eif[A,A,<5] = 5n[A,0] +5gl[A,A,0] +5sc[A,A,(D] + 5mt[A,A,(D] + 5dos[A,A,<1)], (339) 

which we describe in order. 

The action 5n[A, <i>] is the normal metal part of ( 13121 ). which is obtained from S [Q, A, A, (t>] by 
setting Qjc - A and A = 0. It reads as[3 



5n[A,<D] = e\DTT{A 



3C 



03'--% ' 



(340) 



where arrows on top of the time derivative indicate direction of differentiation. Since our starting 
point is the normal saddle point ( I165l l. !}C[<1>] functional is given by Eq. ( 12331 ) and gauged vector 
potential Ajc is defined by Eq. ( 1227b . 

The ScL is the time dependent Ginzburg-Landau part of the action 



5gl[A,A,(D] = 2vTT[A^^(r,t)L-'Aoc(r,t)) 



(341) 



which governs time and space variations of the order parameter under the influence of external 
potentials. The effective propagator L ' has the typical bosonic structure in the Keldysh space 







Ik 



with the components given by 



n 



+5, - + D{dr - lieh'^f - ^|A5^|2 



(342) 



(343a) 



coth—[L-^\co)-L-^\cj)] ^ - , 



(343b) 



2r' " /I V 'J 2 

where u) <^ T ^ Tc and Ginzburg-Landau relaxation time is defined as tgl - ^/^{T - Tc). Note 
that under the assumption T - Tc ^ Tc, GL part of the action acquires a time-local form. 



'^Note that in )340) and throughout the rest of this section we have restored electron charge e accompanying source fields 
A — > eA and d) — > e<S>, such that A and d) are now actual electromagnetic potentials, see Note 8. 
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The S sc part of the action is responsible for the super-current 



5sc[A,A,0] = — Tr{A*3,Im[A-'(5r-2/eA5^)A5^]} . 



(344) 



The abbreviation is due to the fact that S sc, being differentiated with respect to A'', provides standard 
expression for the super-current in terms of the order parameter 11241 . 

The Maki-Thompson part of the action, Smt, is responsible for renormalization of the diffusion 
coefficient in the normal action S'n owing to the superconductive fluctuations. It reads as 



6D 



(345) 



(346) 



5mt[A,A,<5] = e\Tx[ll^{r,t)tsD(.t,t')l3c(X,t')] , 

where the operator Tsoit, t') is given by 

-d^6DfJ, 
-6D^,,t 2iT{SD^J,,+6Dfl) J/ 

The diffusion coefficient correction 5D^"^{A'x\ is the non-local functional of the fluctuating order 
parameter 

5D^, = ^ J dr'dr" C-^A^' (r', r) A5^(r", r) C;;- , (347) 

where t = (f + f')/2. The retarded C"\, ~ 0{t - f) and advanced Cl''], ~ 9{t' - t) Cooperon 
propagators are Green's functions of the following equations: 

d, - /e05|-(r, T+) + /e05|-(r, t_) -D[dr- ieA'^ir, t+) - ieA'^ir, t_)]'} C^f,, = 6r-r'S,-, , 

(348a) 

{-5, + /e05i(r, r+) - ie<:>'^(r, t_) -D[dr- ieA'^(r, t+) - ieA'^jcir, t-)]'} C%, = 6r-r'Sr-,' , 

(348b) 

with T± - T + t/2. Note that MT action (1345b . has exactly the same structure as the normal action 
5n- It therefore can be incorporated into (1340b by adding a renormalization of the normal diffusion 
constant D6,-,' — > DS,-i' + SD^J^, that is non-local in time. 
Finally, ^dos has similar structure to Smy in (1345b 



I DOS 



[A,A,(D] ^e^UYx 



X5,(r,f) 



-d, 
-'d, 4-iT 



A3c(r, f) 



with locally renormalized density of states 



6v 



,DOS 



-v^|A5^(r,f)P. 



(349) 



(350) 



Each term of the effective action ( 13391 ) admits a transparent diagrammatic representation, shown in 
Figure [TtI 

An equivalent way to display the same information, which is encoded in the effective action (13391 ), 
is to use the set of stochastic time dependent Ginzburg-Landau equations. To derive those one 
needs to get rid of terms quadratic in quantum components of the fields: A^ in Sch, and A^ in 
5n + 5mt + Sdos- For the first one, this is achieved with the Hubbard-Stratonovich transformation 

exp(-^Tr{|A^^|2)) = j D[^a] exp(-|^Tr{^ _ /^^A^c " '^aA^^}) ■ (351) 

As a result, the effective action S^ff in ( 1339b acquires the form linear in quantum components of 
the order parameter. Integration over the latter leads to the functional delta-function, imposing the 
stochastic equation of motion. This way the TDGL equation is derived 



n 1 



A5^(r,r) = ^A(r,r). 



(352) 
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A(r,0 



Figure 17: Diagrammatic representation of the effective action 5eff[A, A, O]. a) Conventional 
Ginzbirg-Landau functional Sgl, see ( 1341b . b) Anomalous Gor'kov-Eliashberg coupling between 
the scalar potential and the order parameter (see ( I357l i and discussions below), c) Paramagnetic and 
d) diamagnetic parts of the super-current action Ssc- e) Local DOS term S'dos- f) Nonlocal MT 
term Smt- In the case of diagrams e) and f) there are two possible choices for the vector poten- 
tials: classical-quantum, which is a part of the current, and quantum-quantum, which is its FDT 
counterpart. 



The complex Gaussian noise ^A(r, t) has white noise correlation function 



(^A(r,f)^;(r',f')> = 



167^2 



6(r - r')6(t - t') . 



(353) 



In a similar way one decouples quadratic in A^ terms in the action ( |339t by introducing vectorial 
Hubbard-Stratonovich field ^j(r, f) 



exp(-4rTr{crr,,,r[A^]')) = J D[^j] exp 



-Tr 



+ 2/A'3C^^j 



(354) 



where CTf , ,, - ctd + e^D5vJ?P^ + e^v5D^Ji, is the complete conductivity including both DOS and 
MT renormahzations. The resulting action is now linear in both and A^ fields, allowing us to 
define the charge p(r,f) = (l/2)55eff/(J<l>«(r, f) and current j(r, f) - (l/2)(5'5eff/5A«(r,f) densities. 
It is important to emphasize that the differentiation here is performed over the bare electromagnetic 
potentials {A, <I>), while the action ^eff in ( 13391 ) is written in terms of the gauged ones {Age, •Sac)- 
The connection between the two {O, A) {Age, <5gcl is provided by the functional 3C[<1>], which is 
implicit in ( 12331 ). A simple algebra then leads to a set of the continuity equation 5,£)(r, f)+div j(r, f) = 
0, and expression for the current density 

j(r, = J df' + ] [e^ (v + 5y^°s) E(r, t') - drg{r, t')] 



jievD 
4T 



Im{A;^'(r, t)[dr - 2/eA^3^(r, f)]A5(.(r, f)} + ^j(r, f) , 



(355) 



where E(r, f) - 5,Agc - c^r'l'ac is electric field. The current fluctuations are induced by vector 
Gaussian white noise with the correlator 

{^^(r, t) q(r', t')) = 5,, Te" (2(y + 5v^°^)D6,-,. + v5DfJ,, + v6DfJ) 5{r - r') , (356) 
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guaranteeing the validity of FDT. Equations (13521 ) and ( 1355b together with the continuity relation 
must be also supplemented by Maxwell equations for the electromagnetic potentials. 

It is instructive to rewrite TDGL equation ( 13521 ) back in the original gauge. This is achieved by 
the substitution of the gauged order parameter Aij^ = A''' exp ( - lieOC^') into ( 13521 ). This way one 
finds for the bare order parameter A'^' the following equation 



[d, - 2ied,3C''(r, 0]A"'(r, t) = 



D [dr - lie A''(r, t)f - - ^|A''(r, f)| 

n 1 



A>,0+^A(r,0, 



iA vi, (; -t- gAKi 

(357) 

where we have redefined the order parameter noise as — * exp(2/e3C'''), which, however, does 
not change its correlation function ( 13531 ). Unlike TDGL equations frequently found in the literature, 
the left hand side of ( 13571 ) contains Gor'kov-Eliashberg (GE) anomalous term diOC^'ir, t) instead of 
the scalar potential <l)'^'(r, f), see Figure [TTb. In a generic case ^^^'(r, t) is a non-local functional of 
the scalar and the longitudinal vector potentials, given by ( 12361 ). For the classical component ( 12361) 
provides 

{d, - Ddl)X'\r, t) = 0'^'(r, t) - D divA'''(r, t) . (358) 

Fields dtX'^' and (f' coincide for spatially uniform potentials, however in general they are distinct. 
The standard motivation behind writing the scalar potential <l>'^'(r, f) on the left-hand side of TDGL 
equation is the gauge invariance. Note, however, that a local gauge transformation 

j^ci ^ A''-drX, X''^X''-x, (359) 

leaves ( 13571 ) unchanged and therefore this form of TDGL equation is perfectly gauge invariant. The 
last expression in ( |359t is an immediate consequence of ( 1358b and the rules of the gauge transforma- 
tion for <l)(r, t) and A(r, f)- In the 3C gauge, specified by xii^, = 3C'^'(r, t), the anomalous GE term 
disappears from TDGL Equation ( 13571 ). and one returns back to Eq. ( 1352b . 



8.4 Applications IV: Non-uniform and fluctuating superconductivity 
8.4.1 Proximity effect 

Close to the interface with a superconductor a normal metal acquires partial superconducting proper- 
ties. At the same time the superconductor is weakened by the normal metal. This mutual influence is 
called proximity effect. The quasiclassical Usadel and kinetic equations discussed in the Section [8^ 
give full account of proximity related phenomena for superconductor-normal metal structures. One 
example of this kind is considered in this section. 

Consider a normal diffusive wire of the length L placed between two bulk superconductors, 
forming superconductor-normal metal-superconductor (SNS) junction. We are interested to study 
how the proximity to the superconductor modifies quasiparticle energy spectrum in the normal wire. 
It follows from the Usadel equation ( 13271 ) that the density of states in the wire acquires an energy 
gap €g and exhibits square-root non-analytic behavior ~ ye - eg above it, at e > eg 11321 11581 . 
To see this exphcitly we assume that the wire cross-section dimension is much smaller than the 
superconductive coherence length ^ = y/D/A. In this case the wire may be thought of as being 
quasi-one-dimensional, such that all the variations occur along the x coordinate of the wire. If there 
are no attractive interactions in the wire, A = Q, then according to the self-consistency Equation 
(1334b pair potential A(r) = within the wire -L/2 < x < L/2, and A(r) = A outside this interval. 
If in addition there is no phase difference between the two superconductors, dxX - 0, the Usadel 
equation (1327b simplifies considerably and reads as 

D dl0{x, e) + lie sinh 6{x, e) = . (360) 

At the interfaces with the superconductors, x - +L/1, this equation is supplemented by the boundary 
conditions 0{+L/l, e) - Oscsi^), where tanh Oscsi^) - A/e. It is assumed here that superconductors 
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are very large and negligibly perturbed by the wire, such that one can use coordinate-independent 
0Bcs{^) everywhere inside the superconductors. Having solved (13601 ) one finds density of states as 
v(x, e) - vRe[ cosh 6{x, e)] . 

It is convenient to perform rotation 0{x, e) = ijijl - §(x, e) such that Eq. ( I360l l becomes real and 
allows the straightforward integration 



■T—^ s K(,%.. ■ (361) 

Eth J&bcs Vsinhi?o - sinhi? 



where Erh - D/L?, §o - ??(0, e) and sinh i?bcs = VA^ - e^. Equation ( 1361b defines i?o as a 
function of energy e. Knowing i?o(e) one determines density of states in the middle of the wire as 
v(0, e) - vlm[sinht?o(e)]- 

In the limit of the long wire, ^ <K L, modifications of the density of states occur in the deep 
sub-gap limit, e <sc A. One may thus approximate j?bcs ~ and the function on the right-hand side 
of ( 1361b is essentially energy independent A'(i?o, e) ~ K{d-Q, 0). It exhibits the maximum K^^^ = 
K{^q) ~ 1.75 at X 1.5, whereas the left-hand side of ( 1361b can be larger than K^^^ for e > 
K^^Eji, - eg. Thus, for all the energies e < eg. Equation ( 1361b has only real solution for i?o and 
v(0, e) = 0, since y(0, e) oc Im[sinh)?o]. For e > eg function j?o becomes complex and gives finite 
density of states. Right above the gap, < e - <sc e^, one finds with the help of Eq. ( 13611 ) 



v(e) - 3.75-' " ^ ' ^^^^^ 

where y(e) = ^ / v(x, e)Ax is global density of states, integrated over the volume of the wire {J{ is 
the wire cross-section area, and S - l/{v^L) is its level spacing). Note that since eg ~ Erh ^ A the 
approximation d^Bcsi^ ~ Cg) ~ is well justified. 

In the opposite limit of the short wire, L <K ^, or equivalently, Ejh » A, Equation ( 13611 ) is still 
applicable. However, one must keep the full energy dependence of d-Bcsi^)- One may show that the 
energy gap is given by = A - /d,E^i^ and is only slightly smaller than the bulk gap A. This is 
natural, since the proximity effect for the short wire is expected to be strong. Immediately above the 
induced gap, the density of states again exhibits the square-root non-analyticity. The coefficient in 
front of it, however, is large, v(e) ~ 6^\ETh/^)~ sl^/eg - 1, (see 111591 ). 



8.4.2 Josephson current 

Another example which may be treated with the help of Usadel Equations ( 13261 ) and ( 1327b is the 
Josephson effect. Consider the same geometry of SNS junction, as in the previous section, assuming 
a finite phase difference between the pair potentials on the boundaries of the junction, i.e. xiL/2, e)- 
X{-L/2, e) - (p. Under this condition Josephson super-current Is{(t>) may flow across the junction. 
The aim of this section is to illustrate how Josephson phase-current relation may be obtained from 
the Usadel equations. 

For the model of step-function pair potential, A(x) = A for |x| > L/2 and A = for 1x1 < L/2, 
equations ( 1326b . ( 1327b acquire the form 

D <9,( sinh- 0d,x) = , (363a) 

D d\0 + 2ie sinh B = ^i^xxf sinh 20 . (363b) 

The latter are supplemented by the boundary conditions 0(+L/2, e) = 0bcs(^)^ while boundary con- 
dition for the ;^f-function is determined by the fixed phase across the junction mentioned above, 
X{LI2, e) - x(-L/2, e) = 0. For the short wire, L «; ^, the second term on the left-hand side 
of (1363b ) is smaller than the gradient term by ejEjii ^ 1 and thus may be neglected. Since ( 1363b ) 
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allows for the first integral sinh~ 05 xJl' = JT/i, one may eliminate d^x from Eq. ( 1363b ) and find 
L^d\6 - J'^ cosh 6/ sinh^ 6. This equation may be solved exactly 

cosh e(z, e) = cosh Oq cosh ( .^^ | , (364) 



^ sinh 00 

where Oq = 0(0, e) and z - xjL. Knowing 0{x, e), one inserts it back into the first integral of (1363b ). 

(j) = 'l^^^-^ = ^ f-i/2 "^-Z sinh^ ^(^' e)' to find 

tan((^/2) = tanh ( ."^ | . (365) 



sinh 00 \ 2 sinh 0o 

This last equation along with Eq. (|364| i taken at the NS interfaces, z - ±1/2, constitutes the sys- 
tem of the two algebraic equations for the two unknown quantities: J' and 0o. Such an alge- 
braic problem may be easily solved, resulting in J'{e,(p) - 2 sinh0oarctanh[ sinh0o tan((/i/2)] and 

sinh 00 - sinh0Bcs/ -^1 -n tan^(0/2) cosh" Obcs > where cosh0Bcs = e/ Ve^ - A^. Knowing jr(e, </>) 
one finds Josephson current with the help of 

detanh(^jlmj^(6,0), (366) 

where go is the wire conductance. Using the obtained solution for J'(e, (p) one concludes that 

ttA cos((A/2) 

Im J-(e, 0) = ' = (367) 

Ve^ - A2 cos2(0/2) 

for Acos(0/2) < e < A, and Im J"(e, 0) = otherwise. Employing ( |366l l and ( |367| i. one arrives at 
the result derived by Kulik and Omelyanchuk |1601 for the zero-temperature Josephson current of 
the short diffusive SNS junction 

(0) = cos(0/2) arctanh[ sin(0/2)] . (368) 

e 

In the original work 116011 imaginary time technique was used to derive /s((/i). This result was 
reproduced later in Ill34lll6l1l with the help of real time (energy) Usadel equation. 



8.4.3 Supression of the density of states above Tc 

Superconductor below Tc has an energy gap |A(r)| in the excitation spectrum. Superconductor 
above and far away from Tc has metallic, constant density of states. One of the manifestations of 
superconducting fluctuations in the vicinity of the transition, < T - T^. <K Tc, is the depletion 
of the density of states near the Fermi energy. Fluctuations mediated suppression of the density of 
states increases with the lowering of temperature and eventually transforms into the full gap. In this 
section we calculate the temperature dependence of this effect employing Keldysh formalism and 
compare it to the original works II1621 11631 . where Matsubara technique and analytic continuation 
procedure was used. For comprehensive discussions one may consult the recent book of Larkin and 
Varlamov ifBTl . 

Our starting point is the expression for the density of states given in terms of the Q matrix 

y(£) = ^(Tr{,T,®f,e,,))^ , 

cf. Section |6.6.2| By taking Q - K one finds v(e) - v, as it should be for a normal metal. Expanding 
Q to the quadratic order in the Cooperon fluctuations HV, see ( |336t , one finds for the density of 
states correction 

•^^(^^ ^\Y^\% «c-.a'(q)ck(-q) + -Css'ml'ei'<\)))^M ■ (369) 



81 




Figure 18: Diagram for the density of states correction, (I372l l. in the vicinity of the critical temper- 
ature Tc- Two Cooperon fields c and c*, shown by the ladders, are connected to the order parameter 
A''*''*, shown as a filled triangle, which are paired by the fluctuations propagator. 

The next step is to perform averaging over fluctuating c and c fields. For this purpose one uses ( I435I I. 
which relates Cooper modes c and c with the fluctuations of the order parameter The latter are 
governed by the following correlation functions 

(A'''(q, w)A-'(-q, -co)) ^ = ^ L^iq, , (A^'(q, w)A«'(-q, -w))^ = ^ L«(q, co) , 

(A«(q, w)A«''(-q, -oj)) ^ = ^ L\q, oj) , (A^q, coyA'^i-q, -oj)}^ = , (370) 

which follow from the time-dependent Ginzburg-Landau action ( 1341b . As a result one finds for the 
correlators of the Cooperon fields 

Q,.-.(q)cU.(-q) = - ' (371a) 

((c.,_.(q)C(-q))) = ^ "^1"^:"^'".^ • (371b) 
Inserting these into ( 1369b and summing up the two contributions, one obtains 

,,(.) = im y r ^ ^-(q^-)^^^-^^-(q^-) , (372) 

where terms proportional to F eL/^^'^\(jj) in the averages {{cc*)) and {{cc*}) drop out from ( 13721 ) upon 
cl) integration, as being integrals of purely advanced and retarded functions, respectively. Equation 
( 13721 ) allows convenient diagrammatic representation shown in Figure [18] Using now fluctuation 
propagator in the form of ( 1444b and approximating bosonic distribution function as B^j ~ 2Tc/cl>, 
since the relevant frequencies oj ~ T - T^, the density of states correction ( 13721 ) reduces to 



^v(.) = -l^Re y r ^-^ , , (373) 

^J-<>« [{Dq^ + T^l)^+co^][Dq^-2is + iajY 



where T^l ^^{T -Tc)/n. 

The further analysis of this expression depends on the effective dimensionality of the system. 
We focus on quasi-two-dimensional case: a metal film with the thickness b which is much smaller 
then superconducting coherence length b <k ^{T) - V^tgl- Replace the momentum summation by 
the integration Yjq ^ If introduce the dimensionless parameters x — Dq^/Tc and y - oj/Tc, 
integrates over y using residue theorem and finds 

'^^'^ ' J/(-GL), m = Re r '\ ... , (374) 



16 \ r - r, / Jo (1 + Ji:)(l +2x- 2/z)2 ' 
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where Gi = ti/vDb is the Ginzburg number. For small deviations from the Fermi energy, £Tgl ^ 1, 
the DOS suppression scales as 6v{Q) oc -{T/Tc - 1)"^, while at larger energies stgl 1 DOS 
approaches its its normal value as 6v(s) oc -(Tc/e)^ In(eTGL)- Note also that Jde5y(e) = 0, which 
is expected, since the fluctuations only redistribute states around the Fermi energy. 



8.4.4 Fluctuation corrections to the conductivity 

Superconductive fluctuations above Tc modify not only the density of states, but also transport prop- 
erties. In the case of conductivity, there are three types of the corrections, density of states (DOS) 
^CToos, Aslamazov-Larkin (AL) Sctal and Maki-Thompson (MT) So-mt terms, see. II155II1561I157I 
11621 . Although we have already partially discussed this topic in Section [831 the goal of this section 
is to show explicitly how all of them are obtained within Keldysh cr-model approach. 

According to the definition given by ( 1201b . to find conductivity one needs to know partition 
function 2[A'^', A''] to the quadratic order in vector potential. Using ( 1312b one finds PI 

Z[A^U^]^ J D[e,A] l + ^Tr{SAeSAe}-^^^(Tr{5re[HA,e]))' exp {iS AQ, A]) , 

(375) 

where diamagnetic contribution Tr{BABA) was omitted. As it was demonstrated in the Section l6.6.3l 
by taking Q - A and using ( 12011 ) one finds Drude conductivity ctd- To capture superconductive 
corrections dcr to normal metal conductivity cr^, one has to expand g-matrix in fluctuations 'W to 
the leading (quadratic) order and analyze all possible contributions. 

From the first trace on the right-hand side of ( 13751 ) by taking one of the Q matrices to be A, 
while expanding the other one to "W" order, one finds 

Zdos[A'-', A^'] = ^ ((Tr{A,,,,(,T, ® f,)A,,,3(,T, ® f,ms,s,^s,s, ]))^^ , (376) 
where the current vertex matrix is 

It will be shown momentarily, that Zoos defines density of states type contribution to the conduc- 
tivity in the vicinity of the critical temperature. Indeed, one substitutes ( 13761 ) into ( |201b . carries 
differentiation over the vector potentials, takes the dc limit Q ^ and evaluates matrix traces. As a 
result, one fids 

^troos = X JJ ^'^^-^^^ ((c...4(q)<..(-q) + c....(q)c;,,(-q)»^_^ . (378) 

As the next step, one uses ( 14351 ) and performs A averaging with the help of correlation functions 
Eq. ( 1370b . Changing integration variables S2 s and £4 — > e - w, correction Sctoos becomes 

e^D r+°° L^(q, w) + Fs_„L*(q, w) 

So-Dos = Im ) \ dedcj d^Fe ~ . (379) 

27T Y {Dq^ - 2ie + icof 

By comparing this expression to ( 1372b one concludes that Sctdos oc J dedEFgdvie), which estab- 
lishes connection between Scruos and density of states suppression 6v{s), see also Figure [T9b for 
diagrammatic representation. In order to extract the most divergent part of 60-^05, in powers of 



"since Coulomb interactions do not lead to a singular temperature dependence for kinetic coefficients in the vicinity of 
Tc, we set <^x = and suppress subscript 3C throughout this section. 
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Figure 19: Diagrams for superconductive fluctuation corrections to the conductivity in a vicinity of 
Tc- a) density of states term; b) Maki-Thompson correction; c) Aslamazov-Larkin correction. 

the deviation T - Tc, one only needs to keep in ( 1379b Keldysh propagator. The Fg^c^L'^ term gives 
parametrically smaller contribution. Using (I4441 i one finds 



Sctdos = ^ Re > dedw ; 



(380) 



After remaining frequency and momentum integrations, for the quasi-two-dimensional case, one 
finds 



o-D \T - Tc 

This correction is negative as expected, which stems from the depletion of the density of states by 
fluctuations, and has relatively weak temperature dependence. It is worth emphasizing that dcroos 
can be extracted from the effective time dependent Ginzburg-Landau theory, which was discussed in 
the Section lOl Indeed, one can show that dcroos = e~Z)((5v[?p^)A, where 5vJ?p^ is taken from ( 1350b . 
reproduces ( 1381b . 

Let us return back to ( 1375b and look for different possible contributions. Focusing again on the 
first trace on the right-hand side of ( 13751 ), one may expand now each of the Q matrices to the first 
order in fluctuations "W. This way one identifies 

Zmt[A''', A«] = ^ ((Tr{A,„,((r, ® t^nV,,,, A,,,,((t, ® f,)^,,„ }))^^ , (382) 

which leads to Maki-Thompson correction to the conductivity. After differentiation of -3mt[A'^', A^] 
over the vector potential, and evaluation of the traces, in the dc limit, one finds 



As compared to dcroos in ( 1378b Scr^j consists of products of mixed retarded c and advanced c 
Cooperons, while dcroos contains Cooperon fields of the same causality. Using ( 13701 ) and ( 1435b one 
carries averaging in ( 1383b over A fluctuations, then changes integration variables in the same way as 
in ( 13791 ) and arrives at 

So-MT = > ^ dedw deF^ — ^ ^ .2 — • (384) 
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The corresponding diagram is shown in Figure [T9b . With the same accuracy as earlier, approx- 
imating a! 2Tc/cj, neglecting F^-^ and using ( I444l i for the fluctuations propagator, the latter 
expression for So-mt reduces to 



50-MT = ^^^^-5-^ y rr dsdw- — " ^ ^ . (385) 



Finally, after the remaining integrations for quasi-two-dimensional case, one finds 



In — , (386) 



o-D 8 \r - r,./\l -tgl/t^/ \tgl 

where infrared divergency in momentum integral was cut off' by a dephasing rate Dq^^^^ = r^'. 
This divergency is a well-known feature of the Maki-Thompson diagram. It can be regularized by 
some phase-braking mechanism in the Cooper channel. For example, if magnetic impurities are 
present in the system, then the role of is played by the spin-flip time. In contrast to dcruos Maki- 
Thompson correction ( 1386b is positive and has much stronger (power law) temperature dependence. 
Interestingly, that Sctmt follows from the effective Ginzburg-Landau theory as well. Indeed, defining 
Sctmt - e-v{SD^Jj,)i^, employing ( 1347b and performing averaging over A, one recovers (1386b . 

There is yet another correction to conductivity, called Aslamazov-Larkin contribution. It is 
obtained from the second trace on the right-hand side of (1375b . Indeed, expanding each Q matrix to 
the linear order in 'W, one finds 



ZalW,^'] = -^^^ (((Tr{A,„,(<T,®r,)nV,,,,W,„,}f))^^ . 



(387) 



It is convenient to introduce two vertices, which follows from ( 1387b after differentiation over the 
vector potential 

Vim = --i_Tr{A,,,,((r,®f,)^,,,35rnV,3,,) 

= tr(c£,e3(r)5r<3^,_^^(r) + c;,e3(r)5rCe,s,+n(r) - (c c)} , (388a) 
S 



= -tr{F,3 {r)drcl,^_ (r) + (r)5rQ„, (r) + (c -> c))} . (388b) 

Notice that for it is sufficient to take external frequency to be zero right away, Q = 0, while 
for V^^L it is important to keep finite Q. and take the dc limit, — » 0, only after averaging. 
Performing averaging over Cooperons, one uses ( 1435b . In the case of Vl['VV'], for the product of 
two Cooper fields it is sufficient to retain only contributions with classical components of the order 
parameter, Vl[aV] oc Tr{F[c5rC* + c*5rc]} oc A'^'^rA*"^' - A^'^rA"^'. In contrast, for the Y'^^YW] 
vertex, it is crucial to keep at least one quantum component of the order parameter Isfl, since the 
corresponding contribution with two classical components vanishes owing to causality structure. As 
a result, the leading contribution is ¥^[^1 oc Tx[cdrC* + c*drc] oc A'''5rA*'' + A''dr^*''' - c.c. The 
remaining A averaging of the product (Vl[ny]Vl[nV])^ is done with the help of ( 1370b . Passing 
to the momentum representation and collecting all the factors, Aslamazov-Larkin type correction to 
conductivity in the dc limit takes the form 



2 V J J g I; [coth ^] [lmL\q, ^)]^ . (389) 
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The corresponding diagram is shown in Figure [T9t. Since only Dq^ ~ cj ~ <s; Tc are relevant, 
one may approximate 5„[coth w/2r] ^ -ITdciP' and use ImL'*(q, uj) = —{%iT c(o I n){{Dq- + t'^^^ + 
aP'Y^ to obtain 

> I at,) 



^(TAL = > J dw J . (390) 

Performing remaining integrations, one finds for the quasi-two-dimensional film 



o-D i6\r-r. 

At the level of eff'ective time-dependent GL functional, Aslamazov-Larkin conductivity cor- 
rection dcr AL appears from the ^sc part of the action ( I344l i. The easiest way to see this is to use 
current 

jsc~Im[A«UA^'], 

which follows from ^sc- along with the fluctuation-dissipation relation Scr^i^ cc (j^Q ■ jsc)^ 
Zqoj Dq^\L'^(q, w)|l The latter reproduces (|390l l. 

The technique which was employed within this section allows to reproduce all the results for 
fluctuations-induced conductivity, known from conventional Matsubara diagrammatic approach. 
The simplification here is that no analytical continuation was needed. Although it is not so com- 
plicated for the problem at hand, in some cases avoiding the analytical continuation may be an 
advantage. 

8.4.5 l\inneling conductance above Tc 

Consider voltage-biased superconductor-normal metal tunnel junction, where the superconductor 
is assumed to be at the temperature just above the transition T^, i.e. in the fluctuating regime. It is 
natural to expect that depletion in the density of states, mediated by fluctuations (see Section [8.4.3l l. 
modifies current- voltage characteristics of the junction II1641 11651 11661 . This effect can be studied 
within cr-model, using tunneling part of the action St[Ql, Qr]- 

One starts from (|179t and performs gauge transformation Qa — > e\p{-iEi>a)Qa exp(E<l)a), for 
a - L,R, where ia(t) - f Va(t)dt - [<l>^J{t)d-o + <i>t{t)d-x] ® tq, and Of - Of - eVt, which moves 
an applied voltage V from the Keldysh blocks of the Q matrices, to the tunneling part of the action 

iSAQi, Qr] = ;r-Tr{eie-'^*e«e'^*} , (392) 

here <[) = - 0«, and <E>^(f) serves as the generating field. Indeed, since the phase <t> is quantum 
canonical conjugate to the number of particles .^V = id/d^ the tunneling current is obtained by differ- 
entiating the partition function -37-['l'] - ^'s^P (iS t[Ql, Qr]) with respect to the quantum component 
of the phase 

IHO^iei'-^] . (393) 



Applying this definition to ( 13921 ). using 
5exp(+/£6) 



= +i6{t - t') (&, ® T.) exp [ + ieVtE] (394) 

09=0 



and taking Ql-Qr- A, one finds Ohm's law It - grV, as it should be, for the tunneling junction 
in the normal state. One may account now for the fluctuation effects by expanding one of the Q 
matrices in ( |392t over Cooper modes 'W. This leads to the correction of the form 

TTgy \— 1 rr dede' , ,, 

ShiV) - 2j Jj - ''^-'^^ {{css- (qK^i-q) + -Cee-{q)-c:-si-(l)))^,A . (395) 
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which is physically expected result. Indeed, from the combination of the Cooper modes in ( |395l l 
one recognizes density of states correction 6v{s), see ( 13691 ). The latter is convoluted in ( |395l l with 
the difference of Fermi functions, leading to the correction to the tunneling current of the form 
6It{V) ~ Jds[Fs+eV - FB-ev]SvL{s)vR. Using previous result for 6v(e) from ( 1373b . bringing it 
into ( 13951 ) and transforming to the dimensionless units x = Dq^/T, y - to/T, z = s/2T one finds for 
the tunneling differential conductance correction dgriV) = d6lT(V)/dV the following expression: 



6gT(V) 

gT 



dt/dz 



1 



1 



xRe 



cosh^(z + u) cosh"(z - u) 
1 



(396) 
(397) 



where u - eV/2T and we assumed quasi-two-dimensional geometry. Remaining integrations can 
be done in the closed form, resulting in II 16411 



6gT{V) Gi , 
= ^ In 

gT 



Tc 



T-Tr 



1 



ieV 



(398) 



where if/^^\x) is the second order derivative of the digamma function if/{x). Notice, that although 
having direct relation to the density of states suppression Sv{s), the differential conductance correc- 
tion SgT exhibits much weaker temperature dependence. The sharp suppression in the density of 
states (5v(0) cc {T - T^y^ translates only into the moderate logarithmic in temperature correction 
6gT oc In(rcTGL)- Another interesting feature is that suppression of the Sv{s) occurs at the energies 
s ~ Tq]^ ~ T - Tc, while corresponding suppression of the differential conductance happens at volt- 
ages V ~ Tc, and not at V ~ T - Tc- Finally one should mention, that more singular in {T - Tc) 
MT and AL corrections appear only in the fourth order in the tunneling matrix elements, while the 
discussed DOS effect is linear in gj- (i.e. it is of the second order in the tunneling matrix elements). 



8.4.6 Current noise in fluctuating regime 

Apart from the density of states related effects, there are interesting consequences of superconducting 
fluctuations on the current noise of the tunneling junction II 15 II 11671 11681 1 1691 1 1701 . Assume now 
that both sides of the junction are made from identical superconductors that are kept right above Tc. 
While there is no average Josephson current in this case, the noise power turns out to be sensitive to 
the Jesephson frequency, a>j = 2eV/Ji, and exhibits sharp peak at o) - ujj. The height and shape of 
this peak have a singular temperature dependence near Tc, which makes its experimental detection 
possible. To show this we establish an expression for the fluctuating part of the tunneling current 
dlrit) in terms of the product of fluctuating order parameters A/,(s)(r, f) residing on the different sides 
of the junction, namely Slrit) oc J dr[AR{r,t)A*^(r,t)e\p(-ia)jt) - cc.]. Since (Al(r)) = above 
Tc, it is clear that (Slrit)} - 0. However, the average square of the current {6lT(t)6lT{t')) is not 
vanishing and its Fourier transform displays a peak at the Josephson frequency. In what follows we 
calculate its temperature dependence. 

One starts from the definition of the current-current correlation function 

In the normal state Ql - Qr - ^ and the noise power of the tunneling junction, as it follows 
from (13991 ), is given by the Schottky formula >S7-((:ij) = IgjT v-t coth v±, where v-t - {eV + a))/2T. 
To account for the superconductive fluctuations on both sides of the junction one has to expand 
each of the Q matrices in ( 1392b to the leading (linear) order in Cooper modes. This gives for the 
fluctuation correction to the current 

Sl-rit) = ^ s^^^'^ [e'^^€{{^, ® tJOVl'Z/"' e"'^*'l/(o-; ® fJoVs'l/"' } . (400) 
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To proceed further, one simplifies ( 1400b . exploring separation of time scales between electronic and 
order parameter degrees of freedom. Indeed, one should notice that, as follows from (I444l i. the 



relevant energies and momenta for the order-parameter variations are Dq^ ~ w ~ r^J^, while the 
relevant fermionic energies entering the Cooperons are e ~ e' ~ T » t^}^. As a result, non-local 
relations between Cooper modes (1337b and the order parameter, see Eqs. ( 14351 ). may be approximated 
as 

A«,(r) = A:;'(r,^Jf++Ar'(r,^jf_, a^L,R, (401) 

where Q(f) is the step function. Physically (1401b implies that Cooperon is short-ranged, having 
characteristic length scale = '^D/T^, as compared to the long-ranged fluctuations of the order 
parameter, which propagates to the distances of the order of ^gl - VD^gl ^o- Thus, relations 
( 14351 ) are effectively local, which considerably simplifies the further analysis. Equations ( 14011 ) allow 
us to trace Keldysh subspace in ( 14001 ) explicitly to arrive at 

dlrit) - -^Tr {e{t2 - h)Fr,-,0{t - f2)A^^t-A«,_e'^^<'+'^'^=}^ , (402) 

where we have used ( 1394b and wrote trace in the real space-time representation (note that Tr{. . .} 
here does not imply time f integration). Changing integration variables f i - t - fi and tj, - t - rj, and 
rescaling rj, fi in the units of temperature Trj — » r],Tfi yu, one finds for Eq. ( 1402b an equivalent 
representation, 

W) = -^ [rd,d, ^^^^TrlAf, ,^Af , (403) 

eT J J sinh(;ryu) [ '•' t ' t^' t )n 

where we used equilibrium fermionic distribution function in the time domain F, = -iT/ sinh(;7rrf). 
The most significant contribution to the above integrals comes from ?; ~ // < 1. At this range ratios 
{r],fj.}/T change on the scale of inverse temperature, while as we already discussed, order-parameter 
variations are set by f ~ tgl s> l/T. Thus, performing 77 and// integrations one may neglect {77, yu)/r 
dependence of the order parameters and the exponent. As a result one finds 



(404) 



where the spatial integration runs over the junction area ^ and a>j = 2eV/fi. Finally one is ready 
to calculate corresponding contribution to the current noise. One substitutes two currents in the 
form of ( 1404b into (1399b and pairs fluctuating order parameters using the correlation function, which 
follows from (1370b . 

(A:;'(r, OAf(r', 0)^ = ^S^hL^ir -v',t- t') . 
As a result, superconducting fluctuation correction to the noise power is given by 

<5^r(a>) = -^(^J df[LV,0]' exp(-/a;±0, (405) 

where a;± = co + loj. Performing the remaining integrations one finds first Keldysh component of the 
fluctuation propagator in the mixed momentum/time representation L*- (q, = / (fl- w)e"™'dw/27r, 
which is 

^''(q, f) = - ' = (^GL?)' + 1 . (406) 

1 I r 
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Then insert L*^(r, t) - J L'^{q, t)e"^'dq^ /A-n into (1405b . introduces dimensionless time r = t/rcL, and 
changes from q to Kg integration dq^ = dxg/^^^, which gives ahogether II 1701 

where the spectral function is given by 

dT — exp(-2;^|T| - izT) = - In ^/l+J/4. (408) 

CO Jl ^ 2 

The noise power correction 6St{cl>) is peaked at the Josephson frequency ai - and has strong 
temperature dependence, which makes its experimental detection possible in a vicinity of the super- 
conducting transition. 
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A Gaussian integrals for bosons and fermions 



For any complex N x N matrix Ajj, where i,j - 1, . . .N, such that all its eigenvalues, /I,, have a 
positive real part, ReA, > 0, the following statement holds 



Z[J] 



d(Rez;)d(Imzj) 



exp 



Det(A) 



(409) 

where Jj is an arbitrary complex vector To prove it, one may start from a Hermitian matrix, that 
is diagonalized by a unitary transformation: A = U^AU, where A = diag{/lj). The identity is then 
easily proven by a change of variables (with unit Jacobian) to ly, = UijZj- Finally, one notices 
that the right-hand side of ( 1409b is an analytic function of both ReA,j and ImA,j. Therefore, one 
may continue them analytically to the complex plane to reach an arbitrary complex matrix A,y. The 
identity (14091 ) is thus valid as long as the integral is well defined, that is all the eigenvalues of A/j 
have a positive real part. 

The Wick theorem deals with the average value of a string Zo, . . . Zfl,Z/7, . . . Zb^ weighted with the 
factor exp( - Y^ijZiAijZj). The theorem states that this average is given by the sum of all possible 
products of pair-wise averages. For example. 



{ZaZb) = 



1 6^Z[J] 



{ZaZbZcZd) = 



Z[0] 6Ja6Jb 

1 



ab ■■ 



(410) 



7=0 

S'^Z[J] 



Z[0] 6Ja6JbSJcSJd 



^ac^bd ^r,lA,} 



ad be ' 



and so on. 

The Gaussian identity for integration over real variables has the form 



Z[7] 



X+oo A' 
n 
°° ./=i 



exp 



■ ^ XiAjjXj + 2 ^ XjJj 



ew{Y,f.MA-%Jj) 



VDet(A) 



(411) 



where A is a symmetric complex matrix with all its eigenvalues having a positive real part. The proof 
is similar to the proof in the case of complex variables: one starts from a real symmetric matrix, that 
may be diagonalized by an orthogonal transformation. The identity ( 141 lb is then easily proved by 
the change of variables. Finally, one may analytically continue the right-hand side (as long as the 
integral is well defined) from a real symmetric matrix A,y, to a complex symmetric matrix. 

For an integration over two sets of independent Grassmann variables, and ^j, where j = 
1,2,. . .,N, the Gaussian identity is valid for any invertible complex matrix A 



rr ^ 



A^jA^j exp 



- Det(A) exp 



-Yj^iAij^j + J][lxj+Xj€j] 
U j 

(412) 

Here Xj and Xj are two additional mutually independent (and independent from and ^j) sets of 
Grassmann numbers. The proof may be obtained by, e.g., brute force expansion of the exponential 
factors, while noticing that only terms that are linear in all 2N variables and are non-zero. 
The Wick theorem is formulated in the same manner as for the bosonic case, with the exception that 
every combination is multiplied by the parity of the corresponding permutation. For example, the 
first term on the right-hand side of the second expression of (14 10b comes with the minus sign. 



B Single particle quantum mechanics 

The simplest many-body system of a single bosonic state (considered in Section |2]i is, of course, 
equivalent to a single-particle harmonic oscillator. To make this connection explicit, consider the 
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Keldysh contour action ( fTSl l with the correlator ( fTTb written in terms of the complex field 0(f). The 
latter may be parameterized by its real and imaginary parts as 



0(0 = 



1 



■.{p(t)-iaJoq(t)), 0(0 = 



In terms of the real fields p{t) and q{t) the action ( fTSl l takes the form 



(/?(0 + /wo^(0)- 



S [p, q] 



£dt pq-^{p- + (olq^) 



(413) 



(414) 



where the full time derivatives of p', q^ and pq were omitted, since they contribute only to the 
boundary terms, not written explicitly in the continuum notation (they have to be kept for the proper 
regularization). Equation ( 1414b is nothing but the action of the quantum harmonic oscillator in the 
Hamiltonian form. One may perform the Gaussian integration over the p{t) field to obtain 



S[q] = ^^df[^2_^2^2j 



(415) 



This is the Feynman Lagrangian action of the harmonic oscillator, written on the Keldysh contour. 
It may be generalized for an arbitrary single-particle potential U(q) 



S[qit)] = j^dt 



U{q{t)) 



(416) 



One may split the q{t) field into two components, q+(t) and q-{t), residing on the forward and 
backward branches of the contour, and then perform the Keldysh rotation: q± = q'^' + q''. In terms of 
these fields the action takes the form 



SWKq"] 



f 

%J-cx 



dt 



-2q' 



j'> cl 

, d-q" 
dfi 



U{q^' + q^) + U{q^'-q^) 



(417) 



where integration by parts was performed in the term q'^q'^'. This is the Keldysh form of the Feynman 
path integral. The omitted boundary terms provide a convergence factor of the form ~ iQ{q'')-. 

If the fluctuations of the quantum component q'^{t) are regarded as small, one may expand the 
potential to the first order and find for the action 



S{q' 



', q''] ^ J dt -2 q'l ^ 



df- 



dq'^' 



+ iOiq'^f + 0[(q'')'] 



(418) 



In this limit the integration over the quantum component, q"^, may be explicitly performed, leading 
to a functional (5-function of the expression in the round brackets. This ^-function enforces the 
classical Newtonian dynamics of q'^' 



^q'' _ dU{q'') 
dq'' 



dfi 



(419) 



For this reason the symmetric (over forward and backward branches) part of the Keldysh field is 
called the classical component. The quantum-mechanical information is contained in the higher- 
order terms in q'^', omitted in ( 1418b . Note that for the harmonic oscillator potential the terms denoted 
as Oliq^Y} are absent identically. The quantum (semiclassical) information resides, thus, in the 
convergence term, iO{q'')-, as well as in the retarded regularization of the d^/(df^) operator in ( 1418b . 

One may generalize the single-particle quantum mechanics onto a chain (or lattice) of harmon- 
ically coupled particles by assigning an index r to particle coordinates, qr{t), and adding the spring 
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potential energy, ^{qr+i(t) - qAO)^- Changing to spatially continuum notation, i 
finds for the Keldysh action of the real (e.g. phonon) field 



dt 



1-7 . 



(r, f) = qr(t), one 



(420) 



where the constant has the meaning of the sound velocity. Finally, splitting the field into (0+, (p.) 
components and performing the Keldysh transformation, (f>± = <p'^' + 0*, and integrating by parts, one 
obtains 



dr dt[2f{'d^ + u^^dl 

%J — OO 



'') + U{4>" -4>')]- (421) 



According to the general structure of the Keldysh theory the differential operator ( - 5^ + d^), 
should be understood as the retarded one. This means it is a lower-triangular matrix in the time 
domain. Actually, one may symmetrize the action by performing the integration by parts, and write 
it as 4>''{ - c^? + fj dl)'^<p'^' + <p'^'{ - + d^Y^cffl, with the advanced regularization in the second term. 



C Gradient expansion of the cr-model 

This Appendix serves as the complementary material for Section 16.21 Its purpose is to provide 
technical details hidden behind the transition from ( 11731 ) to ( 11741 ). For the gradient expansion of the 
logarithm in ( |173t one uses Q matrix in the form of ( 11671 ) and finds in analogy with ( |168t : 

iS{Q,K y] = Trhi [i + igndtk'^ + ignypdA^^ + gnv-R-^ + g-kyp^'R'^] ■ (422) 

Expanding this expression to the linear order in Q'Rd/Rr^ and quadratic in Q'tiypd/Rr^ , one repro- 
duces ( 11711) for S {Q\, which leads eventually to ( 1172b . To the linear order in V and A one finds 
from(|422ll 

iS 1 [Q, A, y] = TY[g^V^-^ I - /Trl^C^VfSr'R"' )^(:Rvf A-R-' )) . (423) 

In view of Y,p ^(P> - —iJtvK^, which follows from the saddle point Equation ( |162t . for the 
first term on the right-hand side of ( 142 3 1 one finds, using cyclic property of trace Trf^^V''^"'} = 
-ijivTr['Rr^ MiV]- -ijivTx[VQ]. As to the second term on the right-hand side of ( 14231 ). retaining 
retarded-advanced products of the Green functions 2/, ^'*(P> e)vF^'^(p, e)\F - 2nvD, one finds 

where K o 5r'^ ' = -d^'R o ■^ ^ was used. All together it gives for ( 14231 ) 

iS 1 [Q, A, V] = -invTi{VQ] + i7ivDi:x[kQdrQ] . (424) 
To the second order in V and A one finds 

iSiiaA, V] = -^Tr{^y^y} - iTr{^(!Rvf A5R-i)^(^Vf A^-')} . (425) 

Note that in the term ~ we took = = 1. This is because ~ contribution represents 
essentially static compressibility of the electron gas which is determined by the entire energy band, 
while and ' matrices are different from unit matrix only in the narrow energy strip around the 
Fermi energy. Thus, for the first term on the right-hand side of (1425b one can write Ti{gvgV} = 
Tr{y"Y«^y^}, where 

t"^ = Z r r^^'f^^P' f+)^"'^(P' ' e^^e±cj/2, (426) 
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and trace spans only over the Keldysh matrix structure. Using (11691 ) for the matrix Green function, 
and retaining only retarded-retarded and advanced-advanced products one finds 

p 

(421) 

This result is derived noticing that [^**'*'(p, e)]" = -5e^**'**(p, e), and integrating by parts 

p p 

using 2^ (^*(p, e) -^'^(p, e)) = -2mv and assuming that F^^±oo ±1- An additional contribution 
to T"'', originating from the retarded-advanced products of Green's functions, although non-zero, 
contains an extra small factor wTei <si 1, and thus neglected. 

For the second term on the right hand side of ( 14251 ) one finds Tr{^('^V/7A'^"')^('^V/7A!R"')} = 
7TvDTr{[i+A]'kA'k-^[i-A]'kA'k-^] = 7tvDTi{A- -AQAQ], which finally gives for the S2[Q,A,V] 
part of the action 

iS2[Q,A, V] = ~TT{Vcr,V] + '^Tv{AQAQ - A'] . (429) 

Combining now (1172b together with S \[Q,A,V], and 5'2[2, A, V], and taking into account that 
TT{{drQ)^ -4-iAQdrQ-2{AQAQ- A^)} - TT{{drQ)^], where covariant derivative is defined accord- 
ing to (1175b . one finds the full action in the form of (1174b . 



D Expansion over superconducting fluctuations 

In this section we provide details of the Gaussian integration over the Cooper modes performed 
in ( 1338b . Throughout this section we suppress subscript-3C in Qjc and Ajc for brevity. As a first step 
one expands ( 13121 ) in fluctuations 'W around the metallic saddle point Q - A: S [Q, A] 5 ['W, A]. 
To this end, we take 'W from ( 1337b and substitute it into (1312b ). For the spatial gradient part of 
the action So- one finds in quadratic order Tr{(cJr2)^} - Trl'Wse'd'^'^E'e]- Tracing the latter over 
Keldysh®Nambu space gives 

D TiiidrQf] = 2 Jj ^ Dq^[c*^,{q)c,,A-q) + c:,,(q)c,,,(-q)] , (430) 

where we kept only Cooper modes c and c, while omitting the difiiison modes d and d, since expan- 
sion for the latter was already given in ( |187b . For the time derivative term in the action So- one finds 
Tr{HcJ,2) = -^Tr{e((T, ® fo)'WsE''^e'e], where we took 3, — > -is in the energy space. The latter, 
after evaluation of the trace reduces to 

Tr{Ed,Q} = ^ JJ ^ + e')[cL'(q)c.'.(-q) - C'(q)c.'.(-q)] ■ (431) 

To the leading order in "W the coupling term between Cooper modes and the order parameter. A, 
reads as TrfAQ} = Ti\€(eAe-E'€(-}id-:, ® t,)^^,^} + O(A'W'), where €( is given by Eq. ( |337] |. 
Evaluating traces, one finds 

TrfAei = Z JJ ^ [A^,,(q)c:,,(-q) - AUqKs(-<D ' c-c] , (432) 

where the following form factors were introduced 

A^,,(q) = A^'(q, s - s') + F.A^q, s - s'), A^,,(q) = A'^'(q, s - s') - F,,A?(q, s - s') . (433) 
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It is important to emphasize, that the diffusion modes {d,d} couple to A only starting from the 
quadratic order in 'W. These terms produce non-local and non-linear interaction vertices be- 
tween the order parameter components and will not be considered here, see Ml 5411 for more de- 
tails. Combining now ( I430l )-( l432b . one finds for the quadratic part of the action ^(^[nV, A] - 
S'^VW,A]+SlVW,Al where 



A] = -— tr{c:,,(q)[D^^-/(e+e')]Q'.(-q)+2/A';,,(q)c:,,(-q)-2/A:^,(q)c,,,(-q)} 



(434a) 



iSi[^,A] = -ytr{c:,,(q)[Z)^2+/(e+eO]c.'.(-q)-2/A^^,(q)c:,,(-q)+2/A;^,(q)c,,,(-q)}, (434b) 

and traces stand for energy and momentum integrations tr - Yjq Jf At this stage, one is 

prepared to perform Gaussian integration over the Cooper modes c and c. Quadratic forms in ( 14341 ) 
are extremized by 



:'(q) 



-2/A^,,(q) 
Dq^ - i{s + e')' 



:'(q) 



2'-AL,(q) 
Dq^ + i{e -He') 



(435) 



Similar equations for the conjugated fields, are obtained from (1435b by replacing A — > A* and 
flipping an overall sign. The Gaussian integral JD[^] exp(/5'o-['VV', A]) - exp(/5'o-[A]), where 
5(j-[A] is calculated on the extremum ( 1435b : 



(436) 



Dq^ - lie 



where A^''''^ = A'''^*\+q, +w) and e± - e ± (jj/2. We have also introduced new integration variables 
u)-E-e',e-{E + e')/2 and employed the fact that is an odd function to change variables as 
e — > -e in the contribution coming from c fields. The contribution to /5'o-[A] with the two classical 
components of the order parameter ~ A'^jA*^' vanishes identically after the e-integration as being 
an integral of the purely retarded function. This is nothing, but manifestation of the normalization 
condition for the Keldysh-type action (see Section 1273] for discussions). Adding to /^(^[A] zero in 
the form -Anv\i\A\A*l' l[Dq^ - 2;e]}, which vanishes after e integration by causality, and combin- 
ing ( 1436b with 5a from Eq. ( 1312b ). one finds for 5gl[A] = 5,r[A] + 5 a[A] the foflowing result 

5gl[A] - 2v;^ J ^ [A"'i«'A;' + A-L-^'Al + A^'B.IL^' - L-^']Al] , (437) 



where superconducting fluctuations propagator is given by the integral 

Fel^wjl 



Dq^ - lie 



(438) 



This expression for L(q, w) can be reduced to the more famiUar form. Indeed, adding and subtracting 
right-hand side of ( 1438b taken at zero frequency and momentum one writes 



L«(q,^) = --.J de--,l 



ds 



Dq^ 



■ lis Is 



(439) 



where the second term on the right-hand side is the logarithmically divergent integral which is to be 
cut in the standard way by the Debye frequency ojd- Introducing dimensionless variable x - s/lT, 
and performing the integration in the last term on the right-hand side of ( 1439b by parts with the help 
of the identity dxln(x)sech~(x) = -In^, where y - with C - 0.577 is the Euler constant, 
and using the definition of the superconductive transition temperature Tc - (2ywo/7r) exp(-l//ly), 
one finds for ( |439b 



L„\q, cij) - \n — - - dx 



tanh(x) tanh(jc) 



Dq--iat 
AT 



(440) 
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With the help of the expansion 



tanh(x) = V ^ , x„ = 7r(n + 1/2), 



(441) 



«=o 



one may perform the jc-integration explicitly interchanging the order of summation and integration 

xdx in 



dx _ n_ x_ 

J -co + X„ J_co + x^] 1^ 

Recalling now the definition of the digamma function 



Dq- -ico 
AT 



Dq--iiA) 



+ X„ 



(442) 



«A(x) = -C - 

;i=0 

one transforms ( |438l l to the final result 



1 



1 



n + X n + \ 



T. 

Ls'(q,w) = lny - (A 



Dq^ - ioj 1 



47rr 



—{Dq 



(443) 



(444) 



where t^^\^ = ^{T -Tc)/n:. Since according to the last expression Df^- ~ co ~ t^}^ <s; T, the expansion 
of the digamma function is justified. 

As a result, the time dependent Ginzburg-Landau part of the effective action (I339l l is obtained 
(compare ( 1437b and ( 1444b with (1341b ). The non-linear contribution ~ |Ap in ( 1343b can be restored 
once ~ Trj^^'A) is kept in the expansion of TrIQA) part of the action. Furthermore, for Dq^ 
-Ddl in (1444b . one actually has D{dr - 2/eA'^)^, once the vector potential is kept explicitly in the 
action. 

Let us comment now on the origin of the other terms in the effective action ( 1339b . The super- 
current part of the action S sc emerges from the Ti{drQx['B.Ajc, Qjc]] upon second order expansion 
over the Cooper modes, namely 



5sc[A,A,0] = — Tr{4,(r)Kr Q'r(r) + c;,'(r)3^fr c„(r)}. 



where 



2eD 



-divA'^(r, f) + A^(r, t)[dr - lieA'^ir, t)] 



Deriving 3\r„r one uses an approximation for the equilibrium Fermi function 

iT i»i/T i , 
smh(7rrf) 2T 



(445) 



(446) 



(447) 



which is applicable for slowly varying external fields. Performing integration over the Cooper modes 
one substitutes (1435b into ( 1445b . Noticing that in the real-space representation ( 14351 ) reads as 



c„,(r) = -ie(t-t')A: 
cAr) ^ ie(t - t')K 



t + t' 
t + f 



t + 1 
t + 1' 



(448a) 
(448b) 



Xit) 



de / e \ e-'^' 2 

— tanh — = -arctanh exp(-7rr|r|)) , (448c) 

J_oo 2n \2T / e + lO n 

and keeping contributions only with the classical components of fluctuating order parameter, since 
X^'' is already linear in quantum field A^, one can perform t' integration in (14451 ) explicitly and 
recover S sc in the form given by ( 1344b . 
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The Maki-Thompson part of the effective action Smt emerges from Tr{([EA3c, Gael)"} when 
each Qoc matrix is expanded to the first order in fluctuations 'W: 

5mt[A, a, 0)] = -— Tr{4,(r)?^rc„(r) + cUr)K''c,',(r)] , (449) 



where 



A?c(^0 + ^M5^(r,f) 



Al.(r,f'), (450) 



and we again used (1447b . With the help of (I448l l one should perform now integration over Cooper 
modes in (I449l l. Observe, however, that in contrast to (1445b . where we had product of either two re- 
tarded or two advanced Cooperon fields, which restricted integration over one of the time variables, 
in the case of MT contribution ( |449b . we end up with the product between one retarded and one 
advanced Cooperon and the time integration running over the entire range t > t'. Precisely, this dif- 
ference between (1445b and ( |449b makes contribution Ssc to be local, while Smt non-local. Finally, 
in each of the Cooperon fields c, c, ( |449b . one keeps only contribution with the classical component 
of the order parameter and recovers Smj in the form given by (1345b . 

The remaining density of states part of the effective action Sdos emerges, similarly to Smt, from 
Tr{([SA3<;, Qjc])^}- This time one of the Qjc matrices is kept at the saddle point A, while another is 
expanded to the second order in 'W: 



5dos[A, A, <D] = ^Tr{c;(r)3^°°?c„.(r) + c,V(r)3^°°?c,-,"(r)} , (451) 



where 



A:;,(r,f)[A5^(r,f)-A^3^(r,f")]F,-,. -i- J dt"'Al^(rj)F,^rA''^(rj"')Fr'-," 

(452) 

It is important to emphasize here, that as compared to ( |446b and ( 14501 ), when deriving ?\f'^°^ it is 
not sufficient to take the approximate form of the distribution function ( |447b , but rather one should 
keep full F,. In what follows, we deal with the part of the action ( 14511 ) having one classical and 
one quantum components of the vector potential. The other one, having two quantum fields can 
be restored via FDT. To this end, we substitute Cooperon generators in the form ( 1448b into the 
action ( 1451b . We keep only classical components of Ax (the quantum one produce insignificant 
contributions) and account for an additional factor of 2 due to identical contributions from c and c 
Cooperons. Changing time integration variables t - t" - t and t + t" = Irj, one finds 

5dos[A, a, O] = ine^D Tr [A;^(r, 77 + r/2)[A5^(r, 77 + r/2) - A'-^j^ir, 77 - t/2)]F, 



X 9il^ + T/2 - t')0{t' - 77 + T/2)A*3^' Jr, ^^H^— L J J, 



(453) 



Note that owing to the step functions, integration over t' is restricted to be in the range 77 -t- r/2 > 
f > T] - r/2. Since F^ is a rapidly decreasing function of its argument, the main contribution to the 
r integral comes from the range r ~ l/T <sc 77. Keeping this in mind, one makes use of the follow- 
ing approximations: A^(r, 77 + r/2)[Aij[.(r, 77 + r/2) - Aij[.(r, 77 - r/2)] ^ rA^(r, 77)5,,A|^(r, 77) and 

^3c' '^^^1'' ) '^'^2 '' ) ~ l^ac^*"' ■^'li'^h aflows to integrate over t' expHcidy / dt'e{i] + 

r/2 - t')6(t' -ri + r/2) = t9{t). Using fermionic distribution function (1447b and collecting all factors, 
we find 

5dos[A,A,<I)] =7re2vDrTr[A^(r,f)5,A^3^(r,f)|A5^(r,0l'l (454) 

L J Jo smh(7rrr) 

where we set rj ^ t. Performing remaining integration over r and restoring ^dos ~ ^x'^oc FDT, 
we arrive at ^dos in the form given by ( 13491 ). Additional details of the derivation of the effective 
action ( 13391 ) can be found in 115411 . 
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